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Abstract 



Quantum information theory is a very young area of research offering a lot of challenging 
open questions to be tackled by ambitious upcoming physicists. One such problem is addressed 
in this thesis. Recently several protocols have emerged which exploit both continuous variables 
and discrete variables. On the one hand, outperforming many of the established pure continuous 
variable or discrete variable schemes, these hybrid approaches offer new opportunities. However, 
on the other hand, they also lead to new, intricate, as yet uninvestigated, phenomena. 

An important ingredient of several of these hybrid protocols is a new kind of entanglement: 
The hybrid entanglement between continuous variable and discrete variable quantum systems, 
which is studied in detail in this work. An exhaustive analysis of this kind of entanglement is 
performed, where the focus is on bipartite entanglement. Nevertheless, also issues regarding 
multipartite hybrid entanglement are briefly discussed. The quintessence of this thesis is 
a new classification scheme which distinguishes between effective discrete variable hybrid 
entanglement and so-called true hybrid entanglement. However, along the way, also other 
questions are addressed, which have emerged during the studies. For example, entanglement 
witnessing is discussed not only for hybrid entangled states, but also for fully continuous 
variable two-mode Schrodinger cat states. Furthermore, subtleties regarding entanglement 
witnessing in a certain kind of mixed states are examined. Not only theoretical classification 
and analysis of hybrid entangled states are discussed, but also their generation is presented 
and a few applications are demonstrated. 
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1. Introduction 



What is quantum information theory (QIT) about? 

QIT aims for the combination of information theory with quantum mechanics. Exploiting 
concepts like the superposition principle yields new striking effects, such as entanglement, 
unknown from classical information theory. Going from theory to the experimental regime, 
utilizing these new phenomena, QIT may lead to real-world applications for communication 
and computation. On the one hand, quantum communication deals with the transmission of 
quantum states, that is to say the transmission of information, since information is encoded 
in quantum states in QIT. On the other hand, quantum computation is about the quantum 
mechanical processing of information. Quantum computers might outperform their classical 
analogues concerning certain tasks, such as the factorization of large numbers as well as database 
searches. 1 While the construction of quantum computers is quite complicated and not yet 
realized on large scales, there are already emerging first commercial quantum communication 
systems, which are used for quantum key distribution (QKD). 

It turned out that the impossibility of cloning quantum states as well as the highly intriguing 
and non-classical phenomenon of entanglement are at the heart of QIT and many applications, 
it might bring up. Especially entanglement is a very surprising and non-intuitive effect, which 
was first considered by Einstein, Podolski and Rosen in 1935 in a famous paper, presenting the 
so-called EPR-paradox [EPR35]. Based on these findings, Einstein actually considered quantum 
mechanics being incomplete and regarded entanglement in correspondence with Max Born 
famously as "spukhafte Fernwirkung" or "spooky action at a distance" [Ein71]. During the same 
time, also Schrodinger considered entangled quantum states, as he pointed out the paradoxal 
effects of the quantum mechanical superposition principle, when applied on macroscopic 
objects. 2 In 1964 Bell considered quantum measurements on an entangled two particle system 
and could show that entanglement actually leads to non-local quantum correlations which 
cannot be explained classically [Bel64]. This corresponds to today's image of entanglement 
as a non-local quantum correlation between two or more particles, or, generally speaking, 
between quantum systems held by different parties. Indeed there are still a lot of interesting 
and unanswered questions to be tackled, especially concerning entanglement. In this context 

1 See Shor's algorithm for integer factorization [Sho97], and Grover's algorithm for searching unsorted databases 
[Gro96]. 

2 He considered a cat entangled with a two-level atom, for which the laws of quantum mechanics are valid [Sch35]. 
Due to coupling the atom with a poisoning mechanism, depending on the state of the atom, the cat is dead or 
alive. Since the atom can be in a coherent superposition of ground and excited state, the cat is necessarily in a 
superposition state as well. However, a cat which is somehow dead and alive is contrary to any experience from 
daily life. As shown later, this problem can be settled by considering decoherence effects immediately destroying 
any coherence in macroscopic superpositions. 
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1. Introduction 



also note that QIT is actually a very young area of research, being considered as an own branch 
of physics for only about 20 years. 

When searching through publications regarding QIT, one can find that most ideas and 
schemes are formulated either for quantum systems supported by an infinite-dimensional 
Hilbert space (continuous variables or simply CV), or for systems living in a finite-dimensional 
space (discrete variables or simply DV). Most tasks can be accomplished in both settings, however, 
either regime has its characterizing advantages and disadvantages. For example DV experiments 
are nearly always conditional and hence not very efficient, while their fidelities are quite high. 
On the contrary, CV schemes combine unconditional operation and high efficiencies with lower 
fidelities. Recently, so-called hybrid approaches emerged trying to combine the benefits from 
both regimes. However, these hybrid approaches lead to as yet uninvestigated phenomena, 
such as hybrid entanglement (HE), which is the focus of this diploma thesis. 

The thesis is organized as follows. Chapter 2 introduces QIT and provides fundamental con- 
cepts, necessary for tackling open questions regarding HE and further discussions. It especially 
sets out the basics of entanglement theory. An overview over optical hybrid approaches to 
quantum information is presented in chapter 3, which motivates the following chapter. The 
main part of the diploma thesis is chapter 4 dealing with HE. An introduction and a classifica- 
tion scheme are given in the first section, while the following sections work through different 
kinds of bipartite HE. In particularly, the distinction between hybrid entangled states supported 
by finite-dimensional subspaces and so-called truly hybrid entangled states not being describable 
in finite-dimensional spaces is pointed out and several exemplary states are investigated. The 
final section of the chapter provides some analysis regarding multipartite HE. In chapter 5 the 
generation of HE states is briefly discussed. Furthermore some applications of HE are presented. 
Finally, coming to a conclusion, a brief summary as well as an outlook are given in chapter 
6. The appendices provide additional material. Appendix A summarizes the abbreviations 
and notations used in this thesis, and Appendix B consists of proofs of theorems, ancillary 
calculations, and provides a summary of useful formulas. 
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In this chapter, the reader is provided with an introduction to QIT and its mathematical concepts, 
which are necessary to understand the main chapter of the thesis. First, the Hilbert space 
and the quantum states described therein are introduced from a mathematical point of view. 
Then it is briefly shown, how to describe DV and CV quantum systems and how to deal with 
quantum operations. As the focus of this diploma thesis is on entanglement theory, the concept 
of entanglement is explained and entanglement quantification as well as entanglement detection 
are discussed. 



2.1. Quantum States and Hilbert Spaces 

When physicists first investigated quantum phenomena, they had to think about how to actually 
describe quantum states. There are basically two approaches. On the one hand complex 
wavefunctions can be used, which are typically functions of position or momentum and give 
probability amplitudes. Their moduli squared yield probability densities. Therefore, these 
wavefunctions are easy to interpret and they will be useful if the probability distribution of 
position, momentum, or some other variable of a particle is of interest. However, wavefunctions 
are typically quite complicated and if more subtle properties of the quantum system, such as 
entanglement, are to be examined, it will probably be cumbersome to work with them. Hence, 
on the other hand, there is the so-called Dirac notation, which employs more abstract vectors in 
a Hilbert space to describe quantum states. 

So, what is a Hilbert space? A Hilbert space over the field of the complex numbers H(C) is a 
complete vector space on which an inner product is defined. The inner product allows length 
and angle to be measured and a norm to be defined. Completeness is defined in the following 
way: 

Definition 2.1. A metric space M with a defined norm is called complete if and only if every Cauchy 
sequence in M converges in M: 

M complete: (a n ) net ^ Cauchy sequence in M => lim a„ E M. (2.1) 

Definition 2.2. (a„) neK is called Cauchy sequence if and only if 

Ve>03N: \\a m - a n \\ < e Vm,n > N. (2.2) 
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More intuitively speaking, a space is complete if there are no points "missing" from it. For 
example Q is not complete, since the irrationals are missing, while IR is complete. 

Hilbert spaces may have any dimension, even infinite. This becomes clear, when the physical 
quantum states to be described by the Hilbert space vectors are considered. For example, a 
general spin system may have any finite number of spin states, while a quantum state described 
by the continuous variables position and momentum lives in an infinite-dimensional Hilbert 
space. Let the elements of the Hilbert space %{€) be denoted by so-called "ket"-vectors \tp). 
Then the according dual space V.'(C) contains the "bra"-vector elements (ip\ : H(C) — > C. 
Therefore, a map into complex numbers between a vector (ip\ and a vector \<p) can be defined 
as (ip\cp) E C. 

Consider a continuous basis of a Hilbert space of infinite uncountable dimension {\x} : x,y E 
R, (x\y) = 5(x — y)}. 1 The completeness relation reads: 

oo 

j dx\x) (x\ =1. (2.3) 

— OO 

So, writing 

oo oo 

\tp) = J dx(x\ip)\x)= J dxip(x)\x) (2.4) 

— oo — oo 

yields the connection to the previous wavefunction approach, explicitly (x\ip) = tp(x). The 
continuous {|x)}-basis may for example be the position or the momentum basis and ip(x) the 
corresponding wavefunction. If the basis of the Hilbert space is countable or even finite, we 
have the basis {|A:} : k, I E No, (k\l) = S k i} and the integral becomes an infinite or finite sum, 
but the argumentation is analogue. The wavefunction, in this case discrete, is just written as 
(k\ip) = ip k . 

Again, using the completeness relation, we obtain 

CO 

= J dxt(x)<p(x), (2.5) 

— OO 

or in the countable or even finite-dimensional case 

k 

Hence, every physical Hilbert space is actually the space of square-integrable L 2 (uncountable 
dimension) or square-summable I 2 (countable dimension) functions. Fortunately, I 2 is its 
own dual space, and hence, {\n} : n E N} being a basis, for any state \ip) — YLn tyn \n) the 
dual is given by (ip\ = ip* (n\ (the analogue is valid for L 2 and states in Hilbert spaces of 

Continuous variable and discrete variable states, as well as infinite-dimensional and finite-dimensional bases are 
discussed in more detail in section 2.2. 
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uncountable dimension). Therefore, the inner product between a vector \ip) and a vector \<f>) is 
given by (ip\<p), which is sometimes also called the "overlap" between \ip) and \<p). Using this 
definition of the inner product, the norm is defined by 



lll*>ll = V^>- ( Z7 ) 

In conclusion, every quantum system (in a pure state) can be described by its normalized 
({ip\ip) = 1) "ket"-vector \ip) in a Hilbert space V,^. Using a discrete orthonormal basis {\k} : 
k, I G N , (k\l) = S k i}, it can be written 

= E = E ^k\k), (2.8) 

k=0 k=0 

where xpk G C is the discrete wavef unction in the { | k)} -basis. If a continuous basis is chosen 
{\x) : x,y G R, (x\y) = S(x — y)}, every quantum state can be written as 

CO oo 

\ip) = J dx\x) (x\ip) = J dxTp(x)\x), (2.9) 

— oo — CO 

with the continuous wavefunction ip(x) G C, as previously already shown. 

In the following, {\k} : k,l G No, (k\l) = 4/} W1 H be abbreviated by {\k}}, and {\x) : x,y G 
M, (x|y) = ^(x — y)} analogously by for the sake of simplicity. 

An expectation value of a Hermitian operator A of a pure quantum state is given by 

(A) = (lp\A\lp) = a mPm, (2.10) 
m 

where the spectral decomposition for the operator A, 

A = J2 a m \m) (m\, (2.11) 

m 

was used [NC00]. a m are real eigenvalues, {|m)} is the corresponding eigenbasis, and p m = 
| (m\tp) | 2 are the probabilities for obtaining the result a m G IR. Note that Hermitian operators 
are measurable observables. Analogously, for a continuous basis the spectral decomposition is 
given as 

oo 

A — J dxa x \x) (x\ (2.12) 

— CO 

and the rest follows accordingly. 

If quantum mechanics is extended to quantum statistical mechanics, quantum systems may 
be statistical mixtures of several individual quantum states. Then the whole system is said to be 
in a mixed state and is no longer described by a Hilbert space vector, but by a so-called density 
operator p. A density operator is an operator, exhibiting the following properties: 
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• Linear (A linear <=> V \ip) , \cp) e H ; x,y 6 C : A(x |tf) + y |tf>)) = xA \ip) + yA \<p)) 

• Bounded (A bounded <=> 3 M > Osuch that V \xp) e H : < M (i/?|t/>)) 

• Positive semidefinite (A positive semidefinite <=> V \tp) E H : (tp\A\ip) > 0) =>• also 
Hermitian (A Hermitian A = A + ) 

• Trace-class (A trace-class <=> (k\A\k) < oo and independent of basis {\k}}) with trace 
tr[/3] = 1 

Note that the definitions of linearity, boundedness, positive semidefiniteness, and trace-class 
have been given in the context of the relevant Hilbert spaces with norm || \xp) || = ^/ (ip\ip). 
Any density operator can be written as a convex combination of rank-1 pure state projectors: 

P = EpihP«Xfcl; p.- >ovi, £ p,- = i. (2.13) 

This decomposition is what actually distinguishes pure states from so-called mixed states. 
While pure states cannot be written as a convex combination of two or more projectors, mixed 
states always contain two or more pure state projectors in their convex combination. 2 Hence, 
the density operator of a pure state \tp) is given by pp Ure = and thus: 

• For pure states: tr[p 2 ] = tr[p] = 1. 

• For mixed states: tr[p 2 } < tr[p] = 1. 

Furthermore, it is important to mention that there is no unique pure state decomposition for a 
given density operator. 

Expectation values translate into the density operator formalism simply by 

(A) =£ Pi(ipi\A\iPi) =tr\pA], (2.14) 

i 

which can be easily seen by inserting YLk W) (^1 = ^ ror an arbitrary orthonormal basis {\k}} 
into the term (tpi\A\ipi) and some rearrangements. 



2.2. Discrete Variable vs. Continuous Variable Quantum States 

So far, quantum states are defined as abstract density operators. However, these operators are 
difficult to handle, since for proper calculations, practical objects such as matrices are required. 
Therefore, this section discusses how to actually handle quantum states when working with 
them, and how to experimentally implement them. Two cases have to be distinguished: Either 
the quantum system is finite- or infinite-dimensional. 

2 Pure states are therefore extrema in the convex set of states. 
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2.2.1. Discrete Variable Quantum States 

First, we focus on the finite-dimensional case. The finite-dimensional Hilbert space is 
spanned by a discrete and finite set of orthonormal basis vectors { \n) : n — 0, . . . ,d — 1, (n \m) = 
$nm}- This is the reason why these quantum states are called discrete variable (DV) quantum 
states. A pure state of this kind looks like 

\if>) = Co |o)+ Cl |l) + ... + c d _ 1 |d-l) , £|c,-| 2 = l. (2.15) 

i=0 

For two dimensions this becomes 

|t/>)=c |0)+ Ci |1), (2.16) 

which is also called a qubit in analogy to the classical bit (cbit). This is actually the point, when 
it becomes clear, why it is interesting to do quantum information and why to do research 
on quantum computers: In contrast to the cbit, the qubit can be in a superposition of the 
zero- and the one-state. Therefore, a qubit contains more information than a classical bit 
and in one computation step more information may be processed. For certain tasks quantum 
computer algorithms, exploiting the quantum mechanical superposition principle, outperform 
their classical counterparts. Examples are Peter W. Shor's algorithm for integer factorization 
[Sho97], and Lov K. Grover's algorithm for searching unsorted databases [Gro96]. Shor's 
algorithm factorizes large numbers in polynomial time, while the best classical algorithm only 
achieves sub-exponential time [Pom96]. Its experimental demonstration was accomplished 
for N = 15 [VSB + 01, LBYP07, LWL+07]. Likewise Grover's algorithm provides a quadratic 
speedup, requiring 0(N 1 ^ 2 ) time for a database with N items, while classically this problem 
would only be solvable in linear time. However, quantum computation outperforms classical 
computation only for some special tasks and furthermore a large-scale quantum computer 
could not be realized yet. 

Back to the description of DV quantum states: Motivated by the term for a 2-dimensional 
quantum state, qubit, a 3-dimensional state is called qutrit and a d-dimensional, audit. For such 
a general qudit in a mixed state p consider p nm = (n\p\m) G C. This defines a proper d x d 
matrix containing complex numbers, which is called the density matrix of the quantum state p? 
These complex-valued density matrices are extremely suited for describing DV quantum states 
in a manageable way. Several important quantities, such as entropy or, in the multipartite case, 
entanglement measures, can be calculated from them. 

To illustrate the concepts introduced so far, consider the mixed qubit state 

^ = ^(|0)(0| + |l)(l|+ce^|l)(0|+ce-^|0)(l|), c < 1. (2.17) 



3 Often the terms "density operator" and "density matrix" are used interchangeably and mostly it is clear from the 
context whether the actual operator or the matrix is meant. However, they are actually not the same. 
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It possesses the density matrix 

1 / 1 ce~'^\ 
Pnm = - 2 ^ 1 J- (2.18) 

A pure state decomposition is given by 

P = ^ I'M (tf+l + ^ I?-) <*-| , |?±) = ^ (|0> ± e* |1>) . (2.19) 

How can DV quantum states be experimentally realized? More precisely, how can logical 
qubits or qudits be experimentally encoded? An apparent example is simply to use the ground 
and excited states of a 2-level atom, \0) L = \g) and |1) L = \e), to implement a qubit. Furthermore 
there are a lot of quantum optical encoding techniques. For example, in a finite-dimensional 
subspace of the infinite-dimensional Fock space, photons in a single mode can be utilized 
{single-rail encoding). For a qudit this becomes {|0) L = \0} Fock , |1) L = |1)f oc ; c , ... ,\d— 1) L = 
\d — l)f oc j.}- 4 However, it is more practical to use encodings, where each state holds the same 
number of photons. In the multiple-rail encoding a qudit is represented by a single photon in d 
different modes: {|0) L = |10...0), |1) L = |010...0), ... ,\d-l) L = |0...01)}. Unfortunately, 
this kind of encoding is not scalable, if several qubits or qudits are used. A compromise is 
offered by the dual-rail encoding which uses only two modes and a constant number of photons 
distributed between these two modes: {|0) L = \0 r d-l) , \1) L = \l,d - 2) , . . . , \d - 2) L = 
\d — 2, 1) , \d — l)r = \d — 1,0)}. These modes can be realized as spatial modes or as modes of 
orthogonal polarization. 

Additionally there are optical encodings utilizing coherent states instead of single photons. A 
logical qubit can be encoded with even and odd Schrodinger-cat states (also called coherent-state 
superpositions, CSSs). 

• Even cat state: \tp+) = + 

• Odd cat state: \tp-) - -^=(|«) - |~«))- 

J\f± are normalization constants. The names even and odd rely on the fact that the even cat 
state contains only even photon number states, while the odd cat state contains only odd 
photon number states. Since (n\m) = S nm for Fock states \n) and |m), the overlap between 
the even and odd cat state is zero: (ip+\ip-) = 0. Hence, a qubit can be represented as 
{\0) L = \ip+) t \l) L = \iP-)}\LBHO8\. 



2.2.2. Continuous Variable Quantum States 

This thesis focusses on QIT using quantum optical systems. It is straightforward to show that 
the quantized electromagnetic field is described by a set of quantum harmonic oscillators. The 

4 Fock states and coherent states are introduced in more detail in the next subsection, when the quantum harmonic 
oscillator is discussed. 
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Hamiltonian for a single mode is (unit mass) 

H =^(p 2 +cv 2 x 2 ^ =hcv(a f a + ^y (2.20) 

Here, £ and p, originally the oscillator's position and momentum operators, are now the 
so-called x- and p-quadrature operators of the field, a and a + are the annihilation and creation 
operators of the electromagnetic field, also called ladder operators or sometimes mode operators, 
which destroy or create excitations of the field, the photons. They are given by 

1 

{cox + ip), (2.21) 



V2tuv 

These operators define the so-called Fock space of Fock states or photon number states {\n) : n e 
No, (n\m) = S nm }. Fock states are eigenstates of the number operator h = a f a. 

a \n) = i/n \n - 1), (2.23) 
a f \n) = v / n~+T|n + l), (2.24) 
h \n) — n \n) , (2.25) 

where \n) denotes a state which contains n photons, while the energy of a single photon is 
given by hco. Therefore, the form of the Hamiltonian becomes comprehensible: hhco is simply 
the energy of n photons and additionally there is the vacuum energy \hco. 



However, also multimode systems can be considered. Photons in different modes stand for 
photons of different "fashion", what makes them distinguishable, in contrast to photons in 
the same mode, which are indistinguishable. These modes might correspond for example to 
different energies, different polarizations, or to spatial modes. For two modes and k the above 
defined operators satisfy the commutation relations 

[xj, p k ] = iUp (2.26) 

[*j>*k] = \Pj,Pk]=0> ( 2 - 27 ) 
=S jk , (2.28) 
[a jr a k } = [a/, 4] = 0. (2.29) 

Equations (2.26) and (2.27) are the well-known canonical commutation relations from general 
quantum mechanics. Equations (2.28) and (2.29), the bosonic commutation relations, simply 
follow from the definitions of the ladder operators. 
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Rewrite equations (2.21) and (2.22) into 

&=\J^{a + a*), (2.30) 

p = -i\j^{a-a'). (2.31) 

Hence, up to normalization factors, the x- and p-quadratures are just the real and imaginary 

parts of the annihilation operator. To get rid of these factors, define new dimensionless 
quadratures 

x' = Mx= -Ua + fl 1 "), (2.32) 
V « V2 

P' = {lP = ^-n (2.33) 



which obey the commutation relation 



[Xj,p k ]=i6j k . (2.34) 



Hence, this definition of dimensionless position and momentum operators corresponds to 
setting h — 1. Throughout this thesis, the dimensionless quadratures x.' and p' will be used and 
in the following the prime will be omitted. So, x. and p always stand for a pair of conjugate 
dimensionless quadratures obeying [x, p] = i. 



The eigenstates corresponding to the x- and p-quadratures are the position and momentum 
eigenstates \x) and \p): 

X\x)pcs. =*\x)pos. . P\P)mom. = P\P)mom. ( Z35 ) 

They are orthogonal, 

(xi\x 2 ) = &{xi -x 2 ) , (Pi\P2) =fi(pi-p2), (2-36) 

and complete, 

jdx\x)(x\=l, J dp\p) (p\ =1. (2.37) 

— oo — oo 

Hence, they form bases. However, since they are not normalizable, \x) and \p) are unphysical 
states. Nevertheless, they are very important and useful to calculate position (x\xp) = ip(x) and 
momentum wavefunctions {p\ip) = ip(p), which are indeed well defined. The relation between 
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the position and the momentum bases is given by a Fourier transformation [Leo97]. 

l*W = i / A V*- lxv \v) mom , (2-38) 



— CO 
OO 




(2.39) 



— oo 



The Hilbert space of the quantized harmonic oscillator is an infinite-dimensional "Hoo- States 
living in this Hilbert space are called qumodes. They can be either represented in the infinite, 
but countable Fock basis, or in the continuous x- or p-basis. Note that by qumode only states 
are denoted which are actually infinite-dimensional. Of course, a single Fock state or a single 
position eigenstate is also supported by the infinite-dimensional Hilbert space Hoo, but they 
only use a finite-dimensional subspace and could hence be characterized as DV qudits. 

One of the most prominent examples of qumode states are probably the coherent states \a), 
introduced by Roy Glauber in 1963 [Gla63]. 

U|2 00 a « 

\a) = e -^r V -—= \n) , a e C. (2.40) 
n=o vn! 

Coherent states are eigenstates of the annihilation operator a\m) = a \ a). They are particularly 
important because their field dynamics most closely resemble classical sinusoidal waves, such 
as continuous laser waves. Furthermore, coherent states form an overcomplete basis of the 
Hilbert space, which is another valuable property, 

— [ d 2 a la) (cc\ - 1 , d 2 a = d3m(a)d$He(a). (2.41) 
Note, however, that the coherent state basis is not orthogonal: 

(01*) =e z^l 2 +M 2 - 2 ^). (2.42) 



As previously shown, DV states can always be expressed by density matrices containing only 
complex numbers. Unfortunately, these matrices are not practical for describing qumodes. Due 
to the infinite-dimensionality of the Hilbert space, the density matrix would be of infinite size. 
Instead, the so-called phase-space representations (also denoted as quasi-probability distributions) 
can be employed. They are functions of the real quadrature variables x and p, from which 
any expectation values of the quadrature operators x and p can be calculated. There are 
different kinds of phase-space representations, which allow the calculation of different types of 
expectation values [Leo97, SZ97]. The most famous one is the Wigner function, which is suited 
for calculating expectation values of symmetrically ordered operators [Wig32]. Furthermore, the 
Glauber-Sudarshan P -representation allows for expectation values of normally ordered operators, 
while the Husimi-Q distribution yields expectation values of antinormally ordered operators 
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[Sud63]. 5 Note that from these phase-space representations, the density operator can always be 
retrieved. Hence, phase-space representations fully describe quantum states and the quantum 
phase-space picture is equivalent to the density operator approach. Qumodes living in an 
infinite-dimensional Hilbert space, described by the continuous quadrature variables x and p 
using phase-space representations, are therefore also called continuous variable (CV) quantum 
states. 

As an illustration of the phase-space representations the Wigner function is briefly discussed 
now: 

oo 

W(x, p) := I f dy (x - y\p\x + y) e 2 W. (2.43) 



7T 

— oo 



Its Fourier transform is the characteristic function, which is defined as 

X(fr,fc):=tr[%>]/ (2.44) 

making use of the Weyl operator 

ty = e i(£iP-&*). (2.45) 

The expectation value of a symmetrically ordered operator S S y m {x n , p m ) can be calculated then 
as 

oo oo 

(S sym (x n ,p m )) = I J W(x,p)x n p m dxdp, (2.46) 

— oo — oo 

where S sym denotes symmetrization. 

The most outstanding feature of the Wigner function is that it can become negative in contrast 
to classical probability functions (which is the reason why these representations are only 
called "quasf'-probability distributions.). Negativity of the Wigner function is an indication for 
nonclassicality of the state [KZ04]. Another important property the Wigner function offers is its 
marginal distributions. Integrating one variable out yields the probability distribution of the 
remaining one: 



p x {x) = J dpW(x,p), (2.47) 

— 00 
CO 

p v {p)= I dxW(x,p). (2.48) 

— 00 

Also note that the Wigner function is normalized, 

CO 

J W(x,p)dxdp = 1. (2.49) 



5 Also generalized, so-called s-parameterized, phase-space representation can be defined. For certain values of s the 
Wigner-, P- or Q-function can then be retrieved. 
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As an example, figure 2.1 shows the Wigner function of the cat state 




(2.50) 



with a = 6 and (p = ^ . Its negativity is clearly visible, which is due to the high nonclassicality 
of the state. The definitions of the Wigner function and the characteristic function are easily 




Figure 2.1.: Wigner function of the cat state of equation (2.50). 



generalizable to multimode states. For N modes the phase space becomes 2N-dimensional and 
the phase space representations become functions of 2N variables. 

On the one hand, these quasi-probability distributions offer a neat tool for describing CV 
quantum states. On the other hand, however, they cannot be used for the calculation of some 
important quantities such as the entanglement of the state. In contrast, in the DV setting, 
density matrices are well suited for entanglement quantification. Therefore, now an important 
subclass of CV states is introduced, the Gaussian quantum states, which are more well-behaved. 

Gaussian states are defined as those whose characteristic function is a Gaussian distribution 
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function, 6 i.e., for N modes, 

Yr ■ (F\ = Yr ■ (C\\ & -\f lllM+if InR fO 5D 

AGaussian\'p_) — AGaussian K u ) c - / 

£ 6 R 2N and D £ R 2N are vectors, 7 is a real symmetric 2N x 2N-matrix and the so-called 
symplectic matrix is defined as 

7n:=07, J--=(\ (2-52) 



;=i 



1 0, 



The vector D contains the first moments of the state, the displacements, while the matrix 
7 contains the covariances and is therefore called covariance matrix. Making use of a vector 
R = pi,---, x n, Pn) = ■ ■ ■ / &2n)/ D and 7 can be defined by 

di = \x\pki], (2.53) 
7, 7 = tripiRikj + %^)] - 2 tr^tr^]. (2.54) 

A consequence of the definition of Gaussian states is that they can be solely described by the 
displacements D and the covariance matrix 7. However, not every real symmetric 2N x 2N- 
matrix corresponds to a valid quantum state, since states also have to obey the Heisenberg 
uncertainty relation. In the covariance matrix formalism the latter translates into the inequality 

7 + iJn > 0. (2.55) 



Here the > stands for positive semidefiniteness (in the following, when talking about matrices 
or operators, > will always mean positive semidefiniteness.) [SMD94]. 

In contrast to Gaussian states, general non-Gaussian CV states also require higher moments, 
i.e. the whole infinite set of moments, for their description. This is the big disadvantage of 
non-Gaussian states compared to Gaussian ones. A convenient way of describing general CV 
states using practical C-numbers instead of operators is not known so far. 



2.3. Quantum Operations 

In this section the formalism behind quantum operations is introduced. Quantum operations 
are described by linear, completely positive (CP) maps. 

Definition 2.3. A map T : S(H\) — > S(^) mapping states S^Hi) of some Hilbert space T-L\ onto 
states S(H2) of some possibly different Hilbert space I-L2 which preserves the positivity of the density 
operator is called positive. Furthermore, if for larger Hilbert spaces also T ® 1„ is positive for all 
n £ N, T is called completely positive. 

6 Since the characteristic function and the Wigner function are related by a Fourier transform, Gaussian states also 
have Gaussian Wigner functions, as Gaussians are Fouriertransformed into Gaussians. 
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This requirement of complete positivity is clear, since every physical quantum operation 
should output a valid (and therefore positive) quantum state. 

Next it is to be distinguished between completely positive trace-preserving (CPTP) maps 
and completely positive trace-decreasing maps. Trace-preserving means that for a map T, 
tr[T(p)] = tr[jo] = 1, while trace-preserving denotes tr[T(p)] < tr[p] = 1 (trace-increasing 
quantum operations do not exist.). Linear CPTP maps are called quantum channels, while 
linear CP trace-decreasing maps always involve for measurement operations. 

While quantum channels correspond to a deterministic evolution of the state, measurement 
operations relate to a conditional and hence typically non-deterministic state evolution. How 
can quantum channels (CPTP maps) actually be described? 

Theorem 2.1 (Stinespring's dilation theorem [Sti55]). Let Y : S^H\) — > SCH2) be a CPTP map. 
Then there exists an ancilla Hilbert space K. of dimension N and a joint unitary evolution U onH® K. 
such that 

Y(p)=tr K [tt(p ® |0)(0|)Q + ], (2.56) 

for all p 6 S(%i). The ancilla space K, can be chosen such that dim[/C] = N < dim 2 [H]. 

Furthermore, M. Choi has shown that every CPTP map can be given in its operator-sum 
decomposition or Kraus decomposition: 

Theorem 2.2 (Kraus decomposition [Cho75, Kra71, Aud07]). Every CPTP map Y : S(H\) -> 
SCH2) can he written as 

N 

Y(p) = ^K,pKl (2.57) 

i=i 

for all p 6 S('Hi). The Kj : "H\ — >■ Hi are called Kraus operators and obey the completeness relation 

The proofs of the theorems can be found in the denoted references. Making use of a basis 
{\n)} of the ancilla Hilbert space /C, the Kraus operators can be obtained with the aid of the 
Stinespring unitary U of equation (2.56) via 

&n = (n\ti\0). (2.58) 

Now the measurement operations Ai m (p) are introduced [Aud07]. Such a measurement 
operation is defined as a complete (£ m MJ^M m = 1) set of operators {M m }. Each operator in 
the set corresponds to a possible measurement outcome m, which is obtained with probability 
p{m) = tr[Mj,M ml 6] < 1. The state after the measurement is 

MAP) = ^f^. (2.59) 

The normalization with pirn) has to be included to compensate the trace-decrease of the actual 
measurement operation M„ ;1 6Mj r 



18 



2. Basics 



To be even more general, each measurement outcome may not just correspond to a single 
operator M m , but to a whole set of Kraus operators {M m/i }. Then YL m ,i ; M m/ ; = 1 and the 
state after a measurement with result m is 

■Haw = — ^) — ' ( 2 - 6 °) 

which is obtained with probability p(m) = tr[X^- ; .M m i( o] < 1. 

An important class of quantum operations are those which can be implemented by using 
only linear optical elements. These unitary interactions are described by Hamiltonians of order 
< 2 in the mode operators and have a linear input-output relation with respect to the mode 
operators. Actually there are only 4 distinct such linear optical elements: 



Displacers: 



Phase Shifters: 



Beam Splitters: 



Squeezers: 



D(ct) =e aat -** fl , (2.61) 
D + (a)aD(a) =a + a. (2.62) 



t{(p) = e"f a \ (2.63) 
f + (<?>)flf(4>) =ae'f (2.64) 



%(6) = e 0(e»«&-e *ai«J) # ( 2. 6 5) 
B\{B)aB^{B) = cos 0a x + e l ^sin 6a 2 . (2.66) 



5(0, r) = ez r ( e ^ 2 - e ' ea+ ), (2.67) 
S\e,r)aS{6,r) = a cosh r - flV e sinh r. (2.68) 



These operations are also called Gaussian operations, since they map Gaussian states onto 
Gaussian states. Note that for the evaluation of the action of the operations on the mode 
operators, the Hadamard lemma has been exploited, which reads 

e x te -x = £ [xA (2 69) 

m=0 

for bounded operators X and Y, [X,Y] m = [X, [X, Y] m -\] and [X, Y] = Y. A proof for the 
theorem is presented in appendix B.l. 
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It is important to know that every pure Gaussian state can be created from the vacuum just 
with a displacement and a squeezing operation: 



For universal quantum computation it is not sufficient to work only with Gaussian states 
and operations. A non-Gaussian element is required for universality. In the previous section, 
it has been already shown that general non-Gaussian CV states are difficult to handle and 
unfortunately in contrast to Gaussian operations, non-Gaussian operations are just as well hard 
to implement efficiently. 

2.4. Entanglement 

The main part of this thesis deals with entanglement. So what actually is entanglement? 

Several interesting answers to this question are collected in [Bru02]. The most significant 
ones come from J. Bell, A. Peres and D. Mermin. Bell writes, "entanglement is a correlation 
that is stronger than any classical correlation". This leads to the fact that entangled states 
violate the well-known Bell inequalities [Bel64]. This violation has the consequence that 
quantum mechanics cannot be a realistic and local theory at once, in line with Mermin, who 
writes, "entanglement is a correlation that contradicts the theory of elements of reality". So 
entanglement cannot be explained or simulated classically and hence quantum correlations 
due to entanglement can be concluded to be stronger than any classical ones. Realism in the 
sense of physical theories means that measurements just read off predetermined properties, 
which are so-called elements of physical reality and also exist if they are not measured at all. 
Locality means that for spatially separated particles a measurement on one of the particles 
cannot instantaneously affect the other particle. Classically such actions propagate at most with 
the speed of light. Note that in the widely established Copenhagen interpretation of quantum 
mechanics, both locality and realism are abdicated. However, there exist alternative approaches 
to quantum mechanics, giving up only one of the two, such as Bohmian mechanics, which is 
realistic, using hidden variables, but also non-local [BH93]. 

Peres says, "entanglement is a trick that quantum magicians use to produce phenomena 
that cannot be imitated by classical magicians". An example of such a "trick" is quantum 
teleportation [BPM + 97, UJA+04], which has its origin in a paper by C. H. Bennett et al. from 
1993 [BBC + 93]. However, entanglement is not just the main resource for quantum teleportation 
but for virtually all applications in quantum information, such as long-distance quantum 
communication and quantum computation. 

To come to a conclusion entanglement can be understood as a strong non-classical quantum 
correlation, which is the basis for most applications in quantum information. Now entanglement 
is defined and discussed in a mathematical way. The focus will be on bipartite entanglement. 7 

7 A very comprehensive and accessible introduction into entangled systems can be found in the book by J. Audretsch 
[Aud07], which was one of the main sources for this chapter. 




(2.70) 
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Definition 2.4. An n-partite quantum state p l - n on the product Hilbert space W} (3 ... ® % n is called 
separable if and only if it can be written as 

p X - n = E ViPi ® ■ ■ ■ ® Pi . Pi > V i, Y,Pi = l - (2-71) 

i i 
Any state which is not separable is called entangled, inseparable or quantum correlated. 8 

Note that, if the convex combination in equation (2.71) consists of only one term, the state is 
in a fully uncorrelated product state p 1 ® . . . ® p n . For a convex combination containing several 
terms there are still classical correlations present [Aud07] . 

An important class of transformations regarding entanglement are the so-called LOCC 
operations (LOCC stands for Local Operations and Classical Communication). Since classical 
correlations can be created via LOCC, entanglement can also be defined as those correlations 
which cannot be produced using LOCC [PV06]. Such a definition is consistent with the 
mathematical definition 2.4, since a quantum state p can be generated from a resource of 
separable states via LOCC if and only if it is separable. It can be also shown that separable 
states can be created from entangled states using LOCC. More generally, there are states which 
can be obtained from entangled states via LOCC which are still entangled. It is reasonable to 
assume that these states are "less entangled". So, what does "less entangled" mean? This leads 
to the problem of entanglement quantification: How much entangled is a given quantum state? 

2.4.1. Entanglement Measures 

A bipartite Entanglement Measure is a function E(p) which quantifies the entanglement of a 
given quantum state. What properties should such a function offer [PV06, HHHH09]? 

1. For quantum states p, E(p) is a mapping from density matrices into positive real numbers 

p^E(p)eR+, (2.72) 

which is vanishing on separable states: 

^separable => E(p) = (2.73) 

Actually, E(p) is only required to be minimal for separable states. However, it is reasonable 
to set this constant to zero. 

2. Monotonicity under LOCC: Knowing that entanglement cannot be created via LOCC, it is 

8 It is worth mentioning that there also exists a classification which goes beyond this two sidedness, introducing the 
so-called quantum discord [OZ01]. Quantum discord is a quantum correlation which can be present in separable 
states but which cannot be simulated classically. It has an operational interpretation in state merging [CAB + 10]. 
However, this is not relevant for this thesis and the two-sided classification between separability and entanglement 
is sufficient. 
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clear that entanglement cannot increase under LOCC. So, for any LOCC operation A: 

E(A(p)) < E(p) (2.74) 

Unfortunately, the mathematical description of LOCC operations is very intricate [DHR02]. 
Instead so-called separable operations can be exploited, which are defined as 

A(p) =Y2A®Bi{p)A* ®Bj , J^A^Ai^BjBi =1®1. (2.75) 

i i 

Every LOCC operation can be cast in the form of separable operations. 

Actually, most entanglement measures also satisfy the stronger condition that they do not 
increase on average under LOCC, 

Y^PiE(Ci)<E(p), (2.76) 

i 

where {p;, g,} denotes the ensemble obtained from the state p via LOCC. Mostly it is 
easier to prove this stronger condition than the actually required condition of equation 
2.74. 

These properties are the essential ones a function has to satisfy in order to quantify entanglement. 
However, there are some more properties, which can be helpful: 

3. Normalization: For maximally entangled states 

E(|?+»=logd. (2.77) 

Maximally entangled states only exist in bipartite systems out of two d-dimensional (d 
finite) subsystems in Hilbert spaces of the form T-L^ ® H^. They are always pure and can 
be written as 

= * (|0,0) + |l,l)+...+ |d-l,d-l)). (2.78) 

The notion of maximal entanglement comes from the fact that any state in (3 can 
be prepared from such a state with certainty via LOCC. 

4. For pure states, E(\xp) (tp\) reduces to the entropy of entanglement Eg, which is defined as 

E s Qip) := S(tr A [\ip) (xp\j) = S(tr B [|?) <</>!]), (2-79) 

where S denotes the von-Neumann entropy S(p) = —tr{p\ogp] and tr^ (trg) denotes partial 
tracing over subsystem A (B). 

It is important to know that in the literature sometimes a distinction between entanglement 
measures and entanglement monotones is made, but the definitions are not always consistent 
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regarding different references. Some authors use the terms interchangeably. As an example 
here the definitions of [PV06] are presented: 

Definition 2.5. A function E(p) is called entanglement monotone if and only if it satifies properties 
1, 2 (no increase on average under LOCC) and 3. 

Definition 2.6. A function E(p) is called entanglement measure if and only if it satifies properties 
1, 3 and 4 and does not increase under deterministic LOCC. 

In this thesis, when generally talking about measures and monotones, it is simply used the 
term "entanglement measures" instead of always explicitly mentioning both. 

There is a whole bunch of further properties which entanglement measures may satisfy. For 
example, there are convex entanglement measures and (fully) additive ones. Then there are the 
properties of (asymptotic) continuity and lockability. Since this is supposed to be a rather brief 
introduction to entanglement measures, these properties are not explained here. They are not 
exploited in this thesis anyway. Instead the reader is again referred to [PV06, HHHH09], which 
provide a very comprehensible introduction to entanglement measures. 

Conceptually three different kinds of entanglement measures can be distinguished. First 
there are the physically motivated operational entanglement measures. Then there are the norm- 
based and distance-based entanglement measures and finally axiomatic entanglement measures can be 
defined. An operational measure results for example from the answer to the question "Having 
a supply of n copies of a quantum state p, how many copies rn of a maximally entangled state 
can be created in the limit n — > oo?". Then the rate r is defined as the distillable entanglement Er). 
An example for distance-based measures are the relative entropies of entanglement E%. Based on 
the quantum relative entropy they compare the correlations in a given quantum state with the 
correlations of the closest state from a set of states X. For X, any set of states can be chosen, 
such as the separable states or the nondistillable states. Finally, an axiomatic measure is just a 
function which has been constructed in such a way that it obeys the necessary requirements of 
entanglement measures. However, some axiomatic and norm/ distance-based measures were 
given operational interpretations subsequently. 

Since the description of quantum states in the DV and in the CV regime proceeds quite 
differently, there are also different ways towards entanglement quantification. While for DV 
states density matrices can be employed, for CV states only covariance matrices can be used. 
Hence for non-Gaussian CV states general entanglement quantification is yet an unsolved 
problem. In the following subsections some relevant measures are presented. 

2.4.2. DV Measures 

One of the most important DV entanglement measures is the pure state measure entropy of 
entanglement E$, which has already been defined in equation (2.79). Operationally interpreted 
Eg specifies the maximal reversible rate r of transforming n copies of the pure state at hand via 
LOCC into rn copies of a maximally entangled state and vice versa for n->». However, for 
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mixed states, these transformations are not reversible anymore. This is the reason, why Eg only 
works for pure states and is no entanglement measure anymore for mixed states. 

When talking about DV pure states, the Schmidt decomposition is a very important tool 
[Aud07]. 

Theorem 2.3 (Schmidt decomposition). Let \ip) AB be a normalized pure bipartite state in the product 
Hilbert space 1-L A <8> with dim [H A ] = n and dim [7^,] = m. Then there exist orthonormal bases 
{\ui) A } and {\wi) B } such that 

'"(</>) 

\<P) AB = E V^iWi) K) B / «i > V i, £>; = 1, (2.80) 

i=l i 

with r(ip) < min{n,m}. Such a decomposition is called Schmidt decomposition of\xp) AB . 

The number r(xp) is called the Schmidt rank of \ip) AB and {a,} are its Schmidt coefficients. What 
can be seen from the form of the Schmidt decomposition is that the reduced density operators 
p A = \x^[p AB ] and p B = tr y\[p AB ] both have the same positive eigenvalues {a/}. Furthermore, 
{\u,i) } and {|if,) B } with suitably chosen phases are the orthonormal eigenvectors of p A and 
p B . They are uniquely determined up to a phase. A proof of the Theorem is provided in 
appendix B.2. Note that the Schmidt decomposition is essentially a restatement of the singular 
value decomposition in a different context. For the sake of completeness, a singular value 
decomposition is defined in the following way: 

Theorem 2.4 (Singular value decomposition). Let M be a n x m matrix over the field K which is 
either the field of the real or complex numbers. Then there exists a factorization of the form 

M = UAV\ (2.81) 

where 

• U is an x n unitary matrix over the field K, 

• V is am x m unitary matrix over the field K, 

• and A is an x m diagonal matrix with nonnegative real numbers on the diagonal. 

Such a decomposition is called singular value decomposition of M. The diagonal entries of A are 
called the singular values of M. 

The singular values correspond to the Schmidt coefficients of the previous theorem 2.3 and 
the unitary matrices U and V correspond to matrices which transform the state into its 

new basis. The proof of theorem 2.4 corresponds to the proof of the Schmidt decomposition, 
which is presented in appendix B.2. 

It can be shown that any pure DV bipartite state is in a separable product state if and only if 
it has Schmidt rank r — 1. Furthermore, for states \ip) AB and \(p) AB in *H A <g> TL^ (n < m) it can 
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be shown that \ip) can be LOCC-transformed into \(p) with unit probability if and only if 
the Schmidt coefficients {af } are majorized by {txf } when taking them in decreasing order, i.e. 
«!>...>«„ [Nie99]. 



Definition 2.7. {af : i — 1, . . . , n} is said to be majorized by {af : i = 1, . . . , n}, {af} -< {&f}, if 

E«f<E a f v fc = 1 ( 2 - 82 ) 

;=1 i=l 

If the latter is the case, \ip} AB can be regarded as more entangled than \(p} AB - This is also the 
reason why states with equally distributed Schmidt coefficients are maximally entangled. Their 
set of Schmidt coefficients is majorized by all other sets of Schmidt coefficients. Maximally 
entangled states in H A ®H B 2 with d\ < d^ look like 

l<> = 4 r (|0,0) + |l,l)+-.. + |di-lA-l>)- (2.83) 



However, a consequence of the majorization criterion is that there exist also incomparable states, 
of which neither can be seen as more or less entangled than the other. Other measures have to 
be employed for such states [PV06]. 

A very important mixed state entanglement measure is the entanglement of formation Ef 
[PV06]. It is defined as 

E F (p) := inf{J>E s (hfc) : P = Ep» l#> }• ( Z84 ) 

i i 

It represents the minimal average entanglement over all pure state decompositions of p, where 
the entropy of entanglement is employed as the pure state measure. Unfortunately, due to the 
variational problem, the entanglement of formation is extremely difficult to calculate. However, 
it is a very important measure for two reasons. On the one hand, some open questions regarding 
Ef have tight connections to major open questions in quantum information. Explicitly full 
additivity of the entanglement of formation is equivalent to the additivity of the classical 
communication capacity of quantum channels [Sho04]. However, recently it has been shown by 
providing a counterexample that these quantities are actually not additive [Has09]. 

On the other hand, Wootters could derive a closed analytical solution for the case of bipartite 
qubit states [Woo98, WooOl]. So, for p a two-qubit state in U A ® U\, 

Ef(W= .(!±^E), ,2.85, 



with 



s(x) := — xlog 2 x — (1 — x) log 2 (l — x). 



(2.86) 
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C is the widely-used Concurrence, which is defined as 

C(p) := max{0, ^ - ^ - - (2.87) 
where are the eigenvalues of pp in decreasing order, p is given by 

p = (fry ® fry)?* (fry (8) fry), (2.88) 

where p* denotes the elementwise complex conjugate of p and fry is the Pauli Y matrix. Since 
the entanglement of formation and the concurrence are monotonically related, many authors 
prefer to simply quantify entanglement using the concurrence instead of the entanglement 
of formation itself. For higher dimensions this connection breaks down, as the concurrence 
is not even properly defined anymore. Note, however, that there actually are approaches for 
generalizing the concurrence to higher dimensions [Gou05]. 

Another very important entanglement monotone is the logarithmic negativity [VW02]. To 
define £jv, first the notion of partial transposition has to be introduced. Consider a bipartite DV 
quantum state in local orthonormal bases: 

\<P) AB = L c ^\i) A (l\®\k) B (l\- (2.89) 

i,jjc,l 

Then the partial transposition Tg with respect to system B is defined as 

(p AB ) TB = (l A ®T B )[p AB ] 

= (lA®T B )C£c ijkl \i) A (i\®\k) B (l\) 

iW ( 2 - 9 °) 

= L c w\i) A (i\®\i) B (k\, 

where Tg denotes "normal" transposition of system B. Partial transposition Tg is a positive but 
not completely positive (PnCP) map. So if applied on an entangled quantum state, the output will 
not necessarily be positive semidefinite. However, if applied on a separable state, the output 
will actually be positive semidefinite, since the subsystems of separable states fully separate 
and hence the map acts as two individual maps 1 A and Tg on the state. These individual maps 
are both positive and hence the output will be positive semidefinite again [Per96]. Therefore, 
partial transposition yields an inseparability criterion, which is also called the Peres-Horodecki 
criterion [HHH96]. 

The concept can be generalized: All PnCP maps distinguish some entangled states from the 
separable and the other entangled states. If the output of a PnCP operation is not positive 
semidefinite, it can be concluded that the input was entangled. States which are not positive 
semidefinite anymore after partial transposition are called NPT (negative partial transpose) states, 
while states which are still positive semidefinite are called PPT (positive partial transpose) states. 
Note that PPT states may nevertheless contain entanglement - so-called bound entanglement. It 
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is just not possible to detect this entanglement via partial transposition. What makes the partial 
transposition particularly interesting is that it can be shown that "PPTness" coincides with 
nondistillability [HHH98]. Nevertheless, it is not clear yet whether all NPT states are distillable, 
or whether there exist nondistillable NPT states. 

The logarithmic negativity is a measure which attempts to quantify the negativity in the 
spectrum of the partial transpose. Therefore, it is not able to measure bound entanglement. It is 
defined as 

E N (p) :=log 2 ||/) r «|| =log 2 [l + £(|A I | -A,)], (2.91) 

i 

where | |X| | := tr[VX + X] is the trace norm and A, are the eigenvalues of p TB . It does not matter 
whether the partial transposition is performed with respect to system A or B. Sometimes also 
the related "normal" negativity J\f is used [VW02], which is defined as 

m) : = Mq=l = UMzM. (2 . 92) 

While J\T is convex, but not additive, is not convex, but fully additive. Since additivity is in 
the majority of cases more desired than convexity, the logarithmic negativity is the more often 
used monotone. The major advantage of the negativity quantities is that they can be calculated 
rather easily. It is just an eigenvalue problem to be solved. 

The last monotone, which is briefly presented in this subsection, is the global robustness of 
entanglement R g [PV06, VT99, Ste03, HN03]. 

R g (p) ■= inf{A : A > 0, a 6 QS such that (1 - X)p + Xcr 6 SEP}, (2.93) 

where QS denotes the set of all quantum states and SEP the set of the separable quantum 
states. It quantifies the minimal amount of arbitrary (global) noise a to be mixed in such that p 
becomes separable. However, also other robustness quantities can be defined by considering 
specific types of noise and drawing a for example from the set of separable states or PPT states. 

Robustness monotones can be sometimes calculated, often at least bounded nontrivially. 
Furthermore they find applications in proofs of theorems and similar argumentations, such as 
in this thesis in subsection 4.4.3. 

2.4.3. CV Measures 

Entanglement quantification in the CV regime is a much harder task than in DV. For general 
states the continuity property of some measures breaks down: It is straightforward to construct 
pure product states which have arbitrarily highly entangled states in their arbitrarily small 
neighborhood [PV06, ESP02]. This problem can be solved by applying an energy bound on 
the states. When only states with bounded mean energy, i.e. Sm = {p £ S : tr|j6J?] < M}, are 
considered (S denotes the set of all quantum states.), continuity can be recovered. However, 
there is still no recipe for entanglement quantification in CV even for these energy bounded 
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states. Nevertheless, for specific measures (strongly superadditive ones) 9 rough bounds may be 
calculated with the aid of Gaussian states and their entanglement. Therefore, the set of states is 
restrained even stronger and only the Gaussian states are considered in the following. Once 
more, even for this class of states entanglement quantification is an extraordinarily difficult task. 
However, there do exist at least a few measures which can be calculated. 

One such measure is the entropy of entanglement, which features a translation into the pure 
state Gaussian CV regime. For Alice and Bob holding and n$ harmonic oscillator systems 
in a Gaussian state, 10 the entropy of entanglement is 

Es^) = t^o g2 ti±l-ti-l logl '^l). (2.94) 

1 = 1 

Hi are the symplectic eigenvalues of Alice's reduced covariance matrix y A , which is the submatrix 
referring to Alice's system, of the overall covariance matrix y AB . On the one hand, these sym- 
plectic eigenvalues are the absolute values of the "normal" eigenvalues of iJ^j A . On the other 
hand, the general definition of symplectic eigenvalues corresponds to the Williamson normal 
form. When employing covariance matrices no Hermitian matrices and unitary transformations 
are used anymore as in the case of density matrices. Instead linear Gaussian transformations are 
described by so-called symplectic transformations. Due to the preservation of the canonical com- 
mutation relations, symplectic transformations S are defined by their action on the symplectic 
matrix J^: 

SJ N S i = J N - (2.95) 

This is verified by the real 2N x 2N matrices S, which form the real, symplectic group Sp(2n,R). 
Covariance matrices are then transformed according to SyS T . Williamson proved that for any 
covariance matrix 7 on N harmonic oscillators there exists a symplectic transformation S such 
that 

Then the diagonal elements ji{ are the symplectic eigenvalues of 7 and the whole set {fi,} is 
called the symplectic spectrum of 7 [Wil36, SMD94]. 

Another entanglement measure which can be also employed in the Gaussian CV case is 
the logarithmic negativity. In contrast to the entropy of entanglement it also works for mixed 
Gaussian states. Just as in the DV case, where it is the one exception of many measures which 
is relatively easy to calculate, in the CV case it is the only mixed state measure which in 
general can be actually calculated at all. Note that for pure states it differs from the entropy of 
entanglement and hence it rather ought to be called a monotone. If Alice and Bob are in the 



Strong superadditivity of an entanglement measure E means E(,of 2 B ) > E(pf B ) + E(pAB), where 1 and 2 refer to two 
pairs of entangled particles, held by systems A and B. See also [PV06]. 
°In quantum information it is common to give two system A and B specific names - Alice and Bob. Especially in 
quantum communication it is very popular to call sender and receiver Alice and Bob. 
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possession of N = n& + n$ quantum harmonic oscillators, the logarithmic negativity is 

N 

En(P) = - £log 2 [min{l,^}]. (2.97) 

{p-i} is the symplectic spectrum of the partially transposed state described by the covariance 
matrix 7 Fb . On the level of covariance matrices, the partial transposition corresponds to a time 
reversal. In a system descibed by the canonical quadratures x and p the time reversal operation 
can be expressed as x. — > x and p — > — p. 

Just as in the DV case partial transposition can be used for the derivation of a necessary 
condition for separability of CV states [SimOO]. Recalling the Heisenberg uncertainty relation 
on the level of covariance matrices (2.55), for every separable state 



Note that this is valid for all CV states, not just Gaussian ones. 

Finally, it has to be pointed out once more that entanglement quantification in the CV regime 
is much more difficult than in the DV regime. Only for Gaussian states, there exist some rare 
examples of measures which can be properly defined and possibly calculated. Entanglement 
quantification for general non-Gaussian CV states is yet an unsolved problem, even for energy 
bounded states. However, all physically relevant states are of course energy bounded. 

2.4.4. Entanglement Witnesses 

There are a lot of cases where entanglement quantification is too ambitious, for example, in 
the non-Gaussian CV regime or simply when matrices have to be evaluated which are too 
big for eigenvalue problems to be solved. Then it can be asked whether the entanglement, if 
not quantifiable, can at least be detected or witnessed. Furthermore, there are tasks where the 
amount of entanglement is actually rather secondary and it is only important that entanglement 
is present at all. For these purposes, so-called entanglement witnesses can be used [PV06]. 

Definition 2.8. A Hermitian operator W is called an entanglement witness if and only if: 



An entanglement witness therefore separates some entangled states from the separable states 
and the rest of the entangled states. Hence it can be interpreted as a hyperplane in state 
space as illustrated in figure 2.2. Arguing with this hyperplane interpretation, it becomes clear 
that a quantum state is separable if and only if it has nonnegative expectation values of all 
entanglement witnesses. It can be shown that for any entangled state an applicable witness 
detecting it can be constructed [HHHH09, Aud07]. 



r 



Xb + iJ N > o. 



(2.98) 




(2.99) 
(2.100) 



- 1 Pent 

si. tr[Wp ent ] < 0. 
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Figure 2.2.: An entanglement witness can be interpreted as a hyperplane in state space. All 
states to the right of the hyperplane, e.g. all separable states, have nonnegative 
expectation values of the witness operator and hence, it cannot be concluded 
whether they are entangled or not. However, the entangled states to the left of the 
hyperplane are detected by the witness [HHHH09]. 



An example of a witness operator for the T-Ld ® T-Ld case is the swap operator 

? = £ 10 01 ® l/> (A • (2-ioi) 

y=o 

It is straightforward to show that this operator yields positivity for all separable states. However, 
it can be also shown that it possesses negative eigenvalues, which qualifies it as an entanglement 
witness [HHHH09]. 

A very efficient way of entanglement detection has already been set out and explained: It 
is not provided by actual witness operators but instead by PnCP maps such as the partial 
transposition (see subsection 2.4.2). Due to the positivity of the actual map and the added l„'s, 
PnCP maps always leave separable states positive. However, their missing complete positivity 
enables them to detect some entangled states - the NPT states in the case of partial transposition. 

In the preceding subsection it has been pointed out that entanglement quantification for 
general CV states is still a hard problem. So, in this regime entanglement witnesses are especially 
desirable. Shchukin and Vogel provided a very elegant way for entanglement witnessing in this 
class of states [SV05a, MP06, MPHH09]. A slightly adapted version of their criteria gives rise 
to one of the central tools used in this thesis. Therefore, their result is discussed in more detail. 
However, their method also is not based on witness operators, but on so-called expectation value 
matrices (EVMs, also matrices of moments). For a set of operators {M/t} me EVM x is defined as 

Xij=(M]Mj), (2.102) 

where the observables are typically of the form {M^} = {^k ® S/ c }. Then it is straightforward 
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to see that for separable states, % also is separable. So when combining the EVM approach with 
a PnCP such as partial transposition, the separability of the state in question can be investigated 
based on the matrix of moments [GT09]. Inseparability, which may result in non-positivity of 
the mapped state, then also leads to non-positivity of the matrix of moments. Hence criteria 
based on subdeterminants of the EVM can be derived. However, this will become clear in the 
derivation of the Shchukin-Vogel inseparability criteria. 11 

Theorem 2.5 (Shchukin-Vogel inseparability criteria). A bipartite CV state p is NPT if and only 
if there exists a negative principal minor of M(p TB ), i.e. det[M(p lB ) r ] < for some r= {r\,...,r^) 
with r\<ri< ... < r N . 



The matrix of moments M is here defined as 

Mij(p) = M ili2hk>hhj3k {p) 



= (a a a a b b b b 



(2.103) 



With a convenient arrangement of the indices, i.e. combining i\, ii, h, to i and j\, ]i, 73, 74 to j, 
M can be explicitly written as 



M(p T ) 



( 1 


(a) 


<« + ) 






(« + fl) 




(a f b f ) 


(a) 




(«« + ) 


(atf) 


(b) 


(ab) 




(Pb) 


V : 









(2.104) 
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The process of the index arrangement is not crucial, however, it is explained in detail in the 
proof. Also note that the matrix of moments M is a matrix of infinite size. What actually are 
principal minors? 

Definition 2.9. Consider the submatrices M r of a square matrix M which are obtained by deleting all 
rows and columns apart from the the ones with the numbers r = (/%, . . .,rjj) ; T\ < ri < ... < r^. The 
determinants det[M r ] of these submatrices are called principal minors of M. 

Proof of Theorem 2.5. [SV05a] Recall that a Hermitian operator A is called nonnegative if 

ft\A\il>)=1x[A\ip)(ip\]>0V\ip). (2.105) 

The nonnegative operator \tp) (ip\ can be rewritten as / + / with f = \vac) (xp\. \vac) = 1 00) 
denotes the vacuum in the bipartite case. Any pure bipartite state | ip) can be expressed as 



\ip)=g f \00), 



(2.106) 



11 For a nice introduction and a detailed overview of entanglement detection, see Guhne and Toth's review article 
[GT09]. 
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for an appropriate operator function g = g(a, b). a and b denote the annihilation operators of 
the first and the second mode respectively. Hence 

/=|00)(00|g. (2.107) 

It is straightforward to show, see appendix B.3, that the two-mode vacuum density operator 
can be expressed in the normally ordered form 

|00) <00| =: e -" H - ph :, (2.108) 

where : — : denotes normal ordering. From equations (2.107) and (2.108), it can be seen that 
the normally ordered form of / exists. Hence, a Hermitian operator A is nonnegative if and 
only if for any operator / whose normally ordered form exists, 

</7) A =tr[A/7]>0. (2.109) 

Now the Peres-Horodecki criterion can be applied. For any separable state p, 

(Pf)?B=*\f B f*f\>0. (2-HO) 

which is satisfied for any operator f whose normally ordered form exists. Due to this existence, 
/ can be written as 



Then 



/= E c i02I3! yvw. (2.ni) 

h,'2,hM=0 



with the moments of the partial transpose 

M hhhk,hiihk = ( a a a a b b b b )p r B • ( 2 -H 3 ) 

The left hand side of inequality (2.112) is a quadratic expression in the coefficients c;w,j 3 j.. 
Sylvester's criterion states that such an inequality holds true for all c !l!2 / 3 / 4 if and only if all 
principal minors of the expression (2.112) are nonnegative [Joh90]. However, first a relation 
between the moments of the partially transposed and the moments of the original state has to 
be derived. This is straightforward: 

M hmkhhhk = ( a a a a b b b b y B 

= {a ¥2 ^a ¥l a il {b Vi b i3 b ¥3 b ii y) p (2.114) 
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Furthermore, it is convenient to establish an ordering of indices and number the ordered 
multi-indices by i and (71,72/ h, ji) by 7. For the inseparability criteria the ordering 

is irrelevant, however, an example is the following: Order the set of multi-indices {(ii,i2>h/H)} 
in such a way that for two multi-indices u = z'2, h, z'4) and v = ]!,]?>, ji) 

!\u\ < \v\ or 
f (2.115) 
\u\ — \v\ and u <' v. 

\u\ = i\ + ii + z'3 + z'4, and u <' v means that the first nonzero difference 73 — 13,74 — 24,71 — 
— h is positive. This is the ordering Shchukin and Vogel applied in their original paper 
and which results in a matrix of moments as in equation (2.104). 

With the relation between the moments of the partially transposed and the original state as 
well as with the multi-index ordering procedure the nonnegativity of all minors of the quadratic 
expression (2.112) can be expressed as 

det[M(/3 r B)], >0 V r_ = {r lr ...,r N ), r x < r 2 < ... < r N , (2.116) 

where the matrix of moments is defined by 

Mijip 1 *) = (^Va+V^W^ 4 )^, (2.117) 

with the just introduced multi-index ordering. Note that i and 7 have been exchanged in order 
to follow the convention that ordinary matrices are denoted with Mjj and not with Mp. 

Summing up, for every separable state, all principal minors are nonnegative, explicitly 
det[MdDs e B p)],- > V r. Furthermore a necessary and sufficient condition for the positivity of 
the partial transpose can be formulated: The partial transpose of a bipartite CV quantum state 
is nonnegative if and only if all principal minors det[(j6 FB )]r are greater or equal zero. This 
also yields the necessary and sufficient condition for NPTness, which has been stated in the 
theorem. □ 

The theorem provides a recipe for determining whether a CV state is NPT (and hence 
entangled) or not. However, it consists of an infinite set of inequalities to be examined. 
Therefore, the problem lies in the search for an appropriate principal minor. Nevertheless, 
for any NPT state such a minor exists. Since this infinite set of criteria is based on partial 
transposition it can only detect NPT entanglement. PPT entanglement (bound entanglement) 
remains hidden. It is also worth pointing out that all occuring moments can be measured 
experimentally [SV05b]. 

The SV criteria (in the following "Shchukin- Vogel inseparability criteria" will be abbreviated 
by "SV criteria"; furthermore "SV determinant" stands for a "principal minor of the matrix 
of moments by Shchukin and Vogel".) are especially outstanding as they provide not just 
sufficient but also necessary conditions for NPTness. There have been proposed several sufficient 
conditions for NPT entanglement before [SimOO, DGCZ00, RFSdQB, MGVT02]. However, all of 
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them were given the same footing and they could be rederived with aid of the SV criteria. There 
are several states whose NPT entanglement can not be detected with aid of these earlier criteria. 
Mostly these are states whose entanglement only becomes visible when considering moments 
of sufficiently high order. However, the criteria in [SimOO, DGCZOO, RFSdG03, MGVT02] only 
employ moments up to second order. The SV criteria, on the contrary, make use of moments of 
any order. 



2.4.5. Entanglement Witnessing in Two-Mode Cat States 

To illustrate the power of the SV criteria, entanglement detection in two-mode Schrodinger cats 
is demonstrated. So consider the two-mode cat state 

= -^(V,a)+e^ (2.118) 

with 

Af<p =2 + 2e- 4|a|2 cos tf>, (2.119) 

and < <p < 2tz. Ideally, the goal is to find a way to detect entanglement in |t/?<*) for any cp 
with aid of some witness or SV determinant. It is unknown whether such a witness or single 
determinant exists. However, a set of "witnesses" can be derived from which a feasible one can 
always be chosen such that entanglement witnessing for any [ips) can be performed. 



Concurrence 




ii 



Figure 2.3.: Concurrence of the two-mode cat state |^), (2.118). 



To show that the considered state actually is entangled for n ^ and arbitrary cp, the 
concurrence can be calculated. It is given by 



e -4kl 2 



C ^ = i 7 -4kR 1' (Z120) 

1 + e 4 l ft l cos(p 



and is shown in figure 2.3. It is noticeable that the set of states defined by \ipf) has two-sided 
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properties. The even two-mode cat \ipo) behaves other than the rest of the two-mode cats \ipq,) 
with <p 7^ 7i with regard to entanglement. \ipn) represents a maximally entangled two-qubit state 
^|(|01) + 1 10)) for any amplitude a. ^ [Wan02, vEHOl]. Hence it reveals highly discontinuous 
behavior at a = 0. On the contrary, all other states display continuous behavior. 

In [SV05a] Shchukin and Vogel demonstrated entanglement detection in \ipn) with aid of the 
principal minor 

l (S + ) («S + ) 
si := {To) {tfb) {afrb) . (2.121) 
(a f b) (a f b f b) (a f ab f b) 

Furthermore, it is known that the inseparability criterion by Duan et al. ([DGCZ00]) is able to 
detect entanglement in \ipo). However, the Duan criterion is just a weaker form of the criterion 
corresponding to the SV determinant 



S2 



(b) {ab) 



(a + S + ) 
{b + b) 



(2.122) 



as Shchukin and Vogel showed. So, what are the witnessing capabilities of these determinants 
with regard to the whole set of states {|i/ty) : < (p < 2n}? Evaluation of S\ and s 2 for general 
<p is straightforward: 



Si 



S2 



4|ft| 6 e 4 l a l 2 

( e 4|«| 2 +coS<^) 3 

4|ffi| 4 e 4 l g l 2 
(e 4 H 2 + cos<£) 3 



(1 



(1 



-e 4 !"! costf>), 



„4|a|* 



cos<p). 



(2.123) 
(2.124) 



Hence, determinant Sj is capable of witnessing entanglement in : j < <p < ^}, while s 2 

witnesses entanglement in {\ipip) '-0<<p<^;^j-<(p< 2tl\. This is graphically demonstrated 
in figure 2.4. So, define a "generalized determinant" 



s„ = nsi + (1 - n)s2 : n e {0,1}, 



(2.125) 



which is capable of entanglement detection in the whole set { \ip^,) : < (p < 2tz}. This is meant 
in such a way that for any there is n = V 1 such that s„ < 0, which is equivalent to 
successful entanglement detection. s"o applies for < (p < -j and ^ < <p <2n, while s~\ is the 
right choice for j < <p < For <p = j and <p = t?, both s"o and s"i are fine. However, the 
generalized determinant can be also directly formulated as a function of the phase <p in 
For any |t/ty), 

sj, = 0i [cos(<£ + 7t)\ si + 0i [cos(<p)] s 2 (2.126) 

performs entanglement witnessing for any <p. Here it has been made use of the Heaviside step 
function 0i , which is defined in appendix B.4. It can be concluded that the SV criteria are a very 
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Relative Phase <p Relative Phase <j) 



Figure 2.4.: The SV determinants Sj and S2 for |t/><*). While sj performs witnessing for j < (p < 
^r, S2 does its job for < <p < j and 3p < <p <2n. Hence, these two determinants 
are sufficient for witnessing in \xps) for any (p. Note that the blue and the red plane 
respectively denote the zero-plane. 

strong tool for entanglement detection, when there is an approximate idea of which principal 
minors to look at. The derived "generalized determinants" yields for the considered states \ip,p) 
of course no new results, since there is no witnessing necessary at all. The concurrence can 
be easily calculated for these states, which sets out all entanglement properties of the state. 
However, when the state is subject to a noisy channel, entanglement quantification is not so 
easy anymore. Then it is good to have a set of powerful witnesses or determinants in reserve to 
at least detect entanglement. 

Continuative questions to be tackled are whether the derived generalized determinants also 
work for more general two-mode CSS states, such as 



or even 




(2.127) 



(2.128) 



for general coherent states \ot), |/5), (7) and \S). 

Furthermore, as just mentioned, when states of this kind are subject to lossy or noisy channels, 
forcing them to become mixed, which are then the appropriate witnesses or determinants? 
Calculations suggest the conjecture that a witness (in this paragraph witness may refer not 
just to actual witness operators but also to determinants, etc.) which detects entanglement in 
a state \xp) is also the right one for entanglement detection in Y(|t^)), where Y denotes some 
channel (see calculations in subsections 4.3.1 and 4.4.1). This would not be so surprising, since 
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physical channels do not essentially reshape or transform entanglement. They can only destroy 
or preserve it. However, if not all entanglement has been detected by the witness in front of 
the channel, then it may of course not be able to detect anything afterwards, if the channel has 
destroyed exactly that part of the entanglement which has been initially used by the witness. 
But in this case there would exist better suited witnesses for the state at hand anyway. Namely 
these which would make use of all entanglement of the initial state. Admittedly, do such 
witnesses always exist? If yes, then the state after the channel should be entangled if and only if 
the witness, which is in this sense optimal, takes effect. The other way around, the question is 
the following: For any state p, does there exist a channel Y, such that a witness which detects 
entanglement in p in an optimal way, such that it makes use of all its entanglement, is not 
able to detect entanglement in Y(p), when there is actually still entanglement present? If not, 
the witness would become a necessary and sufficient tool for entanglement detection behind 
the channel. Or when the witness is in this latter sense not optimal, does there always exist a 
channel which preserves some entanglement, but prevents the witness from detecting it behind 
the channel? 

Concluding, there are a lot of interesting open questions regarding two-mode CSS states as 
well as general entanglement witnessing issues. However, these are not in the focus of this 
thesis. 



3. Optical Hybrid Approaches to Quantum 
Information 



The term hybrid is quite in vogue in quantum information. Several authors use it to describe 
various experiments and protocols. First, there are proposals considering hybrid quantum 
devices which combine elements from atomic and molecular physics as well as from quantum 
optics and also solid state physics [WHR+09]. Second, there is the notion of hybrid entanglement, 
when talking about entanglement between different degrees of freedom, for example, the 
entanglement between spatial and polarization modes [NLDS09, GAZ+10]. Finally, schemes 
which employ both CV and DV resources are also called hybrid [BPvLlO, LMSN08, WanOl, 
LloOO]. This is, how hybrid is defined in this thesis. 

Definition 3.1. A quantum information protocol is called hybrid if it utilizes both resources from 
discrete variable quantum information and continuous variable quantum information. These resources 
may include CV and DV states as well as CV and DV quantum operations and measurement techniques. 

For instance, a protocol is considered hybrid if it uses both qubits and qumodes or if it makes 
use of both DV photon detection and CV homodyne detection. In the preceding chapter it 
has been pointed out that the description of CV and DV states proceeds differently. Hence, 
combining CV and DV resources has its own characterizing and challenging subtleties. So, the 
question arises, why actually aim for hybrid approaches? Note that this chapter is mainly based 
on a recent review on optical hybrid approaches to quantum information by P. van Loock [vLll]. 

3.1. Why go Hybrid? 

In CV quantum computation, Gaussian states as well as Gaussian transformations, such as beam 
splitting and squeezing, are used. 1 Computational universality then is the ability to simulate 
any Hamiltonian which is expressed as an arbitrary polynomial in the mode operators with 
arbitrary precision. This kind of universality is called CV universality . However, to reach such 
universality in CV quantum computation just linear Gaussian elements are not sufficient. At 
least one non-Gaussian component is necessary. Actually, any quantum computer utilizing only 
linear elements could be efficiently simulated by a classical computer [BSBN02] . Unfortunately, 

^n this thesis no introduction into quantum computation is given. Hence for such introductions see for example the 
books by Nielsen and Chuang, as well as by Mermin [NCOO, Mer07]. A good overview over the concepts is also 
contained in [vLll]. 
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this single non-Gaussian element forms the problem with CV quantum computation. It is very 
difficult to efficiently realize non-Gaussian transformations on Gaussian states. 

In DV quantum computation the encoding of information takes place in a finite-dimensional 
subspace of the infinite-dimensional Fock space. Then the weaker DV universality refers to the 
ability to simulate any unitary acting on this working space with arbitrary precision. Just as 
in the CV case, for deterministic processing, a nonlinear interaction is required to realize DV 
universality. But when truncating the Fock space only single or few-photon states are left. The 
drawback on DV quantum computation is that nonlinear interactions on the few-photon level 
are hardly accomplishable. Note that there actually exists an efficient protocol for universal 
DV computation with only linear optics by Knill, Laflamme and Milburn (the KLM scheme 
[KLM01]), which, however, is probabilistic (or near-deterministic at the expense of complicated 
states entangled between many photons). 

A way out of the problems of CV and DV quantum computation may provide hybrid 
approaches. The GKP scheme by Gottesman, Kitaev and Preskill can be considered as one of the 
first hybrid protocols for quantum computation [GKP01]. It makes use of CV Gaussian states 
and transformations in combination with DV photon number measurements. So-called non- 
Gaussian phase states can be created from Gaussian two-mode squeezed states (TMSSs) with 
the aid of photon counting measurements. Additionally, the GKP scheme can be considered as 
hybrid, as it employs the concept of encoding logical DV qubits into CV qumodes. 

Gottesman, Kitaev and Preskill exploit so-called measurement-induced nonlinearities for the 
realization of the non-Gaussian element. When talking about hybrid schemes, it has to be 
distinguished between measurement- induced nonlinearities and weak nonlinearities. The latter is 
an actual nonlinearity, which may be enhanced or even mediated through suffenciently intense 
light fields. A quantum bus {qubus) can be used to mediate the nonlinear interaction between 
two possibly even distant qubits and act for example as an entangling gate [SNB+06, CZ95, 
vLMN + 08, vLLMN08]. A measurement-induced nonlinearity is no actual nonlinearity. The 
state to be nonlinearly transformed is entangled with an ancilla state, which is subsequently 
measured in such a way that the effect on the residual state corresponds to a nonlinear 
interaction. This concept is exploited in the seminal works by GKP and especially KLM. 

Concluding, the goal of hybrid approaches is to circumvent the limitations of the practical 
schemes for quantum computation. Nevertheless, their efficiency and feasibility shall be 
maintained as well as possible. For this purpose 

• hybrid states, operations and measurements, 

• qubus concepts, 

• weak nonlinearities, 

• as well as measurement-induced nonlinearities 
are used. 

Besides the quantum computational aspect, there is one more point to mention. CV and 
DV schemes both have their characterizing advantages and disadvantages. The heralding 



3.2. Applications involving Hybrid Approaches 



39 



mechanisms, necessary in DV schemes, make them highly probabilistic and hence rather 
inefficient. However, the fidelities in case of successful operations are very high, often near 
unity. In contrast, CV schemes are typically deterministic, as no heralding is required. Their 
drawback lies in the fact that CV encodings are very sensitive to losses and noise, which yields 
lower fidelities. Hence, there is typically some kind of trade-off between DV and CV schemes. 
Hybrid approaches may offer possibilities to benefit from the individual advantages while 
suppressing the disadvantages. 

The next section will present a few exemplary applications which involve hybrid protocols. 

3.2. Applications involving Hybrid Approaches 

A first example for a protocol making use of hybrid approaches has been already given in the 
previous section: The GKP scheme for quantum computation, which will be discussed in slightly 
more detail now and compared to the KLM scheme. It is a cluster-based quantum computation 
scheme and makes use of linear Gaussian CV resources and operations in combination with 
DV photon counting. 2 A TMSS is produced with a Cz entangling gate. Then a photon number 
measurement is performed which yields the so-called cubic phase state. With aid of this phase 
state, the cubic phase gate can be realized, the non-Gaussian element necessary for universal 
CV quantum computation. It is important to note that all states and gates can be prepared and 
performed "offline", prior to the actual computation, which then only consists of measurements. 
Actually, it is this offline cluster state preparation, which is the most problematic part in the 
protocol. Growing a sufficienty large cluster state with applicable properties is an extremely 
difficult task and could not be realized on large scales yet. 

As a comparison the KLM scheme is teleportation based. It is a fully DV model, exploiting 
linear optics and measurement induced nonlinearities by photon number measurements. The 
gates are teleported on the processed state, which, however, requires highly nonlinear entangled 
ancilla states. This is the drawback of the KLM scheme. While it is in principle efficient, it is 
highly impractical, as the demanded ancilla states are too complicated to be engineered with 
current technologies. 

Besides quantum computation there are several more elementary tasks which can be 
performed with hybrid approaches more efficiently than with only CV or DV schemes. 
One example is the realization of POVMs for optimal unambiguous state discrimination 
[Ban99, WTC+08]. 3 A 50/50-beam splitter, a coherent-state ancilla |«) and DV photon on/off 
detectors can be used for the optimal unambiguous discrimination between binary coherent 
states {|±a}}, see figure 3.1. 

Another example for a hybrid scheme is the generation of Schrodinger-cat states \ip±) = 

2 Cluster state quantum computation is a kind of "one-way measurement based quantum computation" which makes 
use of so-called cluster states. It has been invented by Raussendorf and Briegel, see [RB01]. For an introduction and 
review see [Nie06]. 

3 For an introduction to unambiguous state discrimination, where also positive operator-valued measures (POVMs) are 
explained, see [CheOO, Hel76]. 
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(on /off] 



(on/ off I 



A/"4(|a) =t |— a)) via DV photon subtraction of CV Gaussian squeezed vacuum. The subtraction 
procedure is performed by mixing the 
squeezed vacuum with one-photon or two- 
photon Fock states in a beam splitter fol- 
lowed by homodyne detection and postse- 
lection depending on the measurement result. 
One-photon subtraction yields even cat states, 
while two-photon subtraction results in odd 
cats [LJG+06, OJTBG07]. Another way of gen- 
erating Schrodinger-cats is based on hybrid 
entangled states [SBW90]. It is presented in 
chapter 5. CSS states find a very important ap- 
plication in fault-tolerant, universal quantum 
computation [LRH08]. 




Figure 3.1.: Optimal unambiguous state dis- 
crimination of equally likely bi- 
nary coherent states {|±a}} using 
a beam splitter, a coherent-state an- 
cilla | a) and DV photon on/off de- 
tectors. The inconclusive event here 
corresponds to the detection of the 
two-mode vacuum e — I* 1 00) with 
probability e -2 '"' 2 [vLll]. This fig- 
ure is a modified version of Fig. 
2.2.3 in [vLll]. The figure was used 
with permission by Peter van Loock. 



More examples for hybrid protocols can be 
found in the areas of quantum state charac- 
terization and quantum communication. In 
their theory paper [FC04] Fiurasek and Cerf 
show "how to measure squeezing and entan- 
glement of Gaussian states without homo- 
dyning". Instead of homodyne detection they 
apply DV photon counting to characterize 
the CV Gaussian states. Another theoretical 
proposal suggests, to measure the Bell nonlo- 
cality of CV Einstein-Podolsky-Rosen states 
(which actually are TMSSs) via measurement 
of the DV photon number parity operator 
[BW98]. However, even the simple CV homo- 
dyne tomography of DV one- ([LHA+01]) or two-photon ([OTBG06]) Fock states can be already 
considered as a hybrid experiment. In quantum communication there are entanglement con- 
centration and entanglement distillation subroutines which can be performed more efficiently 
with hybrid protocols. [OKW00] demonstrates, how to theoretically concentrate entanglement 
of pure CV TMSSs through DV photon subtraction. [ODTBG07, TNNT+09] present the ex- 
perimental execution. Entanglement distillation of noisy CV TMSSs using beam splitters and 
on /off DV photon detectors is theoretically set out in [EBSP04] and experimentally performed 
in [HSD+08, DLH+08]. 

A highly important tool being utilized in several models for quantum computation and 
quantum communication is the hybrid quantum bus concept [SNB+06, CZ95, vLMN + 08, 
vLLMN08]. The idea is, to use an intense qumode to mediate a nonlinear interaction (even 
when it is weak) between two qubits. It will be presented in chapter 5. 
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One relevant ingredient of several of the mentioned protocols is a new intriguing form of 
entanglement, the hybrid entanglement between CV and DV systems. What distinguishes this 
kind of entanglement from ordinary DV or CV entanglement? Can, just as in the CV case, 
where one distinguishes between Gaussian entangled states and non-Gaussian entangled states, 
different kinds of hybrid entanglement be identified? How is quantification and witnessing of 
hybrid entanglement achievable? Is there also multipartite hybrid entanglement? And how can 
hybrid entangled states actually be generated? The rest of this thesis addresses these questions 
and provides answers to most of them. Also some of the already mentioned applications of 
hybrid entanglement are presented. 



4. Hybrid Entanglement 



In this chapter hybrid entangled (HE) 1 states, i.e. entangled states whose constituents are both 
DV and CV systems are discussed and classified, while the emphasis lies on bipartite hybrid 
entanglement. The chapter starts with an introductory section 4.1 focussing on some general 
issues regarding HE. Afterwards, in section 4.2 pure bipartite HE states are investigated. Then 
we move on to mixed states. Section 4.3 discusses mixed states which are effectively fully DV 
while in section 4.4 so-called truly hybrid entangled states are introduced and examined. While 
the first 4 sections only focus on bipartite HE, finally, in section 4.5 also multipartite HE is 
explored. 



4.1. Introduction 

An example for a hybrid entangled state is a state like 

\TP) AB = ^(\0) A \rp ) B + \l) A \M B ), (41) 

where |0} A and \1} A are DV qubit states in V.^ and |i/'i) B as well as |i/>2) B are CV qumode states 
living in the infinite-dimensional Hilbert space H^. Already this state shows a combination 
of DV and CV properties. On the one hand, it can be seen that it has the same normalization 
constant 4| as the comparable, fully DV qubit state (| 00) + |11)). On the other hand, since 
| ipi ) and | fa) are general CV qumode states and hence not orthogonal, there does not directly 
exist a Schmidt decomposition of the state. This is also the case for pure fully CV states. 

However, one can think of many HE states like \ip) Ah ' , but with a possibly higher dimensional 
DV system or even mixed states. So, a rigorous definition of hybrid entanglement is necessary. 

Definition 4.1. Any entangled bipartite state of the form 

N N 
pAB = Y j p n \xp„) AB (Vnl , Pn > V n ; p n = 1, 

I'M = I] c nm \m) A \lp„ m ) B , C nm G C ; \ c nm\ 2 = l, 

m=0 m=0 

with generally non-orthogonal qumode states \ip n m) B which is living in a Hilbert space of the form 
V. AB = U A ® V-^c with finite d, is called hybrid entangled. 
'Depending on the context, HE may denote either "hybrid entangled" or "hybrid entanglement". 
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How can HE states actually be described in a useable way? Can their entanglement be 
quantified? Facing an overall infinite-dimensional Hilbert space Hf ® H.^, no density matrices 
can be employed. So, CV methods have to be attempted. But in CV entanglement theory the 
only conveniently representable states are the Gaussian ones. In the general non-Gaussian CV 
regime, for instance, exact entanglement quantification is not achievable. Therefore, do HE 
Gaussian states exist? 

Lemma 4.1. Any singlepartite d-dimensional quantum state with finite d and d > 1 is non-Gaussian. 

Proof. Consider such a general d-dimensional state in a pure state decomposition and make 
use of the Fock basis of a d-dimensional subspace of the Fock space, 

N N 

P = E Pn \ipn) (>pn\ , Pn > OV n; £ p n = 1, 

n=1 " (4.3) 

\tfn) ~ Cnm , Cnm £ C ; \Cnm\ ~ 1> 

m=0 m 

Non-Gaussianity corresponds to non-Gaussianity of the characteristic function. This is equiva- 
lent to non-Gaussianity of the Wigner function, as (non-)Gaussians stay (non-)Gaussians under 
Fourier transformations [BSMM05]. Hence, consider the Wigner function of p, 

CO 

W(x,p):=i / dye 1 ^ (x-y\p\x+y) 

— CO 

N °? d-1 

= E ~ / d v^ vy E c nl c* nk Ux-y)r k (x + y). 



(4.4) 



n=l 71 J m kl=0 



The position wavef unction of a Fock state |n) is given by 

tp n (x) = (x\n) = H »W e -4, (4.5) 

where H n (x) are the Hermite polynomials [Leo97, BSMM05]. For an overview over some relevant 
properties of the Hermite polynomials, see appendix B.5. Therefore, 

n=i 71 J u=o V2 k +'kll\n 

(4.6) 

CO ii , - v ' 

- 1 f dve 2i PV e-* 2 -y 2 f v Y c ,c* H k(^ + y)Hj(x - y) 

- - I dye*" e * l_ Pn ^ c ni c nk , ■ 

nj n= i a=0 V2 k +ikll\Tz 



Now i J dy e 2l Py is just a Fourier transform operator with respect to y, and e * V 2 is a Gaus- 

„* H k{x+y) H i( x - 

~nk ^2 k +ik\l'.TT 



sian. Hence, for W(x,p) being a Gaussian, the sum P := YLn-1 P"^kl=O c nl c n H ^ 1+y ^ H '^ y ^ 
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also has to be Gaussian in x as well as in y. However, P is just a polynomial. Strictly speaking, 
it is a polynomial of finite order. Hence, P can never be an exponential and even more no 
Gaussian. Admittedly, there is one exception. P may be a polynomial of order zero, i.e. a 
constant. 

So, assume that p is Gaussian and show that this is the case if and only if c n \ = V / > 1. 
Showing "<^=" is trivial, since for c n \ = V / > 1, p = |0) (0|, which is obviously Gaussian. For 
"=>" perform an induction with respect to the dimension d. 

• Inductive Basis: For d = 2 there is only one term x 2 ^ 1 ' = x 2 in the sum P, which comes 
from H d _i(x + y)H d ^ l (x - y) = H^x + y)H 1 (x - y). Its coefficient is c\Y^=\ Pn\c n ,i\, 
where c\ is a real constant due to the Hermite polynomials. If the state is Gaussian, 
Hn=i Vn\c n ,\\ = 0. Since p„ > V n and |c n< i| > V n, from J^ =1 p n \c n> i\ = follows 
c n,l =0V«. Hence, c n \ = V I > 1 (Z can be only or 1 for d = 2) is proved as a necessary 
condition for Gaussianity for d = 2. 

• Inductive Step: Assume validity of the statement for d = m and consider d = m + 
1. Then, there is only one term x 2 (( m+1 ) -1 ' = x 2m in the sum P, which comes from 
H (m+1) _ 1 (x + y)H( m+1 )_ x (x - y) = H m (x + y)H m (x - y). Again with c m being a real 
constant the coefficient is c^J^Lj Pn\c n ,m\- If the state is Gaussian, Pn\c n ,m\ = 0. 
Since p n > V n and \c ni , n \ > V n, from T,n=iPn\c n ,m\ = follows c„ im = V n. 
However, from the inductive hypothesis it is known that c n ;=0Vl</<d — 2 = m — 1. 
Hence, c n i=0\/l<l<d — 1 = m, which proves the validity of the statement for 
d = m + 1 and hence for all finite m. 

Concluding, a rf-dimensional quantum state with finite d and d > 2, written as (4.3), is Gaussian 
if and only if c„; = V / > 1, which corresponds to p = |0) (0|. Therefore, it can be said that p 
is Gaussian if and only if p = |0) (0|, which is only one-dimensional. Hence, any rf-dimensional 
quantum state with finite d and d > 2 is non-Gaussian. □ 

Theorem 4.1. Any bipartite hybrid entangled or classically correlated state is non-Gaussian. 

Proof. If a multipartite quantum state is Gaussian, all its subsystems have to be Gaussian. 
So, consider a state of the form (4.2) and trace out the CV subsystem. What is left over is a 
d-dimensional singlepartite system with finite d, which can be described in its Fock basis. It 
is denoted by p. Due to lemma 4.1, p is Gaussian if and only if p = |0) (0|, which is a pure 
state. However, for any entangled or even classically correlated state, the reduced state cannot 
be pure [Aud07]. Hence, for Gaussian p, the overall system cannot have been entangled or 
classically correlated. Therefore, every bipartite hybrid entangled or classically correlated state 
is non-Gaussian. □ 

The proof of theorem 4.1 basically relies on lemma 4.1, which states that any DV state with 
dimension > 2 is non-Gaussian. It is straightforwardly generalized to multipartite systems. 
Arguing that for Gaussianity all subsystems have to be Gaussian, it is sufficient for the non- 
Gaussianity of any multipartite system which also possesses DV constituent(s) that at least one 
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DV subsystem is of dimension > 2. However, if there shall be entanglement between the DV 
subsystem and the rest, the state is necessarily non-Gaussian, since for entanglement dimension 
> 2 is required. So, any multipartite quantum state which involves DV entangled subsystems 
is non-Gaussian. Even more generally, only 1- or infinite-dimensional systems can be Gaussian. 
This is of course not surprising, however, it is worth pointing out that theorem 4.1 is not trivial. 
Since the overall Hilbert space of HE systems is very well infinite-dimensional, it is not a priori 
clear that there do not exist Gaussian hybrid entangled states. 

As shown, HE states live in the non-Gaussian, infinite-dimensional Hilbert space regime, 
which is quite problematic as already known from conventional CV entanglement theory. The 
states can be neither described with proper density matrices, nor with covariance matrices. 
Phase space representations are also not very suitable, since one of the subsystems is DV. 
They are not qualified for investigations regarding entanglement properties anyway. The only 
known quasi-probability distribution which may in some cases make statements about the 
entanglement of the state is the Glauber-Sudarshan P-representation. However, it can be easily 
shown that this function does not even exist for these highly nonclassical states. Of course, 
it may still be possible to perform entanglement witnessing, but what about quantification? 
Is there a way out of the dilemma? Yes, there is! For some HE states the unique Hilbert 
space structure can be exploited in such a way that the states can nevertheless be described by 
density matrices. These states are therefore called DV-like hybrid entangled. This gives rise to a 
classification scheme of HE states. 

4.1.1. Classification 

Consider again a general HE state of the form (4.2). Depending on the number of mix terms N 
and the dimension of the DV subsystem d, there can be maximally N x d linearly independent 
CV qumode states \ip n m) B hi p AB . d is always finite due to the definition of HE. However, N 
may be either finite or infinite and, hence, the number of linearly independent qumode states is 
either finite or infinite. Furthermore, if N — 1, the state is pure. 

If the number N x d of linearly independent CV qumode states is finite, they only span a 
N x (/-dimensional subspace ?{ Nxi j of the initially infinite-dimensional Hilbert space "Hoo- Then, 
the Gram-Schmidt process can be employed to express the qumode states in an orthonormal basis 
of this finite-dimensional subspace. In this case, the state is effectively DV and all the methods 
from DV entanglement theory can be applied. 

However, this shall be discussed in more detail. What actually is the Gram-Schmidt process? 
It is a method for orthonormalizing a finite, linearly independent set of vectors in an inner 
product space [KS98, TI97]. 2 For a linearly independent set of vectors {\ipi) : i = I, . . .,n} (since 
the process is to be exploited in the framework of Hilbert spaces, the inner product space is a 
priorily assumed to be a Hilbert space) a set of pairwise orthonormal vectors {|e,) : i = !,...,«} 



2 Note that the Gram-Schmidt process and the Schmidt decomposition are both named after the German mathematician 
Erhard Schmidt. He developed the Gram-Schmidt process together with Jorgen Pedersen Gram. 
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spanning the same subspace as {|t/>,) : i — 1, . . ., n} is given by 



ki> = Ifc) , 

14) = |^2>-<ci|^2>kl> , 



ki> 
k 2 > 



14) 



(4.7) 



n-i 



14) = I'M - E teliM k«'> ' kn) = 



Making use of this, any finite set of linearly independent qumode states can be expressed in an 
orthonormal basis { |e,) : i — 1, . . . , n}. But strictly speaking, equation (4.7) only tells, how to 
express the new orthonormal basis in terms of the old non-orthonormal one. What is actually 
required, is the other way around. To express the qumode states in terms of {|e,-) : i — 1, . . . , n} 
try the following approach, which can be considered as some kind of inverse Gram-Schmidt 
process. 

Theorem 4.2. n normalized, in general non-orthogonal, linearly independent states {\ipi) : i = 
l,...,n; < | (rpi\ipj) | < lVz,/} can always be expressed as \xpj) — YL l j=i a ij\ e j)> where : 
i — 1, . . . , n} form an orthonormal basis of the space spanned by { | tpj) } and a^ e C. 

Proof. Writing out \tpi) = Ylj=\ a ij \ e j), theorem 4.2 states that the {\~ipi}} can be always written 
as 

|tfl) = «ll kl> , 

\lp 2 ) = fl 2 l kl> +«22 k2) , 

IV3) = fl 31 kl> +«32 |e 2 ) +«33 k3> , , s 

(4.8) 



I V") = a nl l e l) + •■■ + %« l e n) = J2 a "' ki) 



i=l 



To prove it, it has to be shown that 1.) \ipi) = Ylj=\ a ij \ e j) corresponds to a valid basis 
transformation and 2.) it actually performs the right mapping. 
1.) Write the transformation as 

|^>=E^/k;)' (4-9) 



with the transformation matrix 



U = 



( 011 

«21 a 22 

\a„i a n2 



\ 



(4.10) 
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Due to normalization of the initial and the new vectors Ylj—i \ a ij\ 2 = 1 is known and due to the 
linear independence of the {|i/>;)} also flj, ^= Vz is necessary. 

=> det[U]=n? = i«ii^O. 

Lf is invertible. 

=>■ U is a valid basis transformation. 

2.) To show that the lower triangular structure of the basis transformation U in combination 
with the orthonormal basis {|e,)} is sufficient to actually express the accurately in terms 

of {|e,}}, it is demonstrated that the N (^ +1 ) parameters fly can be chosen such that all overlaps 
(ipi\ipj) are preserved when the transformation is applied. 

On the one hand, there are n 2 such overlaps in total and ones with potentially differing 

absolute values. On the other hand, there are complex parameters fly. From the structure 

of the basis transformation and the fact that {fly} are complex, it is clear that if fly can be 
chosen such that the -element set of : i < /'} can be preserved, also the rest of the 

overlaps is preserved, since they are only complex conjugates of the former. Hence, it already 
becomes reasonable that the fly can be chosen appropriately. 

However, a proper proof is performed by induction in n: 

• Inductive Basis: 

n = 1: There is only one overlap to be preserved. 

(VilVi) = 1 = (ei|fln«n|ei) = |«n| 2 . (4.11) 
Hence, choose a n = 1, which preserves the overlap (ipi\tpi). 

• Inductive Step: 

Assume fly have been calculated for i < n — 1 such that all overlaps { {ipi\ipj) ■ i,j < n — 1} 
are preserved. We show that then also all a n j can be chosen such that the overlaps { (ipj\ip n ) : 
i = 1, . . . , n} are preserved. The complex conjugated overlaps follow automatically as 
argued before. 

Applying the basis transformation, the overlaps { (xpj\ip n ) : i — 1, . . . , n — 1} are 

{ipt\ip„) = a* n a nl , 
(ip 2 \tpn) = a* lx a n \ +al 2 a n 2, 

(4.12) 

{lpn-\\1pn) = 0* n -l,l a nl + ■ ■ ■ + a *n-l,n- 
As fly have been calculated for i < n — 1 due to the inductive hypothesis, this is just a 
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system of linear equations, which can be written as an augmented matrix: 



*21 



"22 



o \ 



a *n-l,n-\/ 



(4.13) 



Since a ti ^ V i, this system of equations is exactly solvable. Hence, {a n j : j — 1, . . . , n — 1} 
can be chosen such that {{ipi\ipn) '■ i = 1, ■ ■ ■ ,n — 1} are preserved. Therefore, there are 
only one free parameter a nn and one overlap (rp n \ipn) to be preserved left: 



N 



N-1 



(lpn\lpn) = 1 = E l fl N/| 2 = \ a Nj\ 2 + I«Nn| 2 - 
j=l j=l 



(4.14) 



From Y!j=\ \ a nj\ 2 = 1 an d a nn ^ 0, which is already known, Ejl/ |«N;'| 2 < 1 follows, and 
hence a nn can be chosen as 



inn 



1 



N-1 

i- E M 2 . 

7=1 



(4.15) 



In the end, also {\p n \ip n ) can be preserved. Therefore, the theorem is valid for n under the 
assumption of validity for n — 1. 



Concluding, with the inductive basis it is valid for all n. 



□ 



The proof has been performed in such a great detail, because it also sets out, how to actually 
compute {fly} for a given set of qumode states. Recalling equations (4.8) and (4.13), the {fly} 
can be calculated successively one after another by considering successive overlaps. A parameter 
fl/l is directly obtained from the overlap (ip\\ipj). Then, a\i follows from (1^2 l^i) together with 
the known a^. Likewise, a& is calculated from (ifclipi), an and F° r me other parameters 
just go on like this. Hence, the inverse Gram-Schmidt process can be efficiently implemented 
and computed. 

As an example consider the normalized qumode states {\ipi) : i — 1, . . .,3} with overlaps 



(</>ll</>3) = c 2 , 
(^2 1^3) = c 3 . 



(4.16) 
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They can be expressed in the orthonormal basis { |e,) : i — 1, . . . ,3} as 

1^2) = ci kx) + \Jl- \c\\ 2 \e 2 ), 

Ifo) =c 2 | ei )+ y^'p 1*2) + ^1 

The Gram-Schmidt process and the inverse Gram-Schmidt process respectively are possibly the 
most important tools with respect to HE quantum states. 

So, once again, consider a state of the form (4.2). For N — 1 the state is pure. Then it 
contains only d linearly independent qumode states, which span a d-dimensional subspace H^. 
Therefore, with aid of the inverse Gram-Schmidt process these qumode states can be expressed 
in an orthonormal basis. Then density matrices can be employed for the description of the 
overall state, and also a Schmidt decomposition can be performed or pure state measures such 
as the entropy of entanglement can be calculated. The state is effectively DV. For 1 < N < 00, 
the state is mixed. Nevertheless, it possesses only N x d qumode states, which can be again 
casted in an orthonormal basis. Therefore, also states of this kind are effectively DV and the 
density matrix formalism can be exploited. However, they are not pure anymore and neither 
pure state measures nor a Schmidt decomposition can be applied. Finally, there is N — 00. Here, 
by N = 00 states those states are meant which can be expressed with infinite N only and not 
also with finite N. Remember that pure state decompositions are not unique. Anyway, these 
states hold an infinite number of qumode states, which has the effect that the Gram-Schmidt 
process cannot be applied anymore. Hence, they are not describable with density matrices and 
therefore not effectively DV. Nevertheless, one subsystem is DV. These are the states which are 
called truly hybrid entangled. 

Summing up, for bipartite HE states in pure state decomposition with N denoting the 
number of mix terms in the convex combination of pure state projectors, there is the following 
classification scheme: 

• N — 1 : Pure hybrid entangled states 

- Supported by finite-dimensional subspace. 

- => DV-like entanglement. 

- Schmidt decomposition compatible. 

- DV Pure state measures (Entropy of Entanglement) applicable. 

• 1 < N < 00 : Mixed hybrid entangled states 

- Supported by finite-dimensional subspace. 

- => DV-like entanglement. 

- DV mixed state measures applicable. 



(4.17) 



F2 



|c 3 - qc 2 | 

1-kil 2 



k 3 > 
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• N = oo : Mixed hybrid entangled states 

- No support by finite-dimensional subspace. 

- => True hybrid entanglement. 

- No measures applicable. 

- Witnesses adaptable from CV entanglement theory. 

At first, it has to be pointed out that the possibility of applying DV methods on such a 
wide class of hybrid entangled states is quite remarkable. An initially non-Gaussian infinite- 
dimensional quantum state, which seems rather awkward at first sight, can finally be described 
with neat density matrices and in the pure state case even a Schmidt decomposition can be 
performed. Nevertheless, there is also the class of states which stay truly hybrid entangled and 
cannot be transformed with the Gram-Schmidt process. As mentioned earlier, in this infinite- 
dimensional, non-Gaussian regime exact entanglement quantification is up to now impossible. 
Therefore, witnesses which detect hybrid entanglement are here especially demanded. In the 
next subsection an adapted version of the SV criteria is presented, which can be used for this 
purpose. The other sections of the chapter will investigate relevant states from each class of HE 
states. 



4.1.2. Adapting Shchukin-Vogel Inseparability Criteria 

In subsection 2.4.4 the SV inseparability criteria have been introduced. They are well defined 
for bipartite CV states. However, HE states are not fully CV, but contain a DV subsystem. 
Nevertheless, the SV criteria can be applied when appropriately adapted to the new situation. 
Actually, there are two ways for the adaption. 

The first way is to simply interpret the DV subsystem as living in a subspace of an infinite- 
dimensional Hilbert space. The DV qudit is interpreted as a CV system supported by "Hoo and 
being encoded in a Fock basis. However, it only makes use of a finite number of the basis 
vectors. Then the SV criteria can be applied just as usual. 

The second way is not to adapt the state, but the criteria. The mode operators for CV states, 
which SV make use of, are defined by equations (2.23) and (2.24). Now, however, assume the 
system which the operators a and a + belong to to be d-dimensional. Then the orthonormal 
Hilbert space basis vectors \m) can be written as column vectors with a "1" in row m with 
m — 0, . . .,d — 1: 

\m) d = (0 ,...,0 m _ 1 ,l m ,0 m+1 ,...,0 d _ 1 ) T . (4.18) 
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With this notation the new audit mode operators a d can be defined as proper d x d matrices: 



fo VT ••• 

^/2 




v° 

They have the properties 



Vd^l 




/o 

VT o 

: y/2 







Vd^l OJ 



a d \n) = \fn \n — 1) , 
flj |n) = (1 - S dm )\fn + l\n + 1) , 
(«d) rf = 0, 

(4) d = o. 



Furthermore, the commutator changes: 



A 



V 



dxd 



-(d-l)J 



dxd 



(4.19) 



(4.20) 
(4.21) 
(4.22) 
(4.23) 



(4.24) 



Since any d-dimensional qudit \ip) d can still be written as \\p) d = |0) rf and also the 

projection operator on |0) d (0| nevertheless is : e~ a ^ ad :, the SV criteria can be also derived for 
hybrid systems with aid of these new operators. Structurally the same criteria are obtained, 
just that for certain hAiih'h me moments Mjj(p) = (fl + ' 2 a 1 a +;i di 2 b Vi V 3 i f ' 3 V A ) p are zero. Due 
to the properties (4.22) and (4.23) any combination of i\,i2>h>h sucn t na t (^d)* or i^dY with 
k > d occurs in the moments, nullifies them. In the end, in the matrix of moments the rows and 
columns corresponding to these combinations of i\, i2,j\>h are simply missing. 



When actually applying the SV criteria it rarely makes a difference whether the first or the 
second approach to the adaption of the criteria is chosen. Only for moments involving terms 
like (a d ) k (a^y with k,l £ N the approaches may yield diverging results. However, for all 
determinants and moments used throughout this thesis, they lead to the same results. 

As an example for the application of the SV criteria on HE, entanglement detection in the HE 
state 

\ip) = j=(\0)\«) + \l)\-*)) (4.25) 
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is shortly demonstrated. In subsection 2.4.5 the determinants 



si 



1 

(b) 
(a f b) 



(b f b) 
(a f h f h) 



(aP) 
(aPb) 
(a f aPh) 



and S2 = 



(a) (P) 



(d f a) 



t£t\ 



(a*b 



(b) (ah) (Ph) 



(4.26) 



were employed to detect entanglement in two-mode CSSs, which are structurally comparable. 
Hence, it is reasonable to try the same criteria on the HE state (4.25). The mode operators a and 
« + shall correspond to the CV subsystem, while b and P shall act on the qudit. This results in 



si 



kl 2 e- 4 H 2 



<0 Va^O, 



s 2 = ^L(l- e - 4 H 2 ) >0 V«. 



(4.27) 
(4.28) 



As S2 does not perform any entanglement detection, it can be concluded that second moments 
are probably not sufficient for the detection of entanglement in (4.25). However, s\, which goes 
up to fourth order moments, does the job very well, as it detects entanglement for all n^O. 
For a = 0, the state becomes a product state and hence it is not entangled. The determinant S\ 
is actually the main tool for entanglement detection in hybrid entangled states in this thesis. In 
the pure state case it is of course not so relevant, but especially for the detection of true HE s\ 
proves to be very valuable, as seen in section 4.4. 



4.2. Pure States 



In this section the first, pure state class of HE states is investigated. Some representative 
bipartite pure HE states are discussed and the mentioned tools are exploited. A general pure 
bipartite HE state in Uj <g> looks like 

m AB = E* Vf . a e c , E N 2 = i- ( 4 - 2 9) 

i=0 i=0 

Due to the finite number of qumode states, a Gram-Schmidt process can be always applied. 
Afterwards, also a Schmidt decomposition can be carried out. Finally, the state's entanglement 
properties can be also analyzed with aid of a pure state measure such as the entropy of 
entanglement. 

4.2.1. Pure Qubit - Qumode Entanglement 

Consider the case where the DV system is a qubit: 

|^^ = c |0) A |^o> B +ci|l) A |Vi> B . (4.30) 



54 



4. Hybrid Entanglement 



Since only the relative phase in the superposition is important and due to normalization, this 
can be written as 

\^ B E = V~c\0) A \iPo) B +e'1'VT^~c\l) A \ip 1 ) B , (4.31) 
with real constant c 6 [0,1]. Performing an inverse Gram-Schmidt process on system B yields 

l^o> B = |0> B , (4.32) 



= N\0y + ^\-\N\ 2 \\) a , (4.33) 

where the overlap (tpo|l/>i) is denoted by N. Making use of equations (4.32) and (4.33) \ipi) B 
can be expressed in the new orthogonal basis and DV methods can be applied to calculate the 
entropy of entanglement of | V)hE' w hich is shown in figure 4.1. 3 



Es(\m E ) = -A+log 2 A+ - A_log 2 A_, A± = - (l ± y/l - 4c(l - c)(l - |N| 2 )J. (4.34) 




Figure 4.1.: Entropy of entanglement of the state \tf) HE , (4.31). N denotes the overlap (ipo\xpi). 

Maximal entanglement £5 = 1 is achieved for orthogonal qumode states (N = 0) 
and c — h • Note that the entanglement only depends on the absolute value of the 
overlap N, not on its phase. 

Based on this result, compare the HE state \ip) AB to the structurally comparable, fully DV 
and CV states 

\ip) AB v = sfc |0) A |0) B + e**vT^ |1) A |1) B , (4.35) 
\q) A l = -^yi^ |d> ) B + e'VT^|^i> A l*i> B ), (4.36) 

3 Of course the Gram-Schmidt procedure is not really necessary here, since the CV subsystem can be simply traced 
out, and the entropy of the residual DV state can be ordinarily computed. For the calculation of the entropy of 
entanglement it is sufficient that one of the subsystems is DV. 
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where N$ is a normalization constant. It is known that \ip) cv 's entanglement depends on the 

AB 

phase <p, as also seen in subsection 2.4.5. On the contrary \tp) DV * s equally entangled for any 
phase (p. This the hybrid state has in common with the DV state. Its entanglement also shows 
no dependence on <p. A common property of all the states is that the entanglement is highest 
for c = 1. In this case, the DV state becomes a maximally entangled Bell state. Furthermore, 

AB AB AB A 

note that actually \tp) HE and \tf) DV are onr y special cases of \ip) cv '- For orthogonal \\pi) 's the 
HE state is obtained, while both orthogonal |^,} A 's and |<E>j) B 's result in \ip) D y. 

4.2.2. Squeezed Hybrid Entangled States 

In subsection 4.1.2 the state 

\ip) AB = ^(\0) A \«) B + \l) A \-«) B ) (4.37) 

has been investigated and its entanglement has been detected with aid of the SV determinant 
Sj. This actually is an example for a typical HE state which also occurs in QKD experiments 
(see chapter 5). Analyzing this state is pretty trivial, as it is just a special case of the HE state 
considered in the previous subsection. Regarding the SV criteria, so far only their advantages 
have been pointed out. Now it is demonstrated how easy the criteria can be forced to fail in the 
detection of entanglement. To do this, consider again \ip) AB , but this time squeezed in mode B: 

S B (Z) \ip) AB = -L (|0) A S B (Z) \*) B + \1) A S B (Z) \-a) B ), (4.38) 

with the squeezing parameter £ — re w . This state is still a special case of the general qubit- 
qumode states |j/>) H f considered in the last subsection. Its entanglement quantification is 
straightforward. 

The state is obtained by just applying a unitary transformation, i.e. the squeezing operation, 
on \ty) AB ■ As unitary transformations are a special kind of LOCC and also invertible, entan- 
glement has to be preserved under them. Therefore, \ty) AB necessarily has the same amount 
of entanglement as S B (£) \ty) AB ■ Instead of the latter, squeezed state, regarding entanglement, 
it is sufficient to examine the former, unsqueezed state. As the SV determinant s\ detects 
entanglement in \ty) AB for all « ^ 0, the same would be expected, when evaluating s\ with 
respect to S B (£) \ip) AB , for any squeezing £. 

However, 

Si(S B (£) \ip) AB ) = -sinh 2 r- ^ ^— |a| 2 cosh 2 r- ^^sinh 2 r, (4.39) 
4 2 8 

which is obviously not below zero for all r and « ^ 0. Also note the independence of the 
squeezing phase 9. See figure 4.2 for a visualization of s-[(S B (^) \ip) AB ). What is even more 
surprising is that, on the one side, the higher a is the more difficult it is to witness entanglement. 
But on the other side, the entanglement does actually increase as a gets higher, since S(£) B \a) B 
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Squeezing r 



Inconclusive Parameter Region 



Entanglement Verified 



Amplitude a 



Figure 4.2.: The left-hand graphic shows the SV determinant sj of S B (£) \ip} AE . The yellow 
plane denotes the zero-plane. For values (r,a) such that S\(r, oc) < entanglement 
is verified. The two regimes are plotted in the right-hand diagram. 



and S(£) |— a) approach orthogonality in this case. This is of course no contradiction to 
earlier results. Failing to detect entanglement does not mean there is not any. Nevertheless, it is 
remarkable that a simple unitary transformation, which does not change the entanglement at all, 
is able to cause the SV criteria to fail. While for r — unambiguous witnessing for all a. 7^ is 
obtained, for r ^ there are regions, where the SV determinant breaks down. It could be argued 
that the unitary transformation does on the one hand not destroy entanglement, but on the 
other hand possibly rearranges the existing entanglement in such a way that initially successful 
witnesses fail to detect the rearranged entanglement. However, this is just speculation. 

In S B (£) \ip) AB , the terms S B (£) \±a) B can be expressed as S B (£)D B (±a) |0) B . But what 
about states which involve squeezed qumode states in the form D B (±a)S B (£) |0) B instead 
of S B (£)D B (±a) |0) B ? Since S B (£) f D B (a)S B (£) = D B (£) with £ = acoshr + ct*e w smhr, and 
therefore D B (a)S B (£) = S B (£)D B (/3), these states can then be rewritten in the former way 
[GA90] . Therefore, they do not provide any new insights, as both scenarios are equivalent. 



4.2.3. Pure Qutrit - Qumode Entanglement 



As a final example of pure bipartite HE, an entangled state of a qutrit and a qumode subsystem 
is considered. Using this example all the mentioned tools for pure bipartite HE states are 
successively demonstrated. Consider 

WhI = "4 (|0> A \vac) B + \1) A \.) B + \2) A (4.40) 
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• An Inverse Gram-Schmidt process with respect to subsystem B yields 



vac 



- |o) D , 

a) B =x|0) B + Vl-x 2 |l) B , 

oc) B = x|0) B -xVl-x 2 |l) B + y/l -x 2 -x 4 + x 6 |2) B , 



(4.41) 
(4.42) 
(4.43) 



_ 1 I |2 

with x = e 2\ a \ . Hence 



AB 



- -J=(|0) A \0) B +x\l) A |0) B + y/l -x 2 \1) A |1) B +x\2) A |0) 2 
v3 ^ 

x 2 Vl-x 2 \2) A |1) B + V / l-x 2 -x 4 +x 6 |2) A |2) B ), 



(4.44) 



which is the effective DV form of the state. 

The pure-state entropy of entanglement of the state can be calculated. Its analytic form is 
rather lengthy. Hence, only the graphical version is shown in figure 4.3, from which it can 



Entropy of Entanglement 




Towards Maximal Entanglement 



Amplitude a 



Figure 4.3.: Entropy of entanglement of the state \ip) AB , (4.40). For a — » oo the state becomes 
maximally entangled. 



be seen that for a — » oo the state becomes maximally entangled. Note that, while the state 
deals with qutrit entanglement, the entropy of entanglement has been calculated in qubit 
entanglement units, so-called e-bits. 4 If the logarithm in the von-Neumarrn entropy had 
been calculated with respect to base 3, the maximal qutrit entanglement would actually be 
1 in these qutrit entanglement units and in the plot, the graph would converge against 1 . 



Finally, set x — |, which corresponds to a = V21n2 ~ 1.18, and calculate the Schmidt 
decomposition: 



\<p) AB = 0.76 \0} A |0) B + 0.56 |1) A |1) B + 0.33 \2} A |2) 



(4.45) 



4 A maximally entangled two-qubit Bell state \&)bc11 := ^j(I^O) + |H)) nas exactly 1 e-bit of entanglement. E-bits 
are therefore a unit for entanglement quantification, especially for two-qubit entanglement. 
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What can be concluded from not only considering qubits, but also qutrits as the DV subsystem 
of HE states, is that focussing on higher dimensional DV subsystems will probably not bring 
new insights. The methods applied are always the same whether the DV system is 2-, 3- or 
higher-dimensional. Hence, in the following mostly qubit subsystems are employed. 



4.3. Mixed States supported by a Finite-Dimensional Subspace 

In this section mixed HE states which are effectively DV are investigated. That is to say those 
states with 1 < N < oo according to equation (4.2). One the one hand, states are discussed 
which feature the characteristic properties of this kind of HE entanglement. Some typical 
DV tools are then applied on them. For example the very important entanglement measure 
logarithmic negativity is calculated. On the other hand, adequate mixed HE states which are 
of physical relevance are discussed. This leads to the question, how states actually become 
mixed. The answer are quantum channels which cause them to decohere. Therefore, along the 
way, also some interesting properties of quantum channels are presented, especially regarding 
entanglement evolution. 

What can be generally said about bipartite, effectively DV, mixed HE states is that there always 
an inverse Gram-Schmidt process can be performed, which yields a density matrix description. 
This enables the application of DV tools, such as typical DV entanglement measures. 



4.3.1. Amplitude Damped Qubit - Qumode Entanglement 

The state 



is not only a very simple and illustrating example of HE, but also physically relevant, as already 
mentioned. Therefore, simply send this state through some decoherence-causing channel and 
obtain a HE mixed state. Then it depends on the form of the channel whether the output is either 
effectively DV HE or true HE. As bipartite systems are employed, the question arises whether 
both subsystems or only one of them shall be subject to a channel. As the focus here lies not on 
complicated states and deep mathematics, but on the identification of characteristic effects and 
behavior, only one-sided channels p' AB = (1 <g> Y)p AB are considered. Two-sided channels will 
just make the outcome more intricate and the analysis more complicated. Furthermore, in a 
typical physical application, such as particular cases of QKD, it is only the CV subsystem which 
undergoes the action of the channel. Hence, consider channels of the form p' AB = (1 A (g) Y B )p AB , 




(4.46) 
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where subsystem A is DV and subsystem B is CV. 



p' AB = (l A ®Y B )l(\0) A {0\®\ip ) B {rp \ 



+ |1) A <1|® | 

+ \0) A (1\® \<po) B (ipi\ 
+ \l} A (0\®\xp 1 ) B (xp \ 



(4.47) 



Now write the state as a proper density matrix with respect to subsystem A, which is not 

A \p ,AB \r 

Y B {\^) B {^\) Y B (\^) B {^\i 



subject to the channel, i.e. p' AB = (i A \p' AB \j A ). 



p' AB = 



K Y B (\fi) B (rpo\) Y B (|^> B (^|) y 



(4.48) 



What can be seen is that for these one-sided channels on qubit-qumode HE states, a nice block 
structure is obtained when applying the density matrix formalism on subsystem A. This can be 
generalized to universal pure, qudit-qumode entangled states. For a general state of the form 



\<P) AB = X> V) A m B , qec, £ | Cl | 2 = l, 

!=0 i=0 



(4.49) 



the transformed state p ,AB — (t A <g> Y B ) \ip) AB (tp\ AB , in the density matrix formalism for 
subsystem A, looks like 



p' AB = 



I |c | 2 Y B (|^o> B (iPol) 
Vc (i _ 1 c*Y B (|^_ 1 > B (^o| 



c c*_ 1 Y B (|^ > B (^-il) \ 

|c d -il 2 Y B (|^_ 1 > B (^_ 1 |)/ 



or in index notation 



(4.50) 



(4.51) 



One step further, for completely general bipartite mixed HE states of the form (4.2), the output 
of one-sided channels (t A <g> Y B )p AB is 



N 



p'f=j:PnCn<jY B (\lp m ) B (%,W 

n=l 



(4.52) 



Consider a specific, physically relevant example. Choose as the state subject to the channel the 
already known 

i AB 1 



l^ B = ^(|o> A k> B + |i) A |- a > B ), 



(4.53) 
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and as the one-sided channel itself the amplitude damping channel (also called photon loss 
channel) acting on subsystem B. Exploiting Stinespring's dilation theorem 2.1, write 

p AB = tr E [U BE (\<p) AB (4>\ ® |0} E <0|) U BEf }. (4.54) 

The ancilla Hilbert space is the environment E, which is in a vacuum state, and into which 
photons may leak out of the qumode system. The coupling to the environment is decribed 
by a beam splitter U BE = e e ( fl E flB ~ fl B flE ), where the beam splitter transmissivity r] is given as 
cos 2 9 — tj. For a graphical illustration of the channel modeling, see figure 4.4. For the actual 




Figure 4.4.: Illustration of the amplitude damping channel. The output p ol ,t is obtained by using 
a beam splitter to mix the input state \ip)i n with vacuum environment. Afterwards 
the environment is traced out. 



calculation make use of U B ^U BE = 1, U BE D B (a)U BEi = D B {^/tja)D E (^/\ - no.) as well as 
U BE |0) B |0) £ = |0) B |0) E and obtain 



p AB = \{\Q) A {0\®\^) B {^a\ 



1> A (1|®|-^«> B (-^«I 

(4.55) 



+ e -2(l-,)|«| | ) A (1|® \^/t]a) B (-y^a\ 
+ e -2(i-7)|«|* | X )^ ( | ® \-^a) B (y/rfctlj 



This clearly is a mixed state. The channel causes the state to decohere, as can be inferred from 
the damping of its off-diagonal coherence terms. It can be also expressed in its pure state 
decomposition. Defining the damping t := e -2 ' 1- '?'!"! 2 and the overlap A := ( — ^/ffa\ y/rja) — 

~AB 1 + T | \AB I . 1 — T . \ AB / . i 
P = \<P+) (4>+\ + W-) (4>-\> 

, N (4.56) 
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This state contains obviously only two qumode states \±y/tja.) , which can be put into an 
orthogonal basis. 



\^ij*) B = \0) B , 

\-^ja) B =A|0) B + Vl - A 2 |1> B . 
Therefore, p AB is an effective two-qubit state with density matrix 

/ 1 





At 



At tVI - A 2 \ 

A 2 Ky/l - A 2 



\rVl - A 2 AVl - A 2 1 — A 2 J 



(4.57) 
(4.58) 



(4.59) 



When dealing with mixed two-qubit states, the obvious choice for entanglement quantification 
is the concurrence C. Its calculation is straightforward: 



C(p 



AB\ 



-\JlZ e - 4 '7M 2 (\/l +3 e - 4 ( 1 -'7)M 2 - \Jl - e - 4 ( 1 -'7)l a l ; 



(4.60) 



The plot is shown in figure 4.5. It shows an interesting trade-off behavior. On the one side, 




Concurrence 



Figure 4.5.: Concurrence of the state p , (4.56), as a function of transmissivity rj and amplitude 
a. On the one hand, the higher a the greater the initial entanglement, but also the 
more sensitive the state is against photon loss. On the other hand, for low a there 
is only little initial entanglement. However, the state is more robust against losses. 

Hence, there is some trade-off behavior and depending on rj an optimal ol^ exists 
for which the output entanglement is largest. 



a higher a results in a greater initial entanglement. But it also causes the state to be more 
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sensitive against photon loss. On the other side, a low a yields only little initial entanglement, 
but the state is more robust against loss. This leads to an optimal //-dependent ol^ for which 
the largest output entanglement is obtained. It is worth mentioning that a similar trade-off is 
observed in two-mode cat states J^- (\a) \a when one of the modes suffers 

photon loss. 

Finally, note that just as expected in the limit n = 1 the concurrence becomes C(p AB ) = 
y/l — e~ 4 l"l 2 , which is the concurrence of the initial undisturbed state 

\xp) AB = i=(|0> A |a> B + |l) A |-a) B Y (4.61) 



V2 

4.3.2. The Choi-Jamiolkowski Isomorphism and Entanglement Evolution 

This subsection can be regarded as some kind of bonus. As quantum channels play an important 
role in this thesis and have just been used for the first time, an elegant way to work with them 
and study them is presented. Consider a linear map between two finite dimensional Hilbert 
spaces: 

Y : S(H dl ) ^ S(H d2 ). (4.62) 
Then there can be defined an isomorphism between this map and an operator living in 

n Al <g> n d2 by 

p Y = (l^Y B )|cI>j/ B <cl>+|, (4.63) 
where l^^)" 48 denotes the maximally entangled state in ® Ufa, 

l* d + / B = ^T EV|i> B . (4-64) 
V«l i=o 

This Choi-Jamiolkowski isomorphism holds for all linear maps and operators. If, additionally, the 
linear map is completely positive, it is a quantum channel. The isomorphism was first suggested 
by Jamiolkowski in 1972 [Jam72]. 3 years later Choi employed it to derive the operator-sum 
decomposition of CPTP maps (which had been found independently by Kraus already in 1971) 
[Cho75, Kra71]. This isomorphism proves very useful for the study of quantum channels. If the 
physics of quantum channels are to be studied, instead simply the associated quantum states 
can be investigated. Then advantage can be taken of the fact that much is already known about 
quantum states. Their theory is well developed. 

Exploiting this isomorphism, several remarkable connections between maps and the associ- 
ated linear operators can be derived. For example consider the set of Hermiticity-preserving 
linear maps. Then the associated set of linear operators is the set of Hermitian operators. If 
the set of CP maps is considered, i.e. the set of quantum channels, the associated operators are 
the positive-semidefinite operators, i.e. all quantum states. An interesting link occurs when 
considering entanglement breaking maps. Then the associated linear operators are the separable 
quantum states. Furthermore, with aid of the isomorphism, relations between theorems on 
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states and theorems on channels and operations respectively can be derived. For example the 
operator-sum decomposition directly translates into the spectral decomposition of quantum 
states [AP04, HSR03]. 

Another very important application of the Choi-Jamiolkowski isomorphism can be found in 
entanglement evolution theory [KdMT + 08]. Consider a two-qubit state \x) AB = \pw \0) A |0} B + 
i/l — w \1) A |1} B subject to a one-sided channel p' = (t A ® Y B ) \x) AB (x\- Due to the isomor- 
phism this can be also expressed in the dual picture, where the roles of channel and state are 
interchanged: 

1(1* ®Y B ) \ X ) AB ( X \ = ±j(Y A ®l B )p AB . (4.65) 

Y A denotes the channel associated with the state |^)^ B and p AB is the state which corresponds 
to the channel Y B . p' = tr[(l* ® Y B ) \x) AB (x\] an d p" = tr[(Y A ® l B )p AB ] are normalization 
constants due to a possible trace decrease. This dual picture can be given a very descrip- 
tive interpretation using quantum teleportation, which is explained in detail and in a very 
demonstrative way in [KdMT + 08]. With this teleportation interpretation, Y A and p AB can be 
expressed in terms of the known variables |^)^ B and Y B . Finally, entanglement quantification 
is performed in the dual picture with aid of the concurrence and what comes out is a neat 
entanglement evolution equation for pure two-qubit states: 

C[(1 A Y B ) \ X ) AB (X\] = C[(1 A ® Y B ) |d>+) AB (*+|] C[\x) AB j. (4.66) 

As the concurrence is a convex measure, for mixed input states the evolution equation becomes 
an inequality, only providing upper bounds on the output entanglement. 

C[(l A ®Y B )p AB ] <C[(1 A ®Y B )\<S>+) AB (®+\}C[p AB ]. (4.67) 

In the derivation of this formula, the explicit form of the concurrence has been exploited. As this 
monotone is only defined for two-qubit states, a generalization to higher dimensional systems 
is not trivial. Nevertheless it exists, utilizing the G-concurrence by Gilad Gour [Gou05]. 
It is a generalized concurrence monotone for T-Lj ® finite-dimensional systems, which 
coincides with the "classical" concurrence for d = 2. Therewith, a generalized entanglement 
evolution equation can be derived in the same way, again making use of the Choi-Jamiolkowski 
isomorphism [TdMB08]. 

G d [(l A ® Y B ) \ Xd ) AB {Xd\} = G d [{t A ® Y B ) |d>+) AB <d,+ |] G d [\ Xd ) AB ], (4.68) 
G d [{t A ®Y B )p AB ] < G d [(t A ®Y B ) \^j) AB ($+|] G d [p AB ], (4.69) 

where \Xd) AB = ^ T!}j2o &ij 10 l/) B denotes a bipartite d x d-dimensional state with finite d 
and constants A lf £ C. p AB denotes some mixed state in the same Hilbert space. 

These entanglement evolution equations provide an elegant way for entanglement analysis of 
finite-dimensional states which are subject to a quantum channel. Unfortunately, they cannot 
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be employed for HE states, as these live in an overall infinite-dimensional Hilbert space. The 
question may arise whether it at least works for the class of effectively DV HE states. The 
problem is that the Gram-Schmidt process necessary for the DV description has to be performed 
after the channel. Admittedly it could theoretically also be performed prior to the action of the 
channel. However, in this case also the operators describing the channel have to be transformed 
and expressed in the new basis. In practice, this is often a rather cumbersome undertaking. 
Hence, for HE states entanglement evolution equations do not work very well as long as they 
are formulated for finite-dimensional systems. General CV entanglement evolution equations 
could not be derived yet. 



Two of the most important entanglement monotones for DV as well as for Gaussian CV 
quantum states are the negativity J\f and the related logarithmic negativity Ejv- An apparent 
question is whether a comparable entanglement evolution equation can be also derived for 
one of these measures. By a simple counterexample it can be shown that this is not the case. 
Consider again the state 



Apply a photon number encoding, where |0) corresponds to zero photons and |1) to one photon. 
The channel applied is the one-sided photon loss channel, which has been presented in the 
previous subsection. For qubits, its Kraus operators are given by [NCOO] 



\X) 



AB 




(4.70) 



£o = |0> (0|+v^|l) (1 



(4.71) 
(4.72) 



k x = yr^ |o) (i 



The output state p is given by 




(i A ®yf 0SS )\x) AB (x\ 




(4.73) 



+ ^V\0) A <1| ® |0) B <1| + (1-17) |1) A <1| ® |0) B <0| +n \1) A <1| ® |1) B <1|). 



Its partial transpose (PT) has one negative eigenvalue, i.e. 




(4.74) 



This yields the negativity 




(4.75) 
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It is also known that 

mx) AB ] = y/w(l-w), (4.76) 

^[(^®Yf 0SS )|d>+)^(d>+|] = |. (4.77) 

However, 

H{p AB ) £ ^w{l-iv)=H[{t A ®Y B 0SS ) |d>+) AB (®t\W[\x) AB \- (4.78) 

Hence, there does not exist an entanglement evolution equation in a comparable form to the 
one derived in [KdMT + 08]. Instead the output entanglement J\f(p AB ) is actually a complex 
function of N[\x) AB ] and Af[(l A (g> Yf oss ) |0^) AB and this for the rather approachable 

photon loss channel. It follows that £jv also offers no such entanglement evolution equation, as 
it is directly related to the negativity J\f by £^ = log 2 [1 + ZA/] . 



4.3.3. Logarithmic Negativity of a Nontrivial Hybrid Entangled State 

After this short excursion into entanglement evolution theory, HE states shall be discussed 
again. The first example of the class of mixed, but effectively DV HE states in the last but 
one subsection was the state (4.53), subject to photon loss. The output state only contained 
two qumode states and hence, after a Gram-Schmidt process, the state's entanglement was 
quantified with aid of the concurrence. However, the concurrence only works for two-qubit 
states. Therefore, another more representative state is to be discussed, which has to be analyzed 
with a more universal measure. The example is more artificially constructed, though. 
Consider the state 

p AB = p\ ( p + ) AB + (1-p) \<P-) AB (*-|, 
I a \ AB 1 fln\^l,,„„\ B i n I , ^ \B\ (4.79) 

which contains three qumode states \vac) B and |±a) B with «6lR (see figure 4.6 for a visualiza- 
tion). Perform a Gram-Schmidt process and obtain a qubit-qutrit entangled state in M, A ® TL B ■ 
The concurrence does not work for states of this kind anymore. The pure states in the convex 
combination of (4.79), after the inverse Gram-Schmidt process in an orthonormal basis, look 
like 

\<p+) AB = (|0) A |0) B + x \1) A |0) B + Vl - x 2 \1} A |1) B ), (4.80) 

\<P-} AB = (|0>^ |0) B + x\l) A |0) B -xVT^|l) A |1) B + y/l-x 2 -x* + x*\l) A |2) B ), 

(4.81) 

_ 1 I I 2 

with x = e 2 1*' . This time entanglement quantification is tried with the logarithmic negativity 



66 



4. Hybrid Entanglement 



ll> 

— \-a> 



1-p 



I0> 
Ivao — 



ll> 



\a> 



Figure 4.6.: Visualization of the HE state p AB , see (4.79), in phase space. The blue region 
corresponds to the first pure state in the convex combination with probability p, 
while the green region represents the pure state obtained with probability 1 — p. The 
additional |0) and |1) vectors denote the qubit states associated with the qumode 
states \vac) and |±a). 



Ejv[. As, for general bipartite mixed states with higher dimension than 2 x 2, it is basically the 
only actually calculable monotone, it is also the most important one for mixed, effectively DV 
HE states. The analytic expression for Ejs/(p AB ) is very long. It is omitted here. The graphical 
version can be observed in figure 4.7. It can be seen that for given a the lowest entanglement is 
obtained for maximal mixing p = 1. Furthermore, the higher a, the more entangled the state is. 
This is what has been expected. Finally, for p = V 1 and d^co the state approaches maximal 
entanglement. 

Of course, the discussion of this state does not bring much new insight. However, it is a good 
demonstration of the quantification of HE with aid of the important logarithmic negativity. 



4.3.4. A State containing N x d Qumode States 

To round off this subsection of mixed, but effectively DV HE states, a state which actually 
contains N x d qumode states shall be briefly discussed, namely the state 



\^) AB (fr| + (l-p) \<P2) AB <fc|, 



\<Pl) 



i AB 



AB 



\D A \-cc) 



\D A \-iec) 



(4.82) 
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Figure 4.7.: Logarithmic negativity of the effective qubit-qutrit HE state p , (4.79), as a function 
of probability p and amplitude a. For given cc maximal mixing p = \ yields the 
lowest entanglement. Furthermore, higher a also results in greater entanglement. 
For p = V 1 and «->co the state becomes maximally entangled. 



which contains Nxrf = 2x2 = 4 qumode states |±flt) and |±za) . Assume a £ IR. Figure 4.8 
provides a graphical illustration of the state. Of course again an inverse Gram-Schmidt process 
could be performed, which would yield a 8 x 8 density matrix in an effective ® Tif Hilbert 
space. The appropriate measure to be calculated would again be the negativity. Unfortunately, 
already for this 8x8 density matrix, the calculation of the eigenvalues of the PT state is not 
trivial. Since the negativity has already been calculated in the previous subsection and as the 
focus does not rest on heavy eigenvalue computations, this time witnessing with aid of the SV 
criteria is discussed instead. 

Try to witness entanglement with the well-known SV determinant S\. 
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(4.83) 



The graphical evaluation of this result is shown in figure 4.9. It shows that for small a ^ 
entanglement is verified, whereas for large a the determinant fails. For a = the state is 
separable and hence nothing can be detected. Furthermore, the more mixing occurs, i.e. p — > \, 
the harder to observe the entanglement. For p — V 1 the state is pure and the determinant 
witnesses its entanglement for any a 7^ 0. 

Concluding, it can be stated once again that the determinant Si is a great tool for the detection 



68 



4. Hybrid Entanglement 



\-a> 



ll> 



I0> 
\ia> 
1-p 



I— i&> 
ll> 



\a> 



I0> 



Figure 4.8.: Phase space illustration of the HE state p AB corresponding to equation (4.82). The 
blue and green regions correspond to the first and second pure states \(pi) AB and 
|<^>2} obtained with probability p and 1 — p. The |0) and |1) state vectors denote 
the qubit states affiliated to the qumode states. 



of HE, even though it has been effective DV HE in this case. With this statement the section is 
to be finished and finally the remarkable true hybrid entanglement is investigated. 



4.4. True Hybrid Entanglement 

Finally, the third class of bipartite HE states is considered, the truly hybrid entangled states. 
Recall their definition corresponding to equation (4.2) with N — 00: 



]C Vn \%)ab tyri] , p„>0Vn, Yj Vn = X 



n=l 
d-l 



\$n) = Yj Cnm \ m ) 



m=0 



d-l 

C n m 6 C, \Cnm\ 
m=0 



(4.84) 



1. 



As can be seen from these equations, truly HE states possess an infinite number of qumode 
states \ipnm) B ■ Therefore, the Gram-Schmidt process does not work anymore and the states have 
to stay in a Hilbert space of the form 7i A (g) H^,. Hence, they are not effectively DV, but instead 
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Probability p 



Figure 4.9.: The left diagram displays the SV determinant s\ evaluated for the state of equation 
(4.82). While for small « ^0 entanglement can be always detected, large a causes 
the determinant to fail. Note that the yellow plane denotes zero. On the right hand 
side the two regimes are plotted. For p = V 1 the state is pure and the determinant 
witnesses its entanglement for any a. 7^ 0. 



real combined DV (g) CV states and therefore truly hybrid. This is where the name actually 
comes from. Unfortunately this true "hybridness" has the effect that the states live in an overall 
infinite-dimensional Hilbert space in the non-Gaussian regime. This makes exact entanglement 
quantification impossible. Hence, the focus is only on the detection of true hybrid entanglement. 

In this section two examples for truly HE states are discussed, one more physically relevant 
and one artificially constructed state. Finally, general issues regarding entanglement quan- 
tification and witnessing in mixed states with an infinite number of pure state projectors are 
discussed. 

4.4.1. Noisy Qubit - Qumode Entanglement 

Once again, consider the state 

l^ B = ^(|0> A k) B + |l) A |- a > B ). (4.85) 

In subsection 4.3.1, this state was discussed when being transmitted through a one-sided 
amplitude damping channel of the form \ A (g) ^f oss ■ Now the channel is taken as not only lossy 
but also afflicted with thermal noise. Writing this out, 

p' AB = (l A ®Yf hermd )\ip) AB (ip\ (4.86) 

is the state to be investigated. The photon noise channel can be modelled in a similar way as 
the amplitude damping channel. The coupled environment is taken not simply in the vacuum 
state but in a thermal state (see figure 4.10 for a visualization). It is seen that the channel can 
again be written with aid of an ancilla Hilbert space (the environment E), and joint unitary 
evolution on the overall Hilbert space. Finally, a trace operation over the ancilla Hilbert space 
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Figure 4.10.: Visualization of the photon noise channel. The environment mode is in a thermal 
state and coupled via a beam splitter to the input state to be transmitted. Subse- 
quently the environment mode is traced out and the decohered output state is 
obtained. 



has to be performed. Hence, 

?** = tr £ [!i BE (|^ B mpl^J) U BE \ (4.87) 

with 

• beam splitter unitary U BE = e e (4 a B- a B 3 E), 

• environment thermal state pf hmml = E„=0 (l^S,))^ 1 ^ E ^ [SZ97 ^ 

• mean thermal photon number (n t h ), 

• and beam splitter transmissivity rj = cos 2 9. 

The two main parameters characterizing the channel are the mean thermal photon number 
(nth,) an d m e beam splitter transmissivity rj. For (ttf/,,) = 0, the photon loss channel is obtained, 
which is therefore just a limiting case of this channel. 
For the output, 

p' AB = (t A ®Y? hermal )\rp) AB W 

= «ElU BE (\cp) AB (cp\®pL. mal ) U BE '] 
-i oo n / 

= ~ LPn- E fnk{n)Un)(KUn) |o) A <o| ® « + |^) B (W«I^ 

n=0 k,l=0 V 

+ ^" a ( J /) ii> A ai «> a" 1 * |-v^7«> B <-V^*l ^ 

+ ^7('7) |0) A (1| ® « + "|^> B (-^«|a' 



(4.88) 
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is obtained, a and fl + are the mode operators of system B. p„ ' denotes the thermal photon 
distribution . , and /„)t(?7) and A^(t/) are defined as 



:= ^= (Jj^f-A/R)*, (4-89) 

:= (V^V^^'lV^l*) . (4.90) 
For the calculation, the following relations have been exploited. 

U B£f U BE — 1, (4.91) 

U BE {a Ei ) n U BEf = {^ja Ei - ^T^na B *) n , (4.92) 

U BE D B (a)U BEf = D B ( v ^a)D E ( v /l -not), (4.93) 

Q BE |0) B |0) £ = |0) B |0) E . (4.94) 



It can be infered from equation (4.88) that the state p' AB clearly is truly HE, as it contains 
an infinite number of qumode states {a* k \±^/rja) B : k = 0, 1, ...,oo}. Hence, true hybrid 
qubit-qumode entanglement is obtained. 

Furthermore, writing the state in this form, (4.88), illustrates the effects of the thermal 
photon noise channel in a very concrete way. On the one hand, the damping effect due to 
the beam splitter is clearly visible in form of ^/fj in the states a \±^ffjOL). On the other hand, 
there is not only damping, but also thermal photon noise, which becomes manifest in the 
creation operators a in the states a | ± y^a) • Thermal photons "leak into the system" and 
are "created" in the damped coherent states. Finally, it is clearly visible, how each term of 
Pthermal = ^m=0 P<n' I w )e ( n \ results in the creation of maximally n noise photons in the coherent 
states. However, as descriptive (4.88) is, as inapplicable it is for further calculations. To perform 
entanglement witnessing in this state, the moments required for the SV determinant Si have to 
be computed. Unfortunately, when this is tried, making use of the state as written in equation 
(4.88), intractable infinite sums are obtained whose convergence behavior is impossible to be 
worked out exactly. Of course, truncation at some very large n could be performed, which 
would probably result in very accurate outcomes. However, such a procedure is opposed to 
the actual intention of analyzing true HE, since a truncated state is not truly hybrid entangled 
anymore. This is a point which makes the investigation of true HE particularly challenging. 
Infinite sums or integrals emerge, which have to be calculated exactly. 

At least in this case, there is a way out of the dilemma. Instead of utilizing pf hermaj = 
X^n°=0 Pn I m )e ( n \ m me Fock basis, make use of the coherent state basis. Sufficiently classical 
states p can be written as p = f c d 2 a.Pp(a,a*) \ ot) (a|, where Pp{m, oc*) denotes the Glauber- 
Sudarshan P-representation of p. The thermal state is of course sufficiently classical and hence 
its P-representation exists [SZ97]: 
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Therefore, 



Pthermal _ / \ / 



1 /" o H 2 

— r / d 2 ae W |a)(a|. (4.96) 

7T (nth) Jc 



Expressing the thermal state in this form, the action of the thermal channel on an element 
| a) (jS| looks like 



Y«(k> B (PI) = ^ytr £ [LT BE (| a > B </3| ® jf d 2 7 e"<4) | 7 > E ( 7 |) & 
= / c d2 Te"^tr E [!J BE (|«} B </5| ® | 7 > E < 7 |) fj 

Define £ := U BE (Wf (P\ ® |t) E ( 7 |) fi BE+ and calculate 

X = !i BE (k) B (/5|®| 7 > E ( 7 |)!J BE+ 
= !i BE D B (a)D E ( 7 ) f |0) B (0| ® |0) E (0|) D E \-y)t) B \p)U BE ' 



(4.97) 



= (tJ BE D B (a)tJ BE+ ) (!i BE D E ( 7 )!i BE+ ) !i BE (jO) B (0| 

® |0) E (0|) !i BE+ (!i BE D E+ ( 7 ) !i BE+ ) (!i BE D B+ (/3)LJ BE+ ) (4.98) 
= D B (^ a )D E (yr^ a )D E (^ 7 )D B (-yi^ 7 ) (jO) B <0| 

® |0) E (0|)D Bt (-yi^ 7 )D E+ (^ 7 )D Et (y^)D Bt (^/5) 

= Iv 7 ^" V 1 "77> (v 7 ^ " V 1 ~ 77 1 ® I V 71 " 7" + W7> (Vl- ?/jS + ^ 7 |, 

where once again the relations (4.91), (4.93), (4.94) have been exploited. Additionally use has 
been made of !i BE D E ( 7 )ii BE = D E (y^ 7 )D B ( — y/l — tjy). The derived expression for % yields 

g 

• |V7«- V 1 - 77) (^7/3 - V 1 - 77l • 
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Making use of this auxiliary calculation, the overall state p' AB = (t A <g> Y B hermal ) \ty) AB (Vl with 
W) AB = ^ 2 (\0) A \a) B + \l) A \-a) B ) is given by 

ff M = (l A <SVl mua )\rp) AB (1>\ 
l r o - 17,2 / 

|0>^ (0| ® - yi^/7> B (^7<* " v/l^Tl 
+ |1) A (1| ® |-^7« - 71 - J/T) 5 <"^?a - v 71 " 77 1 

+ A a -»T(7/) |0) A (1| ® | - " v/l^Tl 

+ A-"T(^) |1} A <0| ® |-v^« - yi^7> B (v^« " V^rrl), 



(4.100) 



with 



= exp 

+ ( V^VP* + W7* ) ( Vi 3 ^* + Vvi) 



IW^lP + Vnif-lW^VK + Vrill 2 (4.101) 



This form of the state p' AB , (4.100), is not as insightful as (4.88). However, it is mathematically 
much more controllable. Instead of infinite sums, in this form just integrals which can be easily 
calculated occur. As an example, consider the moment (a^ab^b) (actually a and a f are the mode 
operators corresponding to system B and b and fr + those of system A.): 

If _ M 2 _ 

(a f ab f b) = - — y J d 2 je {n th) (—y/tjcc — y/l — nj\a f a\ — y/rja — y/l — nj) 

= 2n(n th ) Ic d2 ^ e ~^ (v\ a \ 2 + 0-- l)\l\ 2 + yfilQ -V)(*l* + «*7) 

= inku) LL dxdy ^ + (1 " » ){X2 + ^ 

+ \Jn{l-n) (oc{x - iy) + a*(x + iy)) 



(4.102) 



r/\oc\ 2 1 — t] 
2 

|2 



+ p 7 / dx x 2 e > 



n\oc\ z 1 — rj 



{n t h), 

where 7 = x + iy has been substituted. Additionally, some integral identities have been 
exploited, which are presented in appendix B.7. 
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The calculation of the other moments in the SV determinant Sj proceeds similarly and finally, 
s tplAB (x, V ,(n th .)) = (n th ) (l - !^L) - ^e- 4 l"l 2 (4.103) 
is obtained. A graphical illustration for t] = | is shown in figure 4.11. It can be observed 




Figure 4.11.: The left hand graph displays the SV determinant s\ IAB (oc,rj = |, (nth.)) ror the 

state j6 /ylB with n = % > . Without loss of generality a £ R has been assumed. 
On the right hand side the the two regimes are plotted. There can be clearly 
identified a region, where the determinant is below zero. Hence, true HE can be 
witnessed. Furthermore, there is again a trade-off behavior, which results in the 
existence of an optimal a opt . This oc opt corresponds to the most robust state p'^,t 
regarding entanglement witnessing for fixed rj = | and varying mean thermal 
photon number (n t j L ). 



that there clearly is a parameter region in which entanglement can be detected. Also note 
the trade-off behavior and the optimal oc ^ which corresponds to the most robust state p 
regarding entanglement witnessing for fixed rj — | and varying mean thermal photon number 
(n t h ). Furthermore, it is worth pointing out that the witnessed entanglement in this case is 
actually true HE, in contrast to the entanglement witnessed in previous sections. It can be 
concluded that the SV determinants provide a suitable tool for the detection of true HE. 
When setting si_ ;Ag (a, rj, (n^ )) < and solving this equation for (n t h.), 



p 

12 



(n th ) < (4-104) 

N ' (1 — ^)(2e 4 l^l 2 — 1) 

is obtained. For parameters (oc,n, (n^ )) satisfying this inequality, entanglement is detected. 
Furthermore, the inequality can be used to define a surface. The parameters (a, rj, (n t h)) en t for 
which entanglement is verified lie below this surface (see figure 4.12). However, it is rather 
cumbersome to read off exact parameters in such a 3D plot. Hence, regions of successful 
entanglement detection are plotted in figure 4.13 for different values of rj. As expected, 
the higher the transmissivity rj the greater the parameter regions of entanglement detection. 
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Amplitude a 



Figure 4.12.: Surface defined by equation (4.104) (a £ IR). For parameter triples (a, n, (n^ )) 
lying below it, entanglement is witnessed. Once again, the trade-off behavior and 
an optimal a. ?* for which entanglement can be detected for the strongest possible 
noise can be recognized. 



Furthermore, the plot shows that the optimal oc ^ seems not to depend on rj. Actually, it is 
straightforward to calculate a P t . Simply solve 

d 4??|a|2 -0 (4.105) 



doc (l-^)(2e 4 l fl; l 2 -l) 



for a. This yields a opt ?s 0.44, independent of r\. Furthermore, also a similar equation to (4.104) 
can be written down when solving s\ f .„ {ot,rj, (n t h)) < for rj: 



p. 



+ ( 2e *l«l 2 -l) 

Calculating ol ^ with aid of this equation yields of course, as expected, the same result 
a°P f m 0.44, independent of (Wf/,.). It is actually quite interesting that the optimal amplitude 
K opt regarding entanglement witnessing with the SV determinant sj does not depend on the 
channel at all. a ^ seems to be only determined by the choice of the SV determinant. 

Compare this to figure 4.5, which shows the concurrence of the state when setting (Wf/ L ) = 
0. There, the optimal a°P f does very well depend on the transmissivity rj: The greater the 
transmissivity, the higher a°P' . This is quite remarkable. It can be infered that the ability to 
detect entanglement in p' AB , depending only on the choice of the SV determinant, does not 
scale as the entanglement of p' AB itself, which of course depends on the thermal channel's 
parameters rj and (n^ ). 
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Figure 4.13.: Regions of entanglement detection for a few different values of transmissivity rj. It 
is interesting to note that the optimal a°P f seems to have a fixed value (gain a 6 R). 

Concluding, a first example of true HE has been presented and entanglement detection has 
been performed. Some interesting properties regarding the comparison between entanglement 
witnessing and entanglement quantification have been identified. Furthermore, the exemplary 
state is additionally also physically relevant, as it occurs in quantum information applications 
such as QKD. 

Finally, it is worth to mention that the calculations regarding the thermal photon noise channel 
have led to the derivation of a set of corresponding Kraus operators. They are presented in 
appendix B.6. 

4.4.2. Another Truly Hybrid Entangled State 

Another example for a truly HE state which admits entanglement detection shall be presented. 
However, this state is rather artificially constructed and probably not of great physical relevance. 
Nevertheless, it is truly HE and enables nice HE witnessing: 



n=l 




(4.107) 



p n = x n , < X < 1, HER. 



As this state contains an infinite set of qumode states { | ± y/na) : n — 1, 2, . . . , oo}, it is clearly 
truly HE, strictly speaking it shows true HE between a qumode and a qubit. 

Now, we try to witness entanglement in this state with aid of the SV determinant S\. For the 
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calculation of the moments, the identities 

n=l 1 X 

00 r 



(4.108) 
(4.109) 



n=l 



which are derived in appendix B.8, are exploited. As an example, the calculation of the moment 
(a^ab^b) is presented. Note that the mode operators a and a + correspond to system B, while b 
and P denote the mode operators of system A. 



net 



n=l 



a 1 1— x 
~2 x 



£>x" 



n=l 



a 2 1 
Yl-x' 



(4.110) 



For the whole determinant S\, 



si(x,a) = 



2a 1 
1-x 

a(l — x 



J- xe n=l " n=l 



V X n=l 7 V X n=\ 1 



n=l 

a 2 / e- 2a2 (l-x) 
1-iV 1 - xe" 2 * 2 



(4.111) 

is obtained. Note that this expression does not depend on the phase (p. This has been expected, 
since the amount of entanglement in states of this form does not depend on (p, as pointed out 
earlier. 

It is clear that 



X>" < J^Vnx" < J^nx", 

n=l n=l n=l 



for < x < 1. Hence, 



xe 
1-ie- 2 " 2 



T^<£^"<(T-^ 

~ 2fl;2 00 „ ? re 

<^(ie- 2t )"< 



«=1 



2a 2 

(l-xe- 2 « 2 ) 2 ' 



(4.112) 

(4.113) 
(4.114) 



Inserting the lower bounds into the sums of equation (4.111) yields 



s[ (x, a) 



8 



2x / 
1-x ~ \\ 



1-x 



xe 



-let 2 



-4^ 



1-X 



(4.115) 



However, since s\(x,a) > s'^{x, a.) Va e R, x s]0, 1[, from si(;c,a) < follows s[(x r ix.) < 0. 
Therefore, si {x, ot) < is a sufficient criterion for entanglement detection. It is plotted in figure 
4.14. The graphics show that entanglement can be verified for sufficiently small x and a. This 
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Fi gure 4.14.: On the left hand side Sj(x, a) is shown, which is derived from the original SV 

2 

determinant S\(x,0i). Note that for the plot the prefactor ^- has been omitted, 
which, of course, does not affect the witnessing region. The yellow plane denotes 
zero. The right hand side graphic presents the cutting line between s'^x, a) and 
the zero plane. It displays the witnessing region. It also has to be pointed out that 
for a = the state is of course not entangled. 



can be understood in the following way. Small x corresponds to low mixing between the pure 
states \xp n ) AB , which themselves are rather highly entangled, depending on the amplitude ^Jna. 
For pure states of the form \xp n ) AB entanglement witnessing can be performed perfectly with 
the SV determinant Sj. Therefore, efficient entanglement detection for low x is no surprise. 
Furthermore, already previous calculations (compare subsections 2.4.5 and 4.2.2) have exposed 
that entanglement detection via Si works best for low a. and may fail for high a, even if the 
entanglement itself may be higher. 

Concluding, a second example for a truly HE state has been presented, whose true HE could 
again be verified with aid of the SV criteria. For completion, note that states of the form 




(4.116) 



can be generated via the interaction G, 

G := e2\ a \ 2{1 - n \a z ^if\ (4.117) 
from separable product states of the form 

l^ = ^(|0> A +e^|l) A )k) B . (4.118) 

Once again, a and a + are the mode operators corresponding to subsystem B, while cr z is the 
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Pauli Z matrix acting on subsystem A. Therefore, 

G |fr>5 = e^l 2 (i-«) (Oi Vn) rt ^ (|0) ^ + e* |«) B 

= -^(|0> A \^a) B + e i <P \1) A \-Vn-af) ( 4 ' 119 ) 
= \%) A \ 

which becomes obvious when expressing the coherent state \a) in the Fock basis. 
4.4.3. Truncated Mixed States and Entanglement Witnessing 

This subsection does not really deal with truly HE states, but with general entangled quantum 
states of the form 

oo 

P=Y i Vn\lpn)(lpn\, Pn ^ OVn. 5 (4.120) 

n=0 

The state may be bipartite or multipartite entangled as well as CV, DV or hybrid. Hence, truly 
HE states are included (note that effectively DV HE states are not included.). The focus here 
lies on the infinite number of terms. In the previous subsections, it emerged that one of the 
main problems with true HE is that the occurring infinite sums or infinite integrals have to be 
calculated without truncation. A truncation would correspond to analysis of a different state, 
which would not have the important infinite number of mix terms character anymore. 

Imagine for example the following case. Consider a separable state E~=o Pn \$n) {tf>n\ for 
which, however, the truncated state = E^=o Pn \tfn) {tyn\ is entangled for some N. Therefore, 
a truncation in occurring sums or integrals may lead to misinterpretations. Entanglement in 
would potentially be detected and wrongly misinterpreted as entanglement of the state with 
N = oo. 

The questions arising are the following: In which cases is, for an entangled state of the form 

oo 

P=Y,Vntyn){lpn\, (4-121) 
n=0 

also the state 

J N N 

PN = — ^Pn\lpn) (lpn\ , Pn ■= E P«- (4.122) 

P N n=0 n=0 

entangled? The prefactor ^ is due to normalization of the truncated state. Or, coming to the 
important point: When successfully witnessing entanglement in 

Pn = -£p»W {ipn\, (4-123) 
Pn „=o 



5 A1so states with some p„ — could be considered as long as the total number of p„ is still infinite. However, 
when omitting the p„ — terms in this case, one ends up exactly at the form (4.120). 
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in which cases also entanglement in 

oo 

P= LPn\lpn){lpn\ (4.124) 

can be deduced? An answer to this question is provided in the following. 

Theorem 4.3. For p = E~=oP« I'M with Vn 7^ OVn define p N = ^E£LoP« I'M (%\ f or 
some N 6 No wz'f/z Pn := E^=o Vn- Then the following statement holds true: IfRg(pN) > 1 — pw, p 
is entangled. 

Proof. 

• ^g(PN) > 1 - PN- 

• => PnPn + (1 — Pn)^ entangled for any quantum state &. 

• Choose CT = En=N+l Vn \^n) (%\- 

• PnPn + (1 - Pn)^ = P entangled. □ 

The answer to the above question is therefore closely related to the global robustness Rg. 
Unfortunately the calculation of R g is mostly quite difficult or even impossible. Nevertheless, it 
has to be calculated for a state with only finite N, which may be possible in certain circumstances. 
If entanglement detection has to be performed for a state with N = oo yielding infinite sums or 
integrals which cannot be calculated without truncation this theorem may be a useful tool if 
robustnesses can be calculated for the truncated states. However, it is not clear which N has to 
be chosen. In the case that the state is actually not entangled, such an N with R g (p^) > 1 — p^ 
does not exist. Proving this is trivial: Assume that there exists such an N for a separable state. 
Then this would indicate that the separable state is actually entangled, which is a contradiction 
to the assumption. 

Since there is an infinite amount of possible Ns, no statement about the state's entanglement 
can be made, if no such N is found. Working through "all" Ns is impossible, since there is no 
finite "all". Nevertheless, for any entangled state of the form 

oo 

P= J2Vn\lpn) {lpn\ , p„^0VM, (4.125) 
n=0 

such an N does actually exist. 

To prove this a new entanglement monotone is required, which is introduced now. The 
entanglement of subtraction E Su h(p) is defined as 

Esufc(p) := inf {b :b>0,Ae QS, O e SEP such that p - bA = (1 - b)Ci}, (4.126) 

where QS denotes the set of all quantum states and SEP the set of the separable quantum states. 
The interpretation of this monotone is clear. It asks how much of a quantum state must be 
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subtracted from the state in question that it becomes separable. However, for Eg u f,(p) to be a 
proper entanglement monotone, first its LOCC monotonia ty has to be proved [PV06]. 

Consider an entangled state p which is supposed to have Eswb(p) = R- Then, its optimal 
decomposition according to the entanglement of subtraction is given by 

p = aCl + RA, (4.127) 

with a = 1 — R. We perform a general LOCC operation on p, which results in the output 
pi = with probability <]{. The output ensemble is given as 

. flA,-(Q) +RA,-(A) _ , A f (n) p A f - (A) 
{q ' ' } = {<7i ' fli tr[AKn)] +Rf tr[A7(A)I } ' ( > 



where 



dtr[Aj(n) 



<7 



(4.129) 



^ = *tr[A,-(A)] ( 
<?i 

Since a, R > 0, and therefore also 5;, > 0, equation (4.128) is a valid decomposition of 
{cji ; A,(jo)} of the form ; 5,0 + RjA}. However, it is not necessarily the optimal one. Denote 
the actual entanglement of subtraction of the output A ; -(j9) with R{Aj{p)). Then the average 
output entanglement is given by YU a iR{^i(P))i- Finally, 

^iR(Ai{p))i < Lli R i = R J2 tl i A i( A )] = R - (4.131) 

i i i 

Hence, the entanglement of subtraction does not increase on average under LOCC. Furthermore, 
it is clear from the construction that Eg u i, is zero if and only if the state is separable. Therefore, 
it is an entanglement monotone. Some properties of Eg„j, shall be listed: 

• < E Sub < I. 

• p separable E Sub = 0. 



For b = 1, every state "becomes" separable: V p, choose A — p and hence, p — A — O E SEP, 
where O denotes the zero operator, which is of course separable. 

Two annotations have to be made regarding this new monotone: The subtraction of quantum 
states from other quantum states does not necessarily result in an output which is a valid 
quantum state [OttOl]. However, in the definition of Esufc the variation is performed in such a 
way that only those subtractions are considered which yield a valid quantum state as output. 
Furthermore, it has to be pointed out that the entanglement of subtraction is a special monotone 
of a whole family of monotones, as presented in [PV06]: 
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Consider two sets of operators X and Y which are convex cones and closed under LOCC 
operations. Furthermore, the members of the sets can be written in the form u X (y) -positive- 
semidefinite operator with trace oc X (yy a x(Y) are fixed real constants. Finally, the sets have to 
be chosen such that any Hermitian operator h can be written as h = atl — bA, where O E X 
and A e Y and a, b > 0. Then a monotone may be defined as 

R x ,y(P) ■= mf {b : a, b > 0, A 6 Y, Cl £ X such that p = aCi - bA}. (4.132) 

For oc x = 1 and ay = — 1, the entanglement of subtraction is obtained. 

Finally, recall theorem 4.3. With the aid of the entanglement of subtraction it can be proved 
that for every entangled state p of this form actually an N exists such that Rg(p^) > 1 — p^. 

Theorem 4.4. For any entangled state of the form p = E~=o Pn IVn) (^nl w ^ Pn 7^ OVn there exists 
an N 6 N such that R g (p N ) > 1 - p N for p N = ± E~=o Pn I'M (ipn\ with p N := E~=o Pn- 

Proof. 

• p entangled => E sub (p) > x > 0. 

• => jSi = i=f (|6 — tpl) entangled for < t < E su j,(p) and arbitrary states pi. 

• Furthermore, E~=o P» = 1 => ]iai n - yoo p n = Ve > 03N e N such that 

E~=N+1 Pn < e - 

• Choose N such that E~=n+i Pn = 1 _ Pn < x < E sub (p). 

• Set £ = 1 - p M . 

• => jSi = — (1 — Pm)Pi) entangled for any pi and M > N. 

• Choose p 2 = E~=m+i Pn I'M 

• =>• pi = Pm = E^lo P« I'M (ifal entangled VM > N. 

• Consider the global robustness of pM- Define R g (p M ) := x M . As pM entangled VM > N, 
% > VM > N. The important point is that Xm converges definitely not against zero: 
Due to poo — p it is known that Xoo — Rg(p) > 0, as the initial state p is entangled. 

Concluding, on the one hand, lim M ^ooX M = Rg{p) > 0, while, on the other hand, 
lim M _ >co (l - p M ) = 0. 

• => 3M> N such that R g (p M ) >l-p M . 

• ^ Pm = Zn= Pn \ipn) (fn I entangled with Rgip^) > 1 ~Pm- d 

The consequence of the theorems 4.3 and 4.4 is that entanglement witnessing problems 
regarding states of the form 

CO 

P=EPn\lpn){lpn\, Pn^0V«, (4.133) 
n=0 
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can be reduced to entanglement quantification problems of states of the form 

N 

PN=YLVn\^n){^>n\, (4.134) 
n=0 

with finite N, where the employed entanglement monotone has to be the global robustness. 
However, the applicability of this approach is quite limited. 

• No recipe is given for finding the right N. Simply several choices of N have to be tried. 
Obviously there is an infinite number of possible choices. Therefore, it is more or less a 
matter of luck to find an N with R g (p^) > 1 — p^. 

• Calculating the global robustness for general states px is an extraordinarily difficult task, 
since it corresponds to a complex variational problem. Additionally, another entanglement 
measure or necessary separability criterion is required to check whether a state (1 — A)j6jv + 
A<7 actually is separable or not. Powerful necessary separability criteria do not even exist. 
However, this problem may be avoided by loosening the constraints and requiring for 
example only PPTness instead of separability for the robustness calculation. 

• Entanglement witnessing is relevant in cases in which it is not known whether a state 
is entangled or not. But in the case that the state is not entangled at all, no N with 
Rg(p^) > 1 — p;v can be found. Hence, it is in the first place unknown whether such an 
N does exist or not. An extensive search for such an N is therefore pointless, as it may 
not exist at all. 

Nevertheless, the approach is not fully senseless. In cases, where for the original N — oo state 
no calculations can be made without truncating the occurring infinite sums or integrals, it is the 
only chance to actually witness entanglement with certainty. Any such truncation during the 
witnessing process for the original N — oo state p would correspond to witnessing entanglement 
in an actually different state. Then, if robustness calculations are possible for the corresponding 
truncated finite N states pN, being lucky and finding an N with Rg(ppj) > 1 — p?j is the only 
chance to unerringly witness entanglement in p. 

4.5. Multipartite Hybrid Entanglement 

In this section some investigations regarding multipartite hybrid entanglement are presented. 
In general a HE quantum system may of course not only be built of two subsystems of finite 
and infinite dimensionality, but also of several finite- and infinite-dimensional subsystems. 
In general, stepping from the bipartite into the multipartite regime yields a wealth of new 
remarkable effects as well as difficulties. These shall not be discussed here in detail. 6 However, 
some interesting phenomena regarding multipartite HE are examined. In the first subsection, 

6 Brief introductions to multipartite entanglement are given in [PV06] as well as in [HHHH09], which is one of the 
most thorough reference works regarding quantum entanglement. 
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some fundamental concepts are briefly introduced which are necessary for further studies on 
multipartite HE. Afterwards, general multipartite HE issues are discussed, and two specific 
states are investigated. 

4.5.1. Multipartite Entanglement Basics 

Several applications in quantum information require the consideration of entanglement between 
more than two parties. For example, when performing multi-user quantum communication 
between more participants than just one sender and one receiver, multipartite entangled 
states have to be established [MJPV99, MPVOO, KKI99]. Furthermore, the already mentioned 
cluster states required for certain quantum computation schemes are multipartite entangled 
states of a very large number of particles [RB01, Nie06]. Additionally, the already mentioned 
qubus approaches for quantum communication and quantum computation actually involve 
multipartite HE states. The CV qubus entangles itself consecutively with two qubits, resulting 
in a tripartite HE state, before it is measured out such that the remaining two-qubit state is left 
entangled [vLLMN08, BPvLlO]. 

One of the most crucial differences between bipartite and multipartite entanglement becomes 
evident when trying to define maximally entangled states from which all other states can 
be obtained via LOCC with certainty (for the sake of clearness and simplicity consider two- 
dimensional qubit subsystems in this introduction). Motivated from the bipartite setting, an 
obvious possibility would be the Greenberger-Horne-Zeilinger (GHZ) state \ GHZ), which looks 
in a tripartite setting like 



It has some very neat properties. The entanglement of any bipartite cut, for example the 
entanglement of party A vs. parties B and C together, takes one ebit. Furthermore, when a 
projection measurement on ^|(|0) + |1}) is performed, always a maximally entangled two-qubit 
Bell state is left over. Therefore, any bipartite state can be obtained from it with the aid of LOCC 
transformations. However, there are tripartite states which cannot be generated via the GHZ 
state and LOCC, as for instance the W-state: 



This impossibility of interconversion of the GHZ- and the W-state holds also when switching to 
an asymptotic setting and considering arbitrarily many identical copies of the states [PV06]. 7 
Therefore, none of the states can be considered as more entangled than the other. Actually, 
there are different kinds of multipartite entanglement. On the one side, a tripartite state may 

7 In the bipartite setting, the consistent picture of entanglement quantification via LOCC transformability is actually 
not obtained until, on the one hand, transformations between arbitrarily many copies of the states are considered, 
and, on the other hand, asymptotic imperfections are allowed, which vanish when the number of copies tends to 
infinity [PV06]. 




(4.135) 




(4.136) 
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just be bipartite entangled between its subsystems, while, on the other side, it may also show 
"genuine tripartite entanglement". The generalization to higher dimensional systems proceeds 
accordingly. 

This motivates to go beyond a simple distinguishment between separable and entangled 
states. For a N-partite state full N-separability is still defined in such a way that the state can be 
written as 

pABC... = £ v .pA gj pB g, pC g, _ _ _ (4.137) 

i 

However, there is also the notion of partial separability when only some subsystems separate 
from the remaining state, which may be inseparable. Furthermore, an N-partite state which 
shows entanglement between maximally k subsystems is said to be k-entangled. A tripartite 
state is for example said to be 2-entangled if it can be written as 

P ABC = E ViPt ® Pf c + X>t ® Pt c + X>f ® Pt B - (4-138) 

/ i 

However, when taking many copies of states with, for instance, k-entanglement, also 1-entangled 
states with / > k may be created via LOCC [PV06]. Therefore, this notion of "k-entanglement" is 
either not closed regarding LOCC transformations or under taking many copies of the states. 



Three-particle systems: 
Potential entanglement distributions 




A) 2-entangled B) 3-entangled C) 3-entangled 

W-like entanglement GHZ-like entanglement 



Figure 4.15.: Three possibilities for entanglement distribution in a pure three-particle system. A) 
corresponds to a state of the form |0) ® (O) , where jO) denotes an entangled 
two-qubit state. B) is the visualization of the W-state \W) ABC , while picture C) 
presents a GHZ-state \GHZ) ABC . It becomes clear that for a system in a W-state 
losing one particle leaves the remaining particles entangled, while for systems 
prepared in GHZ-states particle loss destroys all entanglement. 
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This does actually not solve the problem with the tripartite GHZ- and the W- state, which 
are both 3-entangled. When tracing out one party of the W-state, always an entangled state is 
left over, while in case of the GHZ-state a trace over one subsystem yields a separable state. 
Therefore, the GHZ-state seems to be genuinely tripartite entangled, showing no bipartite 
entanglement, while the W-state seems to offer only bipartite entanglement between its subsys- 
tems. Nevertheless, both states are 3-entangled. Figure 4.15 visualizes some possible scenarios 
of entanglement in pure tripartite states. This discussion demonstrates quite well that the 
characterization of entanglement in multipartite settings is much more involved than in the 
bipartite regime. Actually, there are a lot of open questions, especially when considering not 
only tripartite states but systems of higher order containing much more parties. 

Interestingly, there is a way to quantify how much GHZ-like, genuine tripartite entanglement 
a tripartite state contains. Coffman, Kundu and Wootters (CKW) have shown that 

C 2 (A\B) + C 2 (A\C) < C 2 {A\BC) (4.139) 

holds true for any tripartite state [CKWOO]. C 2 (A\B) and C 2 (A\C) respectively represent the 
entanglement between parties A and B as well as A and C, measured by the squared concurrence. 
For the calculation of C 2 (A\B) subsystem C has to be traced out and the ordinary bipartite 
concurrence has to be applied, analogously for C 2 (A\C). Finally, C 2 (A\BC) denotes the squared 
concurrence of the bipartite cut between party A and parties B and C. CKW could show that 

C 2 {A\BC) = 4det[,o A ] = 4det[tr BC [/ BC ]]. (4.140) 

The relation (4.139) is also known as the mongamy of entanglement, which says that if two systems 
are maximally entangled, they cannot be entangled with another third system. This becomes 
quite clear from the form of the inequality: If A and B are maximally entangled, C 2 (A\B) = 1. 
As all quantities in the formula are between one and zero, C 2 (A\B) = 1 necessarily yields 
C 2 (A\BC) = 1 and C 2 (A|C) = 0. Hence, system A cannot be entangled with system C. A 
similar equation holds also for system B and therefore also system B cannot be entangled with 
another system. Note that there is a generalization to more parties: 

C 2 (A\B)+C 2 {A\C) + C 2 (A\D) + ... < C 2 {A\BCD . . .). (4.141) 

However, stick to the tripartite case and consider (4.139). For the "=" case the overall en- 
tanglement between party A and parties B and C C 2 (A\BC) is just bipartite entanglement, 
namely C 2 (A\B) + C 2 (A\C). However, when the overall entanglement C 2 (A\BC) is greater 
than the amount of bipartite entanglement C 2 (A\B) + C 2 (A\C), of what kind is the remaining 
entanglement? Hence, define the so-called residual entanglement x res , 

r res :=C 2 {A\BC)-C 2 (A\B)-C 2 (A\C). (4.142) 

In the tripartite case, it corresponds to the GHZ-like genuine tripartite entanglement, which has 
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already been mentioned regarding the GHZ-state. And actually, as expected, 

r^ z = 1, (4.143) 
r% = 0, (4.144) 

is found for the GHZ- and the W-state. This shows very nicely and actually quantifies the 
characteristic difference between these two kinds of states. As T res < 1 in general, the GHZ 
state is in the following considered as maximally entangled with regard to its special kind of 
GHZ-like genuine tripartite entanglement. The W-state is of course also 3-entangled, however, 
its 3-entanglement is based on 2-entanglement between its subsystems, which can be inferred 
from the fact that it satisfies the CKW inequality (4.139) for the "—" case (see also figure 
4.15). Therefore, W-like 3-entanglement is not regarded as "genuine tripartite entanglement". 
Nevertheless, it is a special form of 3-entanglement. 

Note that r res does not change, when switching the roles of the systems A, B, and C. In the 
multipartite case with more than three parties, the analogously defined residual entanglement 
corresponds to the amount of entanglement in the state which is not bipartite or based on 
bipartite entanglement. 

In the following section, this important work by CKW is applied on tripartite HE states. 



4.5.2. General Multipartite Hybrid Entanglement 

In general, N-partite HE states live in Hilbert spaces of the form <8> . . . <g> Wd N , where some 
dj are finite and some infinite. As this chapter is primarily supposed to be an outlook, and for 
the sake of simplicity, mainly tripartite HE states are considered. Then there are two cases: 
Either the Hilbert space looks like <g> T-L^ <g> Hoo or % d <g> Hoo <§D Hoo, with finite d\, di- 
Consider the first case, where only one subsystem is CV. A general state in Tif €5 ® 
can be written as 

]ip) ABC = V) A |y) B ]ipij) C f ^ |c .. |2 = ^ (4 145) 

i,/=l »',/' 

where \ipij) represents some qumode state. As there are at most only d\ x di such qumode 
states, a Gram-Schmidt process can be executed to write the state as a multipartite, effectively 
DV HE state. Therefore, if all \rpij) C are linearly independent the state effectively lives in a 
Hilbert space of the form %j <g> U\ ® xd . 

Now move on to the second case, where the initial Hilbert space looks like Hf ® ® "H£. 
Then the most general states are 

W) ABC = E c vl«> A l^> B l^> C ' ( 4 - 146 ) 

where both \<pij) B and \tyij) C denote qumode states, and the Cy are chosen such that (t/>|i/>) = 1. 
Obviously, in this case an infinite number of qumode states is present for an infinite number of 
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Cjj 7^ 0. No Gram-Schmidt process can be performed. Hence, these states show multipartite true 
HE. 

Concluding, in the multipartite regime already for pure states, there are characteristic 
differences between the two different possible configurations. Furthermore, quite remarkably, in 
contrast to the bipartite setting, there is also pure true HE. When dealing with two parties, true 
HE only occurs for specific types of mixed states (compare section 4.4). Actually, this can be 
generalized. True hybrid entanglement can be obtained in two cases: Either the system is mixed 
with an infinite number of mix terms, then it is sufficient that only one subsystem is CV, or two 
or more subsystems are CV, then the state is not even required to be mixed. However, every 
mixed state with a finite number of mix terms, and which contains only one CV subsystem is 
effectively DV. A mixed state of this type in % dl <g> . . . <g> Wd N _ x ® can be always described 
in a Hilbert space of the form <g> . . . (g> %d N _ x ® "He, where 

N-l 

E = M n di, (4.147) 

1=1 

and M is the number of mix terms. In the following subsections two explicit tripartite HE states 
are investigated, one living in a Hilbert space of the first form % dy <g> H d2 ® and the other 
one in a Hilbert space of the form Hd 1 <8> "Hco <8 Hoo, again with finite d\, di. 

4.5.3. A Pure ^20^2® ^oo State 

Consider the state 

M ABC = ^ (V |0) B |^> C + |1) B 1^1)°), (4.148) 

living in a Hilbert space of the form ® 7{| ® "H£. It contains only two qumode states \ipo) C 
and |i/>i) C . Hence, perform a Gram-Schmidt process: 

|^o> C = |0) C , (4.149) 
|^i> C = Q|0) c + yi-|Q| 2 |l) c , (4.150) 

where the overlap between the qumode states is denoted by Q := (tf \ipi). This yields 

\tp) ABC = (\0} A |0) B |0) c + Q \1} A |1) B |0) c + yfl- \Q\ 2 \1) A |1) B |1) C ), (4.151) 

which is a three-qubit state in V.^ ® %f ® ^2 • Hence, it can be easily analyzed with aid of 
the concurrence. This is the reason why such a HE state which is effectively a qubit triple has 
actually been chosen. Only in this case the state can be examined regarding its bipartite and 
residual, genuine tripartite entanglement. For higher dimensional states this would not be 
possible, as the residual entanglement is defined via concurrences, which only work for qubits. 
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Tracing out one subsystem and calculating the squared concurrences of the remaining two- 
qubit states yields the following results: 



C 2 (A\B) = \Q\ 2 = \(ip \ipi)\ 2 , (4.152) 
C 2 (A|C)=0, (4.153) 
C 2 (B|C)=0. (4.154) 



Furthermore, 
and hence, 



C 2 (C\AB) = 1 - | Q| 2 = 1 - | <Vo|^i> I 2 / (4.155) 



T res = l-|Qr = l-K</'ol</>i>r. (4.156) 



It can be seen that the state offers maximal residual entanglement for Q — 0. Since Q = 
corresponds to orthogonal |i/>o) C an d IVi) ' m m i s case tne state is simply the GHZ state 
itself. However, for the other extreme Q — 1, and therefore |i/>o) C = l^i) ' me state i s onr y a 2- 
entangled product state of the form |i/>o) C / which shows maximal bipartite entanglement 

between systems A and B. Therefore, by tuning the overlap Q between the qumode states, it can 
be gradually switched between genuine tripartite entanglement and bipartite 2-entanglement 
between systems A and B. This corresponds to varying between scenario A) and scenario C) in 
figure 4.15. Actually, the total entanglement in the state stays always constant and assumes the 
value of one ebit: 

C 2 {A\B) + C 2 (A\C) + C 2 {B\C) + T res = 1. (4.157) 

The state is therefore a 1-ebit entangled state for any Q. The overlap Q only determines the 
distribution between normal bipartite 2-entanglement and genuine tripartite entanglement. 8 

Tuning the overlap between two such qumode states should not be very difficult experimen- 
tally. The qumode states can be simply realized by two opposed coherent states |±a) , which 
results in Q = e~ 2 l a l 2 . Then the overlap can be tuned by simply adjusting the amplitude a. In 
this way, relatively small overlaps can be achieved already for low amplitudes. The problem is 
rather to establish the entanglement itself between the qumode and the qubits. 



4.5.4. A Pure H 2 <8> ^oo ® ^oo State 

Now the related state, 

m ABC = ^ (\0) A \cPof |^o> C + II)* |<h) B l^i) C ), (4-158) 



8 The W-state actually yields C%,(A\B) + C^(A|C) + C^(B\C) + r„ Sw = \, which is greater than one. However, when 
only permitting GHZ-like 3-entanglement and "normal" 2-entanglement but no W-like 3-entanglement, in general 
C 2 (A\B) + C 2 (A|C) + C 2 (B\C) + r res < I- Hence, in this sense the state \ f) ABC can be regarded as maximally 
entangled. 
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ipir = S ip |0> L + v /l-|Q lf | 2 |l} L , (4.162) 



which contains two CV subsystems and is supported by the Hilbert space ® <8> 
shall be discussed. Defining Q^, := ((po\(pi) and := {ipo\ipi}, a Gram-Schmidt process yields 

l4>o) C = |0> B , (4.159) 

|^1> C = |0) B + - |Q^| 2 |l) B , (4.160) 

|</>o) C HO) C , (4.161) 
ri \ c - n ... m\ c j- . R~L in ,12 n\C 

and therefore, 

\$) ABC = -j= (\0) A |0> B |0) c + QtpQy \1) A \0) B \0) c + V /(1-|Q^| 2 )(1-|Q^| 2 )|1> A |1) B |1) C 
+ Q^l-\Q f \ 2 )\l) A |1> B |0) c + Q^l-\Q^\ 2 )\l) A |0) B |l) c ). 

(4.163) 

Once again, the state is effectively a three-qubit state in <8> "%f ® ^2 ■ Therefore, calculate 
the relevant squared concurrences: 

C 2 (A\B) = 12/(1 - |S/) = | | 2 (1 - | (fclfc) | 2 ) (4.164) 

C 2 (A|C) = |Q/(1 - |S^| 2 ) = | <4> |<fr> | 2 (1 - I (#>|fc> I 2 ), (4.165) 

C 2 (B\C)=0, (4.166) 

C 2 (A|BC) = 1 - IS/IS/ = 1-| <</> |<fr> | 2 | (ifoltfi) W (4.167) 

Hence, 

Tres = (1 " |fi*| 2 )(l " |Q^| 2 ) = (1 " | (tf> |tf>l> | 2 )(1 " I <^o|^l> (4.168) 

which is plotted in figure 4.16. Furthermore, the total entanglement C 2 0tal is given by 

C 2 0tal = C 2 (A\B) + C 2 (A\C) + C 2 (B\C) + r res = 1 - | (ftlfr) | 2 | (Vol^i) | 2 - (4.169) 

These results provide a basis for interesting discussions. Figure 4.16 looks quite unspectacular. 
However, it contains actually interesting information in combination with the other results. 
First, it can be seen that for both overlaps ((po\<pi) and (ipo\tpi) being zero, the state is once 
again exactly the GHZ-state, which shows only GHZ-like 3-entanglement but no bipartite 
entanglement. In the case that only one of the overlaps is zero, the state considered in the 
previous subsection is obtained. It is always 1-ebit entangled, while it can be tuned between 
GHZ-like entanglement and common bipartite 2-entanglement. If the overlap which is not zero 
becomes one, maximal Bell state-like 2-entanglement occurs. Finally, if none of the overlaps 
is zero the state is not 1-ebit entangled anymore and it stands in a mixture of bipartite 2- 
entanglement and GHZ-like 3-entanglement. If both overlaps are one, the state becomes simply 
fully separable, showing no entanglement at all. All this is basically reflected in figure 4.16. 
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Figure 4.16.: The residual entanglement of the state \ip) as a function of the overlaps ((po\cpi) 
and (ipo\ipi), which were assumed to be real without loss of generality. In the 
upper corner, for T res = 1, the state is maximally GHZ-like entangled. It shows no 
bipartite entanglement. In the left and the right corner, the residual entanglement is 
zero. Nevertheless, the state is maximally entangled, but just bipartite 2-entangled. 
Finally, in the lower corner, the state is fully separable. For either (<po\(p\)=0 or 
(ipo|i/?i}=0, the state is always 1-ebit entangled, while for both (<po\(pi) 7^ and 
(ipol^i) 7^ the state has less than 1-ebit entanglement. 

The upper corner of the graph corresponds to a GHZ-state, while the left and the right corners 
denote states which are maximally bipartite 2-entangled. Every point in the graph which is not 
at the upper left or upper right edge corresponds to a less than 1-ebit entangled state, while 
states lying on these edges are 1-ebit entangled. States on the lower edges of the graph are 
only 2-entangled, while all other states show also GHZ-like 3-entanglement. Finally, the lower 
corner corresponds to a fully separable state. 

Coming to a conclusion, it is pointed out once again that it is quite remarkable that just 
by tuning simple overlaps between states, it is possible to gradually switch between all these 
different entanglement scenarios. Furthermore, tuning overlaps is a relatively easy task, 
when the qumode states are simply realized by ordinary coherent states. Then only the 
amplitudes have to be adjusted. However, the experimental preparation of the overall entangled 

~ ABC 

tripartite state | ip) may cause difficulties. Nevertheless, this idea of tuning between various 
entanglement configurations by modification of the overlaps of the participating qumode states 
may be a scheme which could possibly be experimentally realized in future. 



5. Hybrid Entanglement in Experiments 



In the previous chapter hybrid entanglement has been illuminated from a theoretical point 
of view. Different classes of HE were identified. Furthermore, approaches regarding the 
quantification as well as the detection of hybrid entanglement have been discussed in-depth 
on a theoretical level. Additionally some peculiarities of multipartite HE have been examined. 
This is all very nice, but naturally two questions arise. On the one hand, how can HE actually 
be experimentally generated? And on the other hand, what is it useful for? 

These questions are addressed in this chapter. First, an experimental scheme is introduced, 
which achieves the generation of pure entangled qubit-qumode states. When sending these 
systems through a channel, mixed HE or even true HE may be obtained. Second, three 
experimental applications are briefly described. On the one side, there are schemes regarding 
quantum key distribution as well as the generation of Schrodinger-cat states, which both 
involve pure entangled qubit-qumode states. On the other side, qubus appoaches for quantum 
communication and computation actually utilize multipartite HE. 



5.1. Generating Hybrid Entanglement and Cat State 
Preparation 

Consider a single 2-level system (qubit) interacting with a single mode of the quantized 
electromagnetic field. Such a configuration can be easily achieved by, for example, placing 
an appropriate atom inside an optical cavity which supports a single photonic qumode. The 
required 2-level "atom" may for instance be realized by a quantum dot [BSG + 05]. Then, in 
the Jaynes-Cummings-Paul model (/CM), applying the rotating-wave approximation, the system's 
Hamiltonian is given by 

H JCM = fiD.a'a + hico&z + fig{aa' + a' a), (5.1) 

where CI is the electromagnetic field's frequency and to is the atom's transition frequency. 
Furthermore, a as well as a + are the fields mode operators, while a z is the qubit's Pauli Z 
operator acting on the atom's ground and excited states, |0) and |1) respectively. Additionally, 
a= |0) (1| as well as ^ = |1) (0| are the raising and lowering operators of the atom. Finally, 
the interaction strength is determined by the position dependent vacuum Rabi frequency 

g (f) = ^S , (5.2) 
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with the vacuum electric field inside the cavity S q, the atom's dipole moment p, and the mode 



function u(f) of the resonator [SchOl, JC63, Pau63]. Note that the atom is assumed to be at rest. 

P 

2m 



Otherwise, additionally a center-of-mass motion term would occur in the Hamiltonian. 



The JCM is one of the most fundamental models in quantum optics when considering atom- 
light interactions. It is in principle the most elementary model which is nevertheless able to 
describe most phenomena in cavity quantum electrodynamics. 

The most important term in the Hamiltonian is the interaction term fig{aa^ + <J + fl), which 
describes the dynamics in the JCM. In the interaction picture, 

fil™ = hg(M f e m + <7 + ae- ,At ) (5.3) 

is obtained, with the frequency detuning A :— O — to. Most interesting are the dynamics in 
the limits of this model: The resonant case A = and the far off-resonant case with very large 
A. As the far off -resonant case conserves photon numbers as well as atomic populations, and 
introduces only detuning dependent phase shifts, it is also called the dispersive limit. This is also 
the regime which is of relevance for the preparation of HE states. Hence, we shall no longer 
consider the resonant case. 

More precisely, the dispersive limit is reached if 



with photon number n. In this case the effective Hamiltonian, 

H f// = - h ^(& z a i a+ l -(a z +\)), (5.5) 

can be used to describe the system [SchOl]. Once again, the interesting term, necessary for the 
generation of HE, is only the interaction term. Hence, just consider 

2 

fi$ = -hxte% X-=^- (5.6) 
Then the corresponding unitary evolution operator is given by 

U = expl-^H'K t] = exp[f^ z a + fl f], (5.7) 

which represents a simple rotation in the qumode's phase space [vLll]. Hence, define cp :— xt 
and rewrite this equation as 

U = R(cpcr z ) — ex-p[icpa z a f a]. (5.8) 

R{cp& z ) corresponds to a conditional rotation, acting on the whole Hilbert space of the qubit 
and the qumode. Applying it on a coherent state | a.) yields 
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Figure 5.1.: A coherent state |a) with a 6 R is conditionally rotated by ±<p in phase space, 
depending on the state of the qubit. This procedure results in the qumode states 

|ae ±! >). 



R{cp£r z ) | a) 



exp[icpa z n) e 2 V~! ~j= \ n ) 



_]«r 

e 2 E 

n=0 



(5.9) 



a e' 



The state is rotated in phase space by a value of icp depending on the state of the qubit (see 
figure 5.1). At this point it becomes obvious that in the dispersive limit only phase shifts are 
acquired, but no photon numbers or atomic populations are touched. The most interesting case 
occurs, when the qubit is in a superposition state ^j(|0) + 



H<p**) l*> ® ^(lo> + |1>) = l«e'^> ® -^(|0> 

1 



1 



- ^(ke l >)|0) + | a e-^)|l)). 

Finally, this is a hybrid entangled state between a qubit and a qumode of the form 

\<P)he = c \0) a \tP ) b + Cl \l) A \^) B , 



(5.10) 



(5.11) 



as considered in paragraph 4.2.1. Concluding, it can be summarized that the JCM provides 
an efficient way for the generation of HE between a qumode and a qubit. Of course the 
qumode may also consecutively interact with several qubits, or similarly, several modes may 
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possibly interact with a single or several qubits. This would result in a multipartite HE state. 
Since cp = xt with typically very short interaction times, only small phases can be acquired. 
Nevertheless, for sufficiently large amplitude a, the overlap between the states {|ae ± "P}} can 
become arbitrarily small. 

Finally, from the state of equation (5.10) CSS states can be created by a suitable projection 
measurement on the qubit [SBW90, vLll]. Simply perform a projection on the conjugate 
(j^-basis |±) = -y= ( 1 0) ± |1)) and obtain cat states of the form 



Hence, besides the hybrid cat state generation technique with aid of CV squeezed vacuum 
and DV photon subtraction [LJG+06, OJTBG07], which has been presented in section 3.2, this 
approach constitutes an alternative way that exploits hybrid entanglement. 

Summarizing, in this section an efficient scheme for the generation of qubit-qumode HE has 
been introduced, which also sets out the basis for general multipartite HE creation. Furthermore, 
a first application of HE has been presented - the CSS state generation. 



Another interesting application of hybrid entangled states is found in quantum key distribution 
(QKD) [GRTZ02]. In QKD two parties, Alice and Bob, aim at distributing a shared random 
key, which only they know, and which can be used for encrypting and decrypting messages. 
They make use of the fundamental laws of quantum mechanics to guarantee security. An 
eavesdropper, Eve, must somehow measure the quantum states, which Alice and Bob exchange, 
to gain information. However, in quantum mechanics measurements disturb the system. 
Therefore, Alice and Bob are able to detect Eve. Furthermore, the amount of intercepted 
information can be calculated and bounded. If the shared information between Alice and Bob is 
greater than the information Eve could get, nevertheless a secret key can be distilled via privacy 
amplification, effectively eliminating Eve's amount of information. In privacy amplification the 
shared key, which Eve has partial knowledge about, is used to distill a shorter key, which Eve 
has less knowledge about. This can be repeated continuously to, in principle, completely zero 
out Eve's information about the final key. Finally, the shared secret key can be used for the 
execution of a one-time pad ciphering to interchange unconditionally secure messages. 1 

In principle, there are two different approaches towards QKD. On the one hand, there are 
so-called prepare and measure schemes, where Alice prepares an entangled state first, measures 
her part, and finally sends the remaining part through an insecure domain to Bob. On the 
other hand, there are entanglement-based approaches in which an entangled state is provided 
by an untrusted third party. Afterwards Alice and Bob perform measurements on this state 
and basically quantum teleport the information from Alice to Bob. An example for a prepare 

1 For a detailed introduction to QKD see the review [GRTZ02]. 




(5.12) 



5.2. Quantum Key Distribution 



5.2. Quantum Key Distribution 



97 



and measure approach is the well-known BB84 protocol by Bennett and Brassard, which is 
actually the first suggested QKD protocol at all [BB84]. In 1991, Ekert presented the first 
entanglement-based scheme, the E91 protocol [Eke91]. 

At the end the prepare and measure protocols as well as the entanglement-based schemes are 
more or less similar. The only difference is, where the entangled state is prepared. Hence, for 
a given prepare and measure protocol always the corresponding entanglement-based setup 
can be written down, and vice versa. It is the same key point in both configurations which 
determines whether unconditionally secure QKD can be successfully performed or not: The 
correlated data between Alice and Bob have to violate some Bell inequality [AGM06]. Hence, in 
other words, the data have to correspond to an entangled state. Therefore, security in QKD 
protocols is in principle an inseparability problem. Actually, there are schemes which involve 
hybrid entangled states. Such a protocol shall be presented now [LKL04, RGL06]. 

Suppose Alice prepares randomly entangled states of the form 



| e i(n+2)f 



(5.13) 



with two bases n = V 1 (compare to figure 4.8). Then she measures her part of the state 
and sends the remaining state to Bob through an insecure domain controlled by Eve. Bob has 
to make measurements such that he can distinguish between |e" ! 2 a) and |e^" +2 -' 2 a) (see 
figure 5.2). So he chooses either the n = or n = 1 basis, performs appropriate measurements, 




l+l 0£> 



Quantum channel, 
controlled by Eve 



Heterodyne Detection 



Bob 



Classical, public, and authenticated channel 



Figure 5.2.: The explained QKD scheme: Alice prepares entangled states \ip) n , and measures 
her part of the state. The remaining coherent states are sent through an unsecure 
domain to Bob, who performs for example heterodyne measurements on the state 
[RGL06]. Also an unambiguous state discrimination setup as introduced in section 
3.2 would be reasonable. 



and distinguishes between |e m ?a) and |e ! '" +2 )5a) for the chosen n. Since a is chosen such 
that the states have a significant overlap, the distinguishment between the states can only be 
performed with a certain error probability. This is a very important feature in this scheme 
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and actually in most QKD protocols. After the transmission Alice and Bob talk to each other 
via a classical, public, and authenticated channel, and tell each other, which basis they have 
chosen. 2 They keep only the bits which have been measured for an equal choice. Furthermore, 
all bits corresponding to uncertain measurement results are discarded and only bits are kept 
which Bob knows with high certainty. This process of cheosing only favorable events is called 
postselection and is one of the key ingredients necessary for successful QKD [SRLL02]. 

Consider an eavesdropper Eve. She has mainly two possibilities for attacks. The intercept and 
resend attack as well as the beam-splitting attack. 

An intercept and resend attack works in the way that she simply measures the state, prepares 
a new state corresponding to her measurement result and resends it. However, Eve also has to 
choose a basis for the measurement to obtain a reasonable result. Then there are two cases: On 
the one hand, Eve may choose the right basis, probably learns the bit's value and resends the 
right state (due to the overlap, she may also obtain a wrong measurement result and resend 
a wrong state. However, this may be noticed by Alice and Bob, as will become clear soon). 
Alice and Bob can not notice this. On the other hand, she may choose the wrong basis. Then 
she resends the state also in the other basis and Bob may measure the wrong bit, even though 
he chooses the initially right basis corresponding to Alice's basis choice. In this case Alice 
and Bob notice Eve: When they compare the bases they have been applied and find that they 
have used the same bases but measured too often different bits, they know that they have 
been eavesdropped. These measurement results which they compare are the control bits. The 
statistics of these measurement results will reveal any eavesdropper. Of course, these control 
bits cannot be used for the key anymore, as the classical authenticated channel is public. If Alice 
and Bob perform a sufficiently large number of control measurements and observe anomalous 
statistics in the measurements, they can conclude that they have been eavesdropped and can 
start the key transmission again. Hence, this protocol is secure against intercept and resend 
attacks. Actually, it is not even necessary to begin with a new key transmission as due to 
the postselection Alice and Bob share a greater amount of information than Alice and Eve 
do, anyway This enables key distillation. It will become clear in the next paragraph when 
discussing beam-splitting attacks. 

So, furthermore, to gain information about the key Eve can tap parts of the signal by 
a beam splitter with transmission rj (beam-splitting attack). Then Eve gets an amount of 
1 — t] of the signal and therefore states of the form — j/e m ?a) and | v/l — ne l ( n+2 -^ a) . 
However, this also yields attenuated states for Bob. When performing a large amount of 
measurements Bob notices that the attenuation is higher than expected, and hence can conclude 
that the transmission has been eavesdropped. Nevertheless, especially for very large rj he may 
possibly not notice. Even so, Eve will not be able to take possession of the whole key. For 
her measurements she also randomly chooses the basis n, independent of Bob's later choice. 
It is a crucial point that her measurements are fully independent of Bob's measurements. 

2 It is essential to employ an authenticated channel for the classical communication after the measurements. If this is 
not the case, the protocol will be vulnerable to man-in-the-middle attacks, where Eve impersonates Bob, when talking 
to Alice, and vice versa. 
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Furthermore, she also listens to the classical conversation between Alice and Bob, in which 
they compare the chosen bases. With the control measurements Alice and Bob cannot find 
Eve. She remains hidden. Therefore, the key corresponding to the chosen bits is actually used. 
Due to the postselection, Bob has only kept bits which correspond to measurements with very 
certain results and the right choice of basis. However, it is highly improbable, actually nearly 
impossible for a large number of events, that in case of a "good" measurement of Bob, also Eve 
always performs a "good" measurement, for which she knows the bit with high certainty. The 
postselection of favorable "good" measurements is done by Bob and not by Eve. The postselected 
events do not, in general, all correspond to "good" measurements of Eve. Therefore, Eve has 
more erroneous relevant measurements than Bob. Hence, the amount of shared information 
between Alice and Bob, I^g, is larger than the amount of information held between Alice and 
Eve, Iae- This translates into only partial knowledge of the key. Finally, Alice and Bob can 
perform privacy amplification to eliminate Eve's possible knowledge about the distilled key 
completely. Concluding, the scheme is also secure against beam-splitting attacks. Note that, 
as mentioned, the possibility of key distillation due to postselection also occurs in the earlier 
discussed scenario for intercept and resend attacks. Even though Eve is actually detected in 
this case, a key can nevertheless be distilled. 

In reality, things are slightly more difficult, as the channel from Alice to Bob is not only under 
the control of Eve but also possibly introduces further losses and noise. These effects may 
distort the states in such a way that no data can be obtained which is suited for key distillation. 
However, as noted earlier, QKD can be performed if the data shows nonclassical quantum 
correlations. Hence, denoting the channel by Y B , the security of the key depends on whether 
p' n AB = (t A ® Y B ) \ip) AB (tp\ is entangled or not. Actually, p' n AB is a hybrid entangled state, 
and depending on the channel it may also show true hybrid entanglement. The detection of 
entanglement in states of this form has been discussed in detail in subsection 4.4.1, where 
the channel was considered to be lossy and noisy. Note that there are further approaches for 
witnessing entanglement in this class of states, which are based on special EVMs related to the 
original SV determinants [RGL06, HML08]. 

Concluding, hybrid entanglement and especially its detection is highly relevant for certain 
QKD protocols. 

5.3. Quantum Bus Schemes 

As explained earlier, entanglement is a key ingredient for most applications in quantum 
information. However, how is entanglement between, for example, two qubits experimentally 
established? Nonlinear interactions between the qubits to be entangled are necessary for this 
task. On the one hand, there are measurement-induced nonlinearities, which are for example 
heavily exploited in the quantum computation scheme by KLM [KLM01]. These approaches 
are highly non-deterministic, though. On the other hand, weak nonlinearities are actual 
physical nonlinear effects. They are typically much too weak for the realization of appropriate 
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entangling gates when dealing, for example, with single photons. An elegant way out of this 
dilemma is provided by the already mentioned quantum bus (qubus) approach, which utilizes 
an intense qumode to mediate and enhance the interaction between two qubits to be entangled 
[SNB+06, CZ95, vLMN+08, vLLMN08]. In principle, the qumode interacts subsequently with 
both qubits, resulting in a tripartite qubit-qubit-qumode HE state (see figure 5.3). Afterwards, 
an appropriate measurement on the qumode is performed forcing the qubits to remain in an 
entangled state. 



r , . - t i Qubus Qubus . , . „ 

Qubit 1 v Q u bu 2 




Figure 5.3.: The quantum bus entangles itself with the first qubit, then travels to the second 
qubit and also interacts with it. Afterwards an appropriate measurement on the 
qubus is performed to establish entanglement between the qubits [vLLMN08]. This 
figure is a modified version of Fig. la in [vLLMN08]. The figure was used with 
permission by Peter van Loock. 



This shall be explained in slightly more detail now, according to the references [vLLMN08, 
vLll]. Start with a coherent state |«) as qubus and a first qubit which is in the state 4|(|0) + 11)). 
If they interact dispersively as explained in subsection 5.1, the resulting hybrid entangled state 
has the form 

1 / . . \ . , , 1 / J6 s 



|0) + |l)j®|*) -> ^(J0}k) + |l}|«0| (5-14) 

up to an irrelevant phase space rotation. The qubus subsystem is now sent through a channel 
to the second qubit. However, this channel may introduce losses. These are again modeled by 
coupling to a vacuum environment as illustrated in figure 4.4. Hence, after the channel 



|0> |«) + |1) \*e' e )^ p = \ (JO) (0| ® ly/TfCt) (y/tjd 
+ |1) (1|<8> \^jae w )(^ae w \ 



(5.15) 

+ A\0) {\\®\^){^ae l6 \ 



+ A*\1) (0\®\^ae w )(^a 

where t] denotes the transmissivity of the beam splitter, and A := {yf\ — r]txe ie \^l — not). Now 
the qumode interacts again dispersively with the second qubit, resulting in another controlled 
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phase rotation of the qumode. This interaction then yields a tripartite HE state of the form 



p -> F |Y+) (Y+| + (1 — F) |Y~) (Y~ | , 
1 
2 



|Y±> := 71 1 ^ a> l<1>2±> 111 '^ ae,W) |10> + ^ v ^ oce ~ l0) |01) ' 



m) = 



i 

71 



(|00)±|11», 



cp := rjoi smf 



:= I[i + e -(i-'7)« 2 (i-cose)]^ 



(5.16) 
(5.17) 

(5.18) 
(5.19) 
(5.20) 



Finally, during the last step of the scheme the qubus disentangles from the qubits, projecting 
them on two-qubit entangled states. This proceeds via an appropriate measurement on the 
qubus mode. Homodyne detection along the x-axis is experimentally most efficient, discriminat- 
ing between \-Jfjcc) and \ y/rjae ) [vLLS + 06]. This projects the qubits on approximate versions 
of \<f>2) or TTld-^) However, as the distance on the x-axis between the peaks of the 

coherent states scales as ~ u.8 2 for small 8, the amplitude a has to be very large to compensate 
for the small 8 and to suppress overlap errors. Unfortunately, this increases the decoherence 
induced by imperfect channel transmission, too. A schematic of the whole experimental setup 
is shown in figure 5.4. Note that there are also other ways for measurements. For example 
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Figure 5.4.: Possible experimental realization of the quantum bus scheme, which distributes 
entanglement between two separated qubits. It is important to note that for the 
homodyne detection a strong local oscillator (LO) pulse has to be taken along 
[vLLMN08, vLLS+06]. This figure is taken from [vLLMN08] with permission by 
Peter van Loock. 



homodyne detection along the p-axis on the one hand results in lower overlap errors, yielding 
higher fidelities than in the x-axis homodyning scheme. On the other hand not all outcomes 
of the measurement correspond to a two-qubit state which remains entangled. Hence, several 
outcomes must be discarded, rendering the scheme less efficient. Another approach consists 
in utilizing POVMs and performing unambiguous state discrimination instead of homodyne 
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detection [vLLMN08]. In the end, for every measurement scheme there is some trade-off 
between fidelity and success probability. 

Finally it is worth mentioning that for a qubus implementation inside, for example, a hybrid 
quantum repeater setup, also hybrid entanglement distillation as well as hybrid entanglement 
swapping has to be performed [vLLMN08]. These shall not be discussed here, however, for 
more information, see the reference. 

It can be concluded that not just bipartite hybrid entanglement is experimentally relevant, but 
also multipartite HE. Quantum bus schemes are quite popular, since they are very important 
for the experimental realization of several tasks, especially regarding quantum communication. 
Hence, it is important to also investigate the involved multipartite hybrid entanglement, which 
has been done briefly in this thesis in the previous chapter. 



6. Summary and Conclusion 



The central theme of this thesis has been the investigation of hybrid entanglement with a focus 
on bipartite systems. So far, no rigorous and systematic analysis of hybrid entanglement has 
been performed. This has been the goal of this thesis, and actually, a lot of new insights were 
obtained. However, also several other related questions, which are to some extent associated 
with the subtleties of hybrid entanglement and have therefore arisen during the studies, have 
been discussed and in most cases successfully answered. 

The thesis starts with a brief introduction, chapter 2, explaining some fundamentals of 
quantum information theory. However, not really basic anymore, in subsection 2.4.5 general 
entanglement witnessing with aid of the SV criteria in two-mode cat states is discussed, and 
a generalized determinant is derived which achieves entanglement detection for states of the 
form %=p (|«,fli) + e 1 ^ |— a, — a)) for any phase (p. This is quite interesting, but, nonetheless, 
the witnessing is only demonstrated for pure states and would be much more relevant in the 
mixed state case. Further investigations should concern the performance of the generalized 
determinant for the corresponding mixed states obtained when sending the considered two- 
mode cats through decoherence-causing channels. Additionally, also even more general two- 
mode CSS states, such as the ones in equations (2.127) and (2.128), may be analyzed with the SV 
criteria. Finally, general issues on entanglement witnessing have been discussed, which bring 
out some very interesting questions to be tackled in the future: Are witnesses which work well 
for pure states also the right ones for witnessing in the corresponding mixed states, which are 
obtained by sending the pure states through decoherence-causing channels? Does there always 
exist an optimal witness in the sense that it makes use of all entanglement of a given state? If 
yes, it could be expected that this witness also detects entanglement in the decohered states if 
still present. Concluding, besides the two-mode cat state witnessing, also several interesting 
questions regarding entanglement witnessing have been briefly presented and discussed. It is 
exciting to await future results concerning these open problems. 

In chapter 3, the notion of hybrid as used in this thesis is defined. Furthermore, some 
motivations, why it is reasonable to "go hybrid", are discussed. Concerning this matter, the 
most important point probably comes from quantum computation and the question, how to 
realize universality. Hence, a very important protocol for hybrid quantum computation, the 
GKP scheme, is briefly presented and discussed [GKP01]. Especially important in this regard 
are the different notions of nonlinearity. On the one hand, there are measurement-induced 
nonlinearities, which are heavily utilized by KLM in their teleportation based scheme for 
DV quantum computation, while, on the other hand, weak nonlinearities are actual physical 
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nonlinearities, which are for example exploited and enhanced in the subsequently discussed 
qubus schemes (see section 5.3) [KLM01, SNB+06, CZ95, vLMN+08, vLLMN08]. Additionally, 
chapter 3 summarizes some applications involving hybrid approaches and provides a lot of 
reference material. Concluding, this chapter is mainly a motivation for the analysis of the 
central theme of the thesis, the hybrid entanglement, and serves as an overview over suggested 
applications which make use of hybrid approaches. 

The main topic of the thesis is discussed in chapter 4. It starts with a proof that every 
correlated hybrid state is non-Gaussian. This fact is not a priori clear and therefore not trivial. 
Then, the Gram-Schmidt process is discussed and a modified version, which has been called 
inverse Gram-Schmidt process, is derived. The latter is actually the central tool, when dealing 
with hybrid entanglement, and gives rise to a new classification scheme of hybrid entangled 
states: On the one side, states which contain a finite number of qumode states can, with aid of 
the inverse Gram-Schmidt process, always be expressed in finite orthonormal bases. Therefore, 
they are effectively DV hybrid entangled states. On the other side, for states possessing an 
infinite number of qumode states, the inverse Gram-Schmidt process is not applicable anymore. 
Hence, these states "stay" in their hybrid Hilbert space and are called truly hybrid entangled. 
This classification constitutes the central result of this thesis and it is quite relevant: Since 
hybrid entangled states are initially always non-Gaussian and live in an infinite-dimensional 
Hilbert space, exact entanglement quantification is impossible. However, the effectively DV 
hybrid entangled states can be described by proper density matrices after the application of the 
inverse Gram-Schmidt process. Hence, all the tools, such as several entanglement measures or, 
in the pure state case, the Schmidt decomposition, can be applied. It is quite remarkable that 
many of these complicated non-Gaussian and infinite-dimensional Hilbert space states can be 
handled so easily in the end. It is worth pointing out that the Gram-Schmidt process can be 
also applied on certain fully CV systems for which all subsystems are CV. The only requirement 
is that the superposition sum contains only a finite number of qumodes in each subsystem. 
Then, also these states can be handled with DV tools. Examples for this kind of CV states are, 
for example, the pure two-mode cat states, which have been considered in this thesis. 

In subsections 4.2-4.4, a large variety of exemplary states is presented and analyzed. For 
several states which are either too complicated for DV measures, or which are truly hybrid en- 
tangled and therefore stay in the non-Gaussian infinite-dimensional regime, exact entanglement 
quantification is not possible. However, at least entanglement detection can be achieved with aid 
of the SV criteria, which therefore are one of the most important tools for studying hybrid entan- 
glement. Most important during in investigations are states of the form ^j(|0) \oc) + |1) |— a)), 
as they are relevant in certain QKD schemes (compare section 5.2) [RGL06]. Entanglement 
quantification of states like this, when sent through a lossy channel, has been demonstrated via 
the concurrence. Furthermore, true hybrid entanglement is obtained, when taking the channel 
not only as lossy but also as noisy. Using this example it was shown that the SV criteria also 
work for the new true hybrid entanglement. Also another example of a true hybrid entangled 
state has been considered, (4.107), for which entanglement witnessing was accomplished. 
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When dealing with true HE, infinite sums or integrals always occur. It is important not to 
truncate these and by this obtain approximate results. A truncation would actually correspond 
to the analysis of a different state, strictly speaking not even a truly HE state anymore, but an 
effectively DV HE state. It appeared that this complicates the analysis of truly hybrid entangled 
states considerably. This has led to the examination of general issues regarding entanglement 
witnessing in mixed states with an infinite number of mix terms. It has been shown that, 
when the global robustness of the truncated state (truncation in the convex combination of 
mix terms, i.e. pure state projectors) is higher than a certain amount 1 — p N depending on 
where truncation has been performed, the overall state with an infinite number of mix terms is 
entangled. Additionally, it was derived that, for every entangled state of this form, a truncated 
state exists such that its global robustness is greater than 1 — p^. Due to the complicated 
calculation of robustness monotones these theorems are of limited practical use. Nevertheless, 
they provide insight into the connections between entanglement witnessing and entanglement 
quantification. Furthermore, one may actually think of special scenarios, where no witnessing 
can be performed for the given states, while for any truncated state, robustness calculations can 
be performed. Hence, these theorems can still be considered as a nice result, especially for the 
case that, in the future, efficient algorithms for the computation of the robustness are found. 

In addition, not only bipartite HE has been investigated, but also multipartite HE (section 4.5). 
While in the bipartite case true hybrid entanglement only occurs for certain kinds of mixed states 
(the ones with an infinite number of mixed terms) it has been shown that in the multipartite 
case true HE may be also obtained for pure states. Any hybrid system which contains at least 
two CV subsystems can be truly HE. When the overall system contains only one CV subsystem, 
true HE only appears for the mentioned mixed states. Additionally, investigations regarding 
the entanglement structure of tripartite HE states have been performed. It has been worked out 
that, depending on the overlaps between the involved qumode states, either no entanglement at 
all, bipartite entanglement between two subsystems, GHZ-like genuine tripartite entanglement, 
or a mix of the latter is obtained. This is a very interesting finding, as it is quite easy to 
tune overlaps: Just use coherent states and adjust their amplitudes. Possibly, this continuous 
switching between different entanglement scenarios via overlap tuning is something which can 
be achieved experimentally in the future. 

Finally, to round off the thesis and show that hybrid entanglement can be actually experi- 
mentally realized and is also physically relevant, chapter 5 discusses the generation of hybrid 
entanglement via the dispersive interaction between a qubit and a qumode in the Jaynes- 
Cummings-Paul model. Additionally, some applications, which employ hybrid entangled states, 
are briefly explained. First, via projection measurements on the DV part of qubit-qumode 
entangled states, Schrodinger cat states can be generated. Second, a certain QKD protocol is pre- 
sented, which exploits hybrid entanglement between qubits and binary coherent states {|±a)}. 
The necessary requirement for unconditionally secure key distribution then corresponds to 
inseparability of the mentioned entangled states when subject to the quantum channel, which 
connects Alice with Bob. Fittingly, in this thesis entanglement witnessing in states of this kind 



106 



6. Summary and Conclusion 



is discussed quite in detail (for the thermal channel). These investigations are therefore of 
very high relevance. Third, the quantum bus principle is introduced. A qubus mediates a 
nonlinear interaction between two separated qubits and can therefore serve as an entangling 
gate, important for several applications in quantum computation and quantum communication. 
The approach makes use of tripartite hybrid entangled states, showing that also studying 
multipartite entanglement is very important. 

Altogether, it can be concluded that the thesis provides a lot of new insights into hybrid 
entanglement and associated topics. The motivation of this work has come from questions 
concerning the classification of hybrid entanglement, its detection and quantification as well 
as its generation and its applications. A classification scheme has been derived and the 
detection of hybrid entanglement via the SV criteria has been demonstrated. Exact entanglement 
quantification can be either performed like in the DV case, or it cannot be accomplished at 
all, since the state is non-Gaussian and infinite-dimensional. A generation scheme has been 
presented and three applications have been discussed. 



A. Notations 



A.l. Abbreviations 

Throughout this thesis several abbreviations are used for the sake of appropriate readability. 
Here is an overview. 



Abbreviation 


Meaning 


LKW 


Coffman, Kundu and Wootters 


U 


Completely Positive 


CPTP 


C^nmnlptplv Positive Tr^rp-Prpservinp" 


CSS 


Coherent State Superposition, Schrodinger-Cat State 


CV 


Continuous Variable(s) 


DV 


Discrete Variable(s) 


EVM 


Expectation Value Matrix 


GHZ 


Greenberger-Horne-Zeilinger 


GKP 


Gottesman, Kitaev and Preskill 


HE 


Hybrid Entanglement / Hybrid Entangled 


JCM 


Jaynes-Cummings-Paul Model 


KLM 


Knill, Laflamme and Milburn 


LOCC 


Local Operations and Classical Communication 


NPT 


Negative Partial Transpose 


PnCP 


Positive but not Completely Positive 


POVM 


Positive Operator-Valued Measure 


PPT 


Positive Partial Transpose 


PT 


Partial(ly) Transpose(d) 


QIT 


Quantum Information Theory 


QKD 


Quantum Key Distribution 


SV 


Shchukin-Vogel 


TMSS 


Two-Mode Squeezed State 
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A.2. Mathematical Notation 

In addition to the summary of abbreviations here is also an overview over the mathematical 
notations used throughout the thesis. 



Notation 


Meaning 


R, C 


Field of real numbers, field of complex numbers 


R" 


Set of real n-dimensional vectors 


N, N , Z 


Natural numbers, natural numbers including zero, integers 


H(C) 


Hilbert space over the field of the complex numbers 




d-dimensional Hilbert space of system A, qudit space for finite d 


njA 


Infinite-dimensional Hilbert space of system A, qumode space 


l n , O n 


Identity operator in n dimensions, zero operator in n dimensions 


M* 


Elementwise complex conjugate of matrix M: (M*),y = M*- 


M T 


Transpose of matrix M: (M T )y = Mji 


M + 


Hermitian adjoint of matrix M, complex conjugation 




and transposition: M + = (M T )* 


A t b 


Partial transposition of operator A with respect to system B 


A® B 


Tensor product of operators A and B 


A(BB 


Direct sum of operators A and B 




Pauli matrices 


tr[A] 


Complete trace of operator A: tr[A] = (n\A\n) 


ti B [A],tr U B[A] 


Partial trace of operator A with respect to Hilbert space of 




system B: tr B [A] = (n B \A\n B ) 


det[M] 


Determinant of matrix M 


dim[V] 


Dimension of vector space V 


Sp(2n,M) 


Real, symplectic group on K 2 " 




Symmetrization of operators x n , p m 


SEP, QS 


Set of separable quantum states p se p, set of all quantum states p 


: A : 


Normally ordered form of operator A 


W 


Floor of the real number x: [x\ — max{m £ Z : m < x} 


Yp 


CPTP map Y(jo) on state p 


(l A ®Y B )p AB 


One-sided channel: CPTP map Y B on subsystem B of the state p AB 


x\Jy 


Logical disjunction, "or"-conjunction: x or y 


}a,b[ 


Open interval of real numbers between a and b: ]a,b[:= {x <E K : a < x < b} 


{«/} -< {M 


Set {oLi} is majorized by set {/3,} 



B. Calculations and Formulas 



This part of the appendix provides auxiliary material. For the understanding of the thesis often 
only theorems themselves are important but not their proofs. These proofs are shifted here. 
Furthermore, ancillary calculations are presented and some important definitions and formulas 
are listed. 



B.l. The Hadamard Lemma 



Lemma B.l (Hadamard's Lemma). For bounded operators X and Y, 



e 



[X,Y] 



m=0 

with [X,Y] m = [X, [X,Y] m -i] and [X, Y] = Y. 
Proof. Consider 

F{t) = e^Ye"' 1 , F(0) = Y. (B.2) 

Then 

dF ^± = [X,F]. (B.3) 
This differential equation has the solution 



n=0 

which can be seen by evaluation of the derivative: 

dF(t) nt»-i % f 



F (t)= L^jMn, (B.4) 



n=0 n[ (B.5) 

oo lyi oo .n 

= XE^[X,Y]„-E^[XY]„X 

n=0 n=0 

= [X,F]. 
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Evaluated at t = 1, this yields 



e*Ye-* = f (B.6) 



which completes the proof. □ 



B.2. The Schmidt Decomposition 



Theorem B.l (Schmidt decomposition). Let \ip) AB be a normalized pure bipartite state in the product 
Hilbert space l-i A ® 'Hf n with dm\[H A ] = n and dim[%^] = m. Then there exist orthonormal bases 
{\ui) A } and {\iVi) B } such that 

r(tp) 

W) AB = E V^i \ui) A \™i) B , «« > V i, = 1, (B.7) 

i=\ i 

with r(ip) < min{n,m}. Such a decomposition is called Schmidt decomposition of \ip) AB - 



Proof. [Aud07] Expand \xp) AB in the orthonormal bases {\k) A : 1 < k < n} of U A and 
{\i} . 1 < i < m} of 7if n . Assume n < m without loss of generality. 

n,m 

\q) AB = £ «ki\k) A \i) B . (B.8) 

k,i=l 

Introduce the relative states 

m 

\w k ) B :=J^a ki \i) B , (B.9) 

i=i 

leading to 

\ip) AB =t\ k ) A \^) B - (B-10) 

The relative states {|i%} B } are in general neither orthogonal nor normalized. Orthogonality 
can be obtained in the following way: 

1. For {IA:}" 4 } choose the orthonormal eigenvectors {\u k ) A } of p A : 

p A =Ep k \u k ) A (u k \, p k >0Vk, f> k = l. (B.ll) 

k=\ k=l 

Let p k > for 1 < k < r and p k = for r + 1 < k < n. The eigenvalues p k can be degenerate. 
The {I"*:}" 4 } are determined only up to a phase for non-degenerate eigenvalues. 



B.3. The Vacuum Density Operator 
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2. The relative states {|%} B } are now associated with the {\u k ) A }. Consider subsystem p A : 

p A = tr B lW AB (ip\] 

n 



k,l 
n 

= J^(Wl\w k ) |«fc) A <«/| 

k,l 



(B.12) 



Comparison of equation (B.ll) with equation (B.12) leads to the desired orthogonality of the 
relative states: 

(wi\u>k) B = Pkhi- ( B - 13 ) 

For k > r + 1, the {|*)c} B } ar e null vectors. Combining equations (B.10) and (B.13) as well as 
identifying {p k } with {a^} completes the proof. □ 

B.3. The Vacuum Density Operator 

Lemma B.2. Consider a qumode system with mode operators a and a*. The vacuum density operator 
|0) (0| can then be written as 

|0> <0| = :e-***:, (B.14) 

where : • • • : denotes normal ordering. 
Proof. 1. It is known that 



(* + y)" = E u * y ■ (B - 15) 

k=0 



Withx = 1 andy = -1, 



E(-l) fc (f) = <W (B.16) 



fc=0 

2. The vacuum density operator |0) (0| is defined by its action on a Fock state \n), 



|0)(0|n) =6on\0). (B.17) 



Consider : e fl a : \n): 



: e -« + « : |n> = E i^-U + "> |n> 
ml 



m=0 



= E ( ' 1)W In) 

m! (»-«)! (B.18) 



|n) 



m=0 

= Son |0) , 

where equation (B.16) has been exploited. □ 
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B.4. The Heaviside Step Function 

In subsection 2.4.5, the so-called Heaviside step function 0i has been used. However, varying 
definitions can be found in literature. Therefore, here is the definition used in this thesis: 

x < 0, 

x = 0, (B.19) 
x > 0. 

Further properties of this function are of no relevance for this thesis. Just note that it is related 
to the Dirac delta function via differentiation. However, in contrast to the delta function it can 
be also used as an actual function and is not only distributively defined. Furthermore, smooth 
approximations of 1 are given by so-called logistic functions, for example 

0i (x) = lim (\ + — arctan(fcc)V (B.20) 

2 jt^oo\2 71 J 




B.5. Hermite Polynomials 



The Hermite polynomials H n (x) ([BSMM05, TT96]) are defined by 

H n (x) :-- 

They are polynomials of order n. The first five Hermite polynomials are 



-l)'V 2 ^-e~ x \ 
; dx" 



H (x 
Hi(x 
H 2 (x 
H 3 {x 
H^{x 



1, 
2x, 

Ax 2 - 2, 

8x 3 - 12, 

16x 4 - 48x 2 + 12. 



However, there is also an explicit expression, 

where the floor function has been used. It is defined as 



(B.21) 



(B.22) 
(B.23) 
(B.24) 
(B.25) 
(B.26) 



(B.27) 



[x\ = max{m £ Z : m < x}. 



(B.28) 



B.6. Kraus Operators of the Thermal Channel 
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Additionally, the polynomials can be calculated recursively with the relation 

H n+ i (x) = 2xH n (x) - 2nH„_! (x), (B.29) 

which yields together with Hq(x) and H\{x) an efficient way for the actual calculation of the 
polynomials. 

2 

Furthermore, Hermite's polynomials are orthogonal with respect to the weightfunction e~ x , 

oo 

I H n (x)H m (x)e- x2 dx = 2 n nly/7tS nm , (B.30) 

— OO 

and they form an orthogonal basis of the Hilbert space of functions f(x) for which 

oo 

J \f(x)\ 2 e- x2 dx < co. (B.31) 

—oo 

Beside their stong importance in quantum mechanics, where they are relevant for the 
construction of the the energy eigenstates of the quantum harmonic oscillator, the Hermite 
polynomials find applications in the finite element method for the numerical solution of partial 
differential equations. 

B.6. Kraus Operators of the Thermal Channel 

In paragraph 4.4.1, the thermal photon noise channel has been investigated (recall figure 4.10, 
which explains the modeling of the channel.). Furthermore, a set of Kraus operators, describing 
the channel, has been derived, which is presented now. Consider the thermal channel 

^Lmal(P A ) =tr E [Lf AE (p A ® pf hermal ) U AE '\, (B.32) 

with 

jjAE _ e 9(a\a A -a\a E ) ^ (B.33) 

and 

Plermal = t/t ("I . Pf = (1 + < "^ > ) ) n + l • (B.34) 
This can be rewritten in the form 



oo oo 



Y thermal(P A ) = Pn' ^ninp A ^ mn , ( B -35) 
n=0 m=0 

with 

K mn = (m\U AE \n) E . (B.36) 
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Defining K mn := yj p%-K mn yields 

CO 

Y?hermal(P A ) = E K mn p A Kl ln . (B.37) 
n,m=Q 

Now, merge the indices m and n into an index k — k — (m, n) via an appropriate index ordering 
procedure. For example consider k x = (m,n), k 2 = (m,n) and order such that 

, Ifci I < \k 2 \ or 

fc 1 <fc 2 ^S ( B - 38 ) 
\ki I — fel an d ^1 <' ^2/ 

with \ki\ = m + n, and <' k 2 denotes that the first nonzero difference in — m, n. — n is positive. 
This is by the way the ordering, which has also been applied in the derivation of the SV criteria 
in section 2.4.4. Making use of this ordering, 

00 

^Lrmal(p A ) = L^p A Kt (B.39) 

k=0 

is obtained, which corresponds to an operator-sum decomposition of the thermal channel. 

The Kraus operators K mn are derived in the Fock basis as well as in the coherent state basis. 
The transformation from K mn to the desired operators K mn = is then trivial. 

1. Fock basis. 

Consider (m| U AE \n) E \k) A : 

(m| U AE \n) E \k) A = — L (m| U AE {a E ' ) n {a A " f |0} E |0>^ 
Vklnl 

= -L (m| - ^/— n a A "y{^a A ' + y^V^f \0) E \0) A 

Vklnl 

■ ^f-i{^f-i^,\a Ef y\o) E \o) A 

= 4n t E (") ( k )vv k - i+ 'V—v n - l+i (-rr-> 

Vklnl !=0;=0 \ tj\jj 

\i+j) E yj {i + j)\(n +k-i-j)\\n + k-i- j) A 



(m 



minfn.m} / \ / t \ 

b£_|(. 1 - ,) ^ W' + "-Vir' 

• \Jml(n + k — m)l \n + k — m) A 

ml(n+k- m)l n a 
j^j K mfc „(?7) |n + *-m) , 

(B.40) 



B.6. Kraus Operators of the Thermal Channel 
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with 

min{n,m} 



mini n,m> / \ / ; \ 

**.0j):= E ( n ^(J_i)VV k - m+2j V^ n+m - 2 \-^- (B.4D 
= £ ^ m!(W +^ m)! K mfc „( ?y )|n + fc-^(fc|. (B.42) 



!=0 

Hence, 



k=max{0,m-n} 

2. Coherent state basis. 

Consider (m\ U AE \ n) E \u.) A : 

(m\ U AE \n) E \a) A = -L (m\ U AE {a £i ) n D A (cc) \0} E \0} A 



X — (m\ - jT^a Af ) n b E {Jl - W )b A {JJja) |0) E |0>^ 



1 E (D W(-v/i - n^) n -\ljz^ (m - i\V^v«f \y/n*Y 



'n\ £g w v y m - l r- 

minjn.mj- / \ / ■ 



»! 5 V' 



(m — i)\ 



■e-^-^\ 2{ ^^\^a) A . 
y/{m-i)\ 



(B.43) 



Therefore, 



pCoh. 



-, mmfn.m} / , 

Jc sj{m-i)\ 



(B.44) 



Shorter, 



min{n,m} 

££f = E ^ m /('/)(« A+ ) n -'' / d 2 ^ mi (^,a) («| ^ - a), (B.45) 

i=0 7C 



with 



and 



^■■-{oini^s-^^ (b - 46) 

fiUfrg) :^e-^H)tfUp)r, (B .47) 



The delta distribution 5{Q, which has been made use of, is defined by 



jfd 2 W -&>)/(£) =/(£<>), (B.48) 



116 



B. Calculations and Formulas 



with Z J , £o G C. Furthermore, make use of 

S{f[a}) \cl)=8{J[x\)\x), (B.49) 

for coherent states | a) . 

Finally, note that the Kraus operators ' Coh ' for n = are the Kraus operators of the 
amplitude damping channel, which can be also found in [NCOO]. 



B.7. Integral Identities 

In subsection 4.4.1, some important integrals have been used. These shall be calculated here. 

X 2 

• The integral J„ dx e~ ~* . 

f dxe-^ = J~a [ dye-y\ (B.50) 
with substitution x = yfay. Hence, focus on the so-called Gaussian integral 

[ dxe~ x \ (B.51) 

On the one hand, 




(B.52) 



On the other hand, 




(B.53) 



where dA denotes surface integration, which has been performed in polar coordinates. 
Furthermore, r — i/s has been substituted. With this result and equations (B.50) and 
(B.52), 

/ dxe~^ = y/arz. (B.54) 
Jr 

X 2 

• The integral J K dx xe~ **~ . 

X 2 X 2 

This is easy. While x is an odd function, e^^r is an even function. Hence, xe~^~ is odd 



B.8. Geometric Series 
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and therefore, 

/ dxxe~^ = 0. (B.55) 
Jr 



x 2 

The integral J" R dx x 2 e~ t . 



[ dx x 2 e " — f dxa 2 —e •> — a 2 —( dx e 
Jr Jr o« da Jr 

(B.56) 

2 o , — v a 6 n 



a ' —van 



da 2 



B.8. Geometric Series 



In section 4.4.2 some formulas regarding infinite geometric series have been exploited. These 
shall be proved here. 

• The series Y4L1 x ' with < x < 1. 



n n 



Hence, 



For n = 00, 



Finally, 



(1 - x) E x l = E x l - x E x* = 1 - x n+1 . (B.57) 

i=0 i=0 i'=0 



n -1 _ r n+l 

E*' = ^— . {B - 58) 



i=0 1 x 

1 . x 



CO CO 



L^ = L^-^ = — x -^ = — x - (B.60) 

; I / l! 



• The series J^Lj ix i with < x < 1. 



00 00 ^ 3 00 

E^' = E*^' = ^E* ! 

,=i ,-=1 ax ° x 1=1 

= x 



(B.61) 



dxl — x (1 — x) 2 ' 
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is taken from [vLLMN08] with permission by Peter van Loock 101 
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